/*


!
!  Dalton, a molecular electronic structure program
!  Copyright (C) The Dalton Authors (see AUTHORS file for details).
!
!  This program is free software; you can redistribute it and/or
!  modify it under the terms of the GNU Lesser General Public
!  License version 2.1 as published by the Free Software Foundation.
!
!  This program is distributed in the hope that it will be useful,
!  but WITHOUT ANY WARRANTY; without even the implied warranty of
!  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
!  Lesser General Public License for more details.
!
!  If a copy of the GNU LGPL v2.1 was not distributed with this
!  code, you can obtain one at https://www.gnu.org/licenses/old-licenses/lgpl-2.1.en.html.
!

!

*/
/*-*-mode: C; c-indentation-style: "bsd"; c-basic-offset: 4; -*-*/
/* fun-pw92ac.c:

   Reference: J. P. Perdew, J. A. Chevary, S. H. Vosko, K. A. Jackson, 
              M. R. Pederson and C. Fiolhais, Phys. Rev. B 46, 6671 (1992)

   Correction to PW91c functional due to neglect of \Delta \zeta. 
   This is equation 16 in the original reference, in which case the
   full correlation functional labeled pwGGA-IIA comprises:
   PW91c + PW92ac

   Implemented by Pawel Salek
   (c) Pawel Salek, pawsa@theochem.kth.se, sep 2001
   4th derivatives added by David Wilson (david.wilson@latrobe.edu.au) June 2005 

   Automatically generated code implementing PW92ac functional and
   its derivatives. It is generated by func-codegen.pl being a part of
   a "Automatic code generation framework for analytical functional
   derivative evaluation", Pawel Salek, 2005

    This functional is connected by making following changes:
    1. add "extern Functional pw92acFunctional;" to 'functionals.h'
    2. add "&pw92acFunctional," to 'functionals.c'
    3. add "fun-pw92ac.c" to 'Makefile.am', 'Makefile.in' or 'Makefile'.

    This functional has been generated from following input:
    ------ cut here -------
Cc0:0.004235;

rho:  rhoa + rhob;
grad: sqrt(grada*grada + gradb*gradb + 2*gradab);
zeta: (rhoa-rhob)/(rhoa+rhob);

gzgr: (1-zeta)/rho*grada^2 - (1+zeta)/rho*gradb^2 - (2*zeta/rho)*gradab;
t1: -0.458*zeta * gzgr/((rho*(1-zeta^2))^(1/3)*rho);
gz2: (1-zeta)^2/rho^2*grada^2 + (1+zeta)^2/rho^2*gradb^2 - 2*(1-zeta^2)/rho^2*gradab;
t2: ( (-0.037+0.10*zeta^2)*gz2 / rho^(1/3)*(1-zeta^2) );

K(rhoa,rhob,grada,gradb,gradab):=Cc0*rho*(t1+t2);

    ------ cut here -------
*/
 
/* strictly conform to XOPEN ANSI C standard */
#if !defined(SYS_DEC)
/* XOPEN compliance is missing on old Tru64 4.0E Alphas and pow() prototype
 * is not specified. */
#define _XOPEN_SOURCE          500
#define _XOPEN_SOURCE_EXTENDED 1
#endif
#include <math.h>
#include <stddef.h>
#include "general.h"

#define __CVERSION__
 
#include "functionals.h"
 
/* INTERFACE PART */
static integer pw92ac_isgga(void) { return 1; } /* FIXME: detect! */
static integer pw92ac_read(const char *conf_line);
static real pw92ac_energy(const FunDensProp* dp);
static void pw92ac_first(FunFirstFuncDrv *ds,   real factor,
                         const FunDensProp* dp);
static void pw92ac_second(FunSecondFuncDrv *ds, real factor,
                          const FunDensProp* dp);
static void pw92ac_third(FunThirdFuncDrv *ds,   real factor,
                         const FunDensProp* dp);
static void pw92ac_fourth(FunFourthFuncDrv *ds,   real factor,
                          const FunDensProp* dp);
 
Functional PW92acFunctional = {
  "PW92ac",       /* name */
  pw92ac_isgga,   /* gga-corrected */
   1,
  pw92ac_read,
  NULL,
  pw92ac_energy,
  pw92ac_first,
  pw92ac_second,
  pw92ac_third,
  pw92ac_fourth
};
 
/* IMPLEMENTATION PART */
static integer
pw92ac_read(const char *conf_line)
{
    fun_set_hf_weight(0);
    return 1;
}

static real
pw92ac_energy(const FunDensProp *dp)
{
    real res;
    real rhoa = dp->rhoa, rhob = dp->rhob;
    real grada = dp->grada, gradb = dp->gradb, gradab = dp->gradab;

    real t1, t2, t3, t4, t5, t6, t7, t8, t9, t10;

    t1 = rhob+rhoa;
    t2 = rhoa-1.0*rhob;
    t3 = pow(t2,2.0);
    t4 = 1/pow(t1,2.0);
    t5 = 1.0-1.0*t3*t4;
    t6 = pow(grada,2.0);
    t7 = 1/t1;
    t8 = 1.0-1.0*t2*t7;
    t9 = pow(gradb,2.0);
    t10 = t2*t7+1.0;

   /* code */
    res = 0.004235*t1*((0.1*t3*t4-0.037)*t5*(-2.0*t4*t5*gradab+
        pow(t10,2.0)*t4*t9+t4*t6*pow(t8,2.0))/pow(t1,0.333333333333333)-
        0.458*t2*(-2.0*t2*t4*gradab-1.0*t10*t7*t9+t6*t7*t8)/(pow(t1,
        2.333333333333333)*pow(t5,0.333333333333333)));

    return res;
}

static void
pw92ac_first(FunFirstFuncDrv *ds, real factor, const FunDensProp *dp)
{
    real dfdra, dfdrb, dfdga, dfdgb, dfdgab;
    real rhoa = dp->rhoa, rhob = dp->rhob;
    real grada = dp->grada, gradb = dp->gradb, gradab = dp->gradab;

    real t1, t2, t3, t4, t5, t6, t7, t8, t9, t10;
    real t11, t12, t13, t14, t15, t16, t17, t18;
    real t19, t20, t21, t22, t23, t24, t25, t26;
    real t27, t28, t29, t30, t31, t32, t33, t34;
    real t35, t36, t37, t38, t39, t40, t41;

    t1 = rhob+rhoa;
    t2 = rhoa-1.0*rhob;
    t3 = 1/pow(t1,2.333333333333333);
    t4 = pow(t2,2.0);
    t5 = 1/pow(t1,2.0);
    t6 = 1.0-1.0*t4*t5;
    t7 = 1/pow(t6,0.333333333333333);
    t8 = 1/pow(t1,3.0);
    t9 = 4.0*t2*t8*gradab;
    t10 = pow(grada,2.0);
    t11 = 1/t1;
    t12 = t2*t5;
    t13 = -1.0*t11;
    t14 = t13+t12;
    t15 = pow(gradb,2.0);
    t16 = -1.0*t2*t5;
    t17 = t11+t16;
    t18 = 1.0-1.0*t11*t2;
    t19 = -1.0*t10*t18*t5;
    t20 = t2*t11+1.0;
    t21 = t15*t5*t20;
    t22 = 2.0*t4*t8;
    t23 = t22-2.0*t2*t5;
    t24 = 1/pow(t6,1.333333333333333);
    t25 = -2.0*t2*t5*gradab-1.0*t11*t15*t20+t10*t11*t18;
    t26 = 1/
        pow(t1,3.333333333333333);
    t27 = 1.068666666666667*t2*t26*t7*t25;
    t28 = 1/pow(t1,0.333333333333333);
    t29 = 0.1*t4*t5-0.037;
    t30 = 4.0*t6*t8*gradab;
    t31 = pow(t18,2.0);
    t32 = -2.0*t10*t31*t8;
    t33 = pow(t20,2.0);
    t34 = -2.0*t15*t33*t8;
    t35 = -0.2*t4*t8;
    t36 = -2.0*t5*t6*gradab+t15*t5*t33+t10*t5*t31;
    t37 = -0.333333333333333*t29*t36*t6/pow(t1,1.333333333333333);
    t38 = 0.004235*
        (t28*t6*t29*t36-0.458*t2*t3*t7*t25);
    t39 = t13+t16;
    t40 = t11+t12;
    t41 = 2.0*t2*t5+t22;

   /* code */
    dfdra = 0.004235*t1*(t28*t29*t6*(-2.0*t23*t5*gradab+2.0*
        t15*t17*t20*t5+2.0*t10*t14*t18*t5+t34+t32+t30)-0.458*t2*t3*
        t7*(-2.0*t5*gradab+t9+t21+t19-1.0*t11*t15*t17+t10*t11*t14)+
        t37+t28*(0.2*t2*t5+t35)*t6*t36+t28*t23*t29*t36+t27-0.458*t3*
        t7*t25+0.152666666666667*t2*t3*t23*t24*t25)+t38;
    dfdrb = 0.004235*t1*(t28*t29*t6*(-2.0*t41*t5*gradab+2.0*
        t10*t18*t40*t5+2.0*t15*t20*t39*t5+t34+t32+t30)-0.458*t2*t3*
        t7*(2.0*t5*gradab+t9+t10*t11*t40-1.0*t11*t15*t39+t21+t19)+
        t37+t28*(t35-0.2*t2*t5)*t6*t36+t28*t41*t29*t36+t27+0.458*t3*
        t7*t25+0.152666666666667*t2*t3*t41*t24*t25)+t38;
    dfdga = 0.004235*t1*(2.0*t29*t3*t31*t6*grada-0.916*grada*
        t2*t26*t7*t18);
    dfdgb = 0.004235*t1*(2.0*t29*t3*t33*t6*gradb+0.916*gradb*
        t2*t26*t7*t20);
    dfdgab = 0.004235*t1*(0.916*t4*t7/pow(t1,4.333333333333333)-
        2.0*t29*t3*pow(t6,2.0));


    ds->df1000 += factor*dfdra;
    ds->df0100 += factor*dfdrb;
    ds->df0010 += factor*dfdga;
    ds->df0001 += factor*dfdgb;
    ds->df00001+= factor*dfdgab;
   
}

static void
pw92ac_second(FunSecondFuncDrv *ds, real factor, const FunDensProp* dp)
{
    real dfdra, dfdrb, dfdga, dfdgb, dfdgab;
    real d2fdrara, d2fdrarb, d2fdraga, d2fdragb, d2fdraab, d2fdrbrb,
        d2fdrbga, d2fdrbgb, d2fdrbgab, d2fdgaga, d2fdgagb, d2fdgagab,
        d2fdgbgb, d2fdgbgab, d2fdgabgab;
    real rhoa = dp->rhoa, rhob = dp->rhob;
    real grada = dp->grada, gradb = dp->gradb, gradab = dp->gradab;

    real t1, t2, t3, t4, t5, t6, t7, t8, t9, t10;
    real t11, t12, t13, t14, t15, t16, t17, t18;
    real t19, t20, t21, t22, t23, t24, t25, t26;
    real t27, t28, t29, t30, t31, t32, t33, t34;
    real t35, t36, t37, t38, t39, t40, t41, t42;
    real t43, t44, t45, t46, t47, t48, t49, t50;
    real t51, t52, t53, t54, t55, t56, t57, t58;
    real t59, t60, t61, t62, t63, t64, t65, t66;
    real t67, t68, t69, t70, t71, t72, t73, t74;
    real t75, t76, t77, t78, t79, t80, t81, t82;
    real t83, t84, t85, t86, t87, t88, t89, t90;
    real t91, t92, t93, t94, t95;

    t1 = rhob+rhoa;
    t2 = rhoa-1.0*rhob;
    t3 = 1/pow(t1,2.333333333333333);
    t4 = pow(t2,2.0);
    t5 = 1/pow(t1,2.0);
    t6 = 1.0-1.0*t4*t5;
    t7 = 1/pow(t6,0.333333333333333);
    t8 = 1/pow(t1,3.0);
    t9 = 4.0*t2*t8*gradab;
    t10 = pow(grada,2.0);
    t11 = 1/t1;
    t12 = t2*t5;
    t13 = -1.0*t11;
    t14 = t13+t12;
    t15 = pow(gradb,2.0);
    t16 = -1.0*t2*t5;
    t17 = t11+t16;
    t18 = 1.0-1.0*t11*t2;
    t19 = -1.0*t10*t18*t5;
    t20 = t2*t11+1.0;
    t21 = t15*t5*t20;
    t22 = -2.0*t5*gradab+t9+t21+t19-1.0*t11*t15*t17+t10*t11*
        t14;
    t23 = 2.0*t4*t8;
    t24 = t23-2.0*t2*t5;
    t25 = 1/pow(t6,1.333333333333333);
    t26 = -2.0*t2*t5*gradab-1.0*t11*t15*t20+t10*t11*t18;
    t27 = 1/
        pow(t1,3.333333333333333);
    t28 = 1.068666666666667*t2*t27*t7*t26;
    t29 = 1/pow(t1,0.333333333333333);
    t30 = 0.1*t4*t5-0.037;
    t31 = 4.0*t6*t8*gradab;
    t32 = pow(t18,2.0);
    t33 = -2.0*t10*t32*t8;
    t34 = pow(t20,2.0);
    t35 = -2.0*t15*t34*t8;
    t36 = -2.0*t24*t5*gradab+2.0*t15*t17*t20*t5+2.0*t10*t14*
        t18*t5+t35+t33+t31;
    t37 = -0.2*t4*t8;
    t38 = 0.2*t2*t5+t37;
    t39 = -2.0*t5*t6*gradab+t15*t5*t34+t10*t5*t32;
    t40 = 1/pow(t1,1.333333333333333);
    t41 = -0.333333333333333*t30*t39*t40*t6;
    t42 = t41+t29*t24*t30*t39+t29*t38*t6*t39+t29*t6*t30*t36-
        0.458*t3*t7*t26+t28+0.152666666666667*t2*t3*t24*t25*t26-0.458*
        t2*t3*t7*t22;
    t43 = 0.004235*(t29*t6*t30*t39-0.458*t2*t3*t7*t26);
    t44 = t13+
        t16;
    t45 = t11+t12;
    t46 = 2.0*t5*gradab+t9+t10*t11*t45-1.0*t11*t15*t44+t21+
        t19;
    t47 = 2.0*t2*t5+t23;
    t48 = -2.0*t47*t5*gradab+2.0*t10*t18*t45*t5+2.0*t15*t20*
        t44*t5+t35+t33+t31;
    t49 = t37-0.2*t2*t5;
    t50 = t41+t29*t47*t30*t39+t29*t49*t6*t39+t29*t6*t30*t48+
        0.458*t3*t7*t26+t28+0.152666666666667*t2*t3*t47*t25*t26-0.458*
        t2*t3*t7*t46;
    t51 = 2.0*t3*t30*t32*t6*grada-0.916*grada*t2*t27*t7*t18;
    t52 = 2.0*
        t3*t30*t34*t6*gradb+0.916*gradb*t2*t27*t7*t20;
    t53 = 1/pow(t1,4.333333333333333);
    t54 = pow(t6,2.0);
    t55 = 0.916*t4*t53*t7-2.0*t3*t30*t54;
    t56 = 1/pow(t1,4.0);
    t57 = -12.0*t2*t56*gradab;
    t58 = 2.0*t2*t8;
    t59 = -2.0*t5;
    t60 = t59+t58;
    t61 = -2.0*t2*t8;
    t62 = 2.0*t5;
    t63 = t62+t61;
    t64 = 2.0*t10*t18*t8;
    t65 = -2.0*t15*t20*t8;
    t66 = 1/pow(t6,2.333333333333333);
    t67 = -6.0*t4*t56;
    t68 = 8.0*t2*t8+t67+t59;
    t69 = -3.562222222222223*t2*t53*t7*t26;
    t70 = -12.0*t56*t6*gradab;
    t71 = 6.0*t10*t32*t56;
    t72 = 6.0*t15*t34*t56;
    t73 = 0.6*t4*t56;
    t74 = 0.2*t5;
    t75 = 0.444444444444444*t3*t30*t39*t6;
    t76 = t62+t67;
    t77 = 1/pow(t1,5.0);
    t78 = 0.004235*t51;
    t79 = 2.137333333333334*grada*t2*t53*t7*t18;
    t80 = -0.666666666666667*t27*t30*t32*t6*grada;
    t81 = -2.0*t18*t5*grada;
    t82 = -4.0*t32*t8*grada;
    t83 = 0.004235*t52;
    t84 = -2.137333333333334*gradb*t2*t53*t7*t20;
    t85 = -0.666666666666667*t27*t30*t34*t6*gradb;
    t86 = 2.0*t20*t5*gradb;
    t87 = -4.0*t34*t8*gradb;
    t88 = 0.004235*t55;
    t89 = -2.137333333333334*t4*t7/pow(t1,5.333333333333333);
    t90 = 4.0*
        t2*t8;
    t91 = 0.666666666666667*t27*t30*t54;
    t92 = 4.0*t6*t8;
    t93 = t59+t61;
    t94 = t62+t58;
    t95 = -8.0*t2*t8+t67+t59;

   /* code */
    dfdra = t43+0.004235*t1*t42;
    dfdrb = t43+0.004235*t1*t50;
    dfdga = 0.004235*t1*t51;
    dfdgb = 0.004235*t1*t52;
    dfdgab = 0.004235*t1*t55;
    d2fdrara = 0.004235*t1*(t29*t30*t6*(8.0*t24*t8*gradab-
        2.0*t5*t68*gradab-8.0*t15*t17*t20*t8-8.0*t10*t14*t18*t8+t72+
        t71+t70+2.0*t10*t18*t5*t63+2.0*t15*t20*t5*t60+2.0*t15*pow(t17,
        2.0)*t5+2.0*t10*pow(t14,2.0)*t5)-0.458*t2*t3*t7*(8.0*t8*gradab+
        t65+t64+t10*t11*t63-1.0*t11*t15*t60+t57+2.0*t15*t17*t5-2.0*
        t10*t14*t5)+t75+t69-0.203555555555556*t2*pow(t24,2.0)*t26*
        t3*t66-0.666666666666667*t38*t39*t40*t6-0.666666666666667*
        t30*t36*t40*t6+2.0*t29*t36*t38*t6-0.666666666666667*t24*t30*
        t39*t40+t29*(t74-0.8*t2*t8+t73)*t6*t39+2.0*t24*t29*t38*t39+
        t29*t68*t30*t39+2.0*t24*t29*t30*t36+2.137333333333334*t27*
        t7*t26+0.152666666666667*t2*t3*t68*t25*t26+0.305333333333333*
        t3*t24*t25*t26-0.712444444444445*t2*t27*t24*t25*t26-0.916*
        t3*t7*t22+2.137333333333334*t2*t27*t7*t22+0.305333333333333*
        t2*t3*t24*t25*t22)+0.00847*t42;
    d2fdrarb = 0.004235*t1*(t29*t30*t6*(4.0*t47*t8*gradab+
        4.0*t24*t8*gradab-2.0*t5*t76*gradab-4.0*t10*t18*t45*t8-4.0*
        t15*t20*t44*t8-4.0*t15*t17*t20*t8-4.0*t10*t14*t18*t8+4.0*t15*
        t2*t20*t77-4.0*t10*t18*t2*t77+t72+t71+t70+2.0*t10*t14*t45*
        t5+2.0*t15*t17*t44*t5)+t75-0.458*t2*t3*(t65+t64+t57-2.0*t15*
        t2*t56-2.0*t10*t2*t56-1.0*t10*t45*t5-1.0*t10*t14*t5+t15*t5*
        t44+t15*t5*t17)*t7+t69-0.333333333333333*t39*t40*t49*t6-0.333333333333333*
        t30*t40*t48*t6-0.333333333333333*t38*t39*t40*t6-0.333333333333333*
        t30*t36*t40*t6+t29*t38*t6*t48+t29*t24*t30*t48-0.333333333333333*
        t30*t39*t40*t47-0.458*t3*t7*t46+1.068666666666667*t2*t27*t7*
        t46+0.152666666666667*t2*t3*t24*t25*t46-0.333333333333333*
        t24*t30*t39*t40+t29*(t73-0.2*t5)*t6*t39+t29*t24*t49*t39+t29*
        t38*t47*t39+t29*t76*t30*t39+t29*t49*t6*t36+t29*t47*t30*t36-
        0.203555555555556*t2*t3*t24*t47*t66*t26+0.152666666666667*
        t2*t3*t76*t25*t26+0.152666666666667*t3*t47*t25*t26-0.356222222222222*
        t2*t27*t47*t25*t26-0.152666666666667*t3*t24*t25*t26-0.356222222222222*
        t2*t27*t24*t25*t26+0.458*t3*t7*t22+1.068666666666667*t2*t27*
        t7*t22+0.152666666666667*t2*t3*t47*t25*t22)+0.004235*t50+0.004235*
        t42;
    d2fdraga = 0.004235*t1*(t29*t30*t6*(4.0*t14*t18*t5*grada+
        t82)-0.458*t2*t3*t7*(2.0*t11*t14*grada+t81)+2.0*t3*t32*t38*
        t6*grada+2.0*t24*t3*t30*t32*grada+t80+t79-0.916*grada*t27*
        t7*t18+0.305333333333333*grada*t2*t27*t24*t25*t18)+t78;
    d2fdragb = 0.004235*t1*(t29*t30*t6*(4.0*t17*t20*t5*gradb+
        t87)-0.458*t2*t3*t7*(t86-2.0*t11*t17*gradb)+2.0*t3*t34*t38*
        t6*gradb+2.0*t24*t3*t30*t34*gradb+t85+t84+0.916*gradb*t27*
        t7*t20-0.305333333333333*gradb*t2*t27*t24*t25*t20)+t83;
    d2fdraab = 0.004235*t1*(t29*t30*t6*(t92-2.0*t24*t5)+t91+
        t89-0.458*t2*t3*(t59+t90)*t7+0.916*t2*t53*t7-2.0*t24*t3*t30*
        t6-2.0*t3*t38*t54-0.305333333333333*t4*t53*t24*t25)+t88;
    d2fdrbrb = 0.004235*
        t1*(t29*t30*t6*(-2.0*t5*t95*gradab+8.0*t47*t8*gradab+2.0*t15*
        t20*t5*t94+2.0*t10*t18*t5*t93-8.0*t10*t18*t45*t8-8.0*t15*t20*
        t44*t8+t72+t71+t70+2.0*t10*pow(t45,2.0)*t5+2.0*t15*pow(t44,
        2.0)*t5)-0.458*t2*t3*t7*(-8.0*t8*gradab-1.0*t11*t15*t94+t10*
        t11*t93+t65+t64+t57-2.0*t10*t45*t5+2.0*t15*t44*t5)+t75+t69-
        0.203555555555556*t2*t26*t3*pow(t47,2.0)*t66+2.0*t29*t48*t49*
        t6-0.666666666666667*t39*t40*t49*t6-0.666666666666667*t30*
        t40*t48*t6+2.0*t29*t39*t47*t49+2.0*t29*t30*t47*t48-0.666666666666667*
        t30*t39*t40*t47+0.916*t3*t7*t46+2.137333333333334*t2*t27*t7*
        t46+0.305333333333333*t2*t3*t47*t25*t46+t29*(t74+0.8*t2*t8+
        t73)*t6*t39+t29*t95*t30*t39-2.137333333333334*t27*t7*t26+0.152666666666667*
        t2*t3*t95*t25*t26-0.305333333333333*t3*t47*t25*t26-0.712444444444445*
        t2*t27*t47*t25*t26)+0.00847*t50;
    d2fdrbga = 0.004235*t1*(t29*t30*t6*(4.0*t18*t45*t5*grada+
        t82)-0.458*t2*t3*t7*(2.0*t11*t45*grada+t81)+2.0*t3*t32*t49*
        t6*grada+2.0*t3*t30*t32*t47*grada+t80+t79+0.916*grada*t27*
        t7*t18+0.305333333333333*grada*t2*t27*t47*t25*t18)+t78;
    d2fdrbgb = 0.004235*t1*(t29*t30*t6*(4.0*t20*t44*t5*gradb+
        t87)-0.458*t2*t3*t7*(t86-2.0*t11*t44*gradb)+2.0*t3*t34*t49*
        t6*gradb+2.0*t3*t30*t34*t47*gradb+t85+t84-0.916*gradb*t27*
        t7*t20-0.305333333333333*gradb*t2*t27*t47*t25*t20)+t83;
    d2fdrbgab = 0.004235*t1*(t29*t30*t6*(t92-2.0*t47*t5)+
        t91+t89-0.458*t2*t3*(t62+t90)*t7-0.916*t2*t53*t7-2.0*t3*t30*
        t47*t6-2.0*t3*t49*t54-0.305333333333333*t4*t53*t47*t25)+t88;
    d2fdgaga = 0.004235*
        t1*(2.0*t3*t30*t32*t6-0.916*t2*t27*t7*t18);
    d2fdgagb = 0.0;
    d2fdgagab = 0.0;
    d2fdgbgb = 0.004235*t1*(2.0*t3*t30*t34*t6+0.916*t2*t27*
        t7*t20);
    d2fdgbgab = 0.0;
    d2fdgabgab = 0.0;


    ds->df1000 += factor*dfdra;
    ds->df0100 += factor*dfdrb;
    ds->df0010 += factor*dfdga;
    ds->df0001 += factor*dfdgb;
    ds->df00001+= factor*dfdgab;

    ds->df2000 += factor*d2fdrara;
    ds->df1100 += factor*d2fdrarb;
    ds->df1010 += factor*d2fdraga;
    ds->df1001 += factor*d2fdragb;
    ds->df10001+= factor*d2fdraab;
    ds->df0200 += factor*d2fdrbrb;
    ds->df0110 += factor*d2fdrbga;
    ds->df0101 += factor*d2fdrbgb;
    ds->df01001+= factor*d2fdrbgab;
    ds->df0020 += factor*d2fdgaga;
    ds->df0011 += factor*d2fdgagb;
    ds->df00101+= factor*d2fdgagab;
    ds->df0002 += factor*d2fdgbgb;
    ds->df00011+= factor*d2fdgbgab;
    ds->df00002+= factor*d2fdgabgab;

}

static void
pw92ac_third(FunThirdFuncDrv *ds, real factor, const FunDensProp* dp)
{
    real dfdra, dfdrb, dfdga, dfdgb, dfdgab;
    real d2fdrara, d2fdrarb, d2fdraga, d2fdragb, d2fdraab, d2fdrbrb,
        d2fdrbga, d2fdrbgb, d2fdrbgab, d2fdgaga, d2fdgagb, d2fdgagab,
        d2fdgbgb, d2fdgbgab, d2fdgabgab;
    real d3fdrarara, d3fdrararb, d3fdraraga, d3fdraragb, d3fdraraab,
         d3fdrarbrb, d3fdrarbga, d3fdrarbgb, d3fdrarbab, d3fdragaga,
         d3fdragagb, d3fdragaab, d3fdragbgb, d3fdragbab, d3fdraabab,
         d3fdrbrbrb, d3fdrbrbga, d3fdrbrbgb, d3fdrbrbab, d3fdrbgaga,
         d3fdrbgagb, d3fdrbgaab, d3fdrbgbgb, d3fdrbgbab, d3fdrbabab,
         d3fdgagaga, d3fdgagagb, d3fdgagaab, d3fdgagbgb, d3fdgagbab,
         d3fdgaabab, d3fdgbgbgb, d3fdgbgbab, d3fdgbabab, d3fdababab;
    real rhoa = dp->rhoa, rhob = dp->rhob;
    real grada = dp->grada, gradb = dp->gradb, gradab = dp->gradab;

    real t1, t2, t3, t4, t5, t6, t7, t8, t9, t10;
    real t11, t12, t13, t14, t15, t16, t17, t18;
    real t19, t20, t21, t22, t23, t24, t25, t26;
    real t27, t28, t29, t30, t31, t32, t33, t34;
    real t35, t36, t37, t38, t39, t40, t41, t42;
    real t43, t44, t45, t46, t47, t48, t49, t50;
    real t51, t52, t53, t54, t55, t56, t57, t58;
    real t59, t60, t61, t62, t63, t64, t65, t66;
    real t67, t68, t69, t70, t71, t72, t73, t74;
    real t75, t76, t77, t78, t79, t80, t81, t82;
    real t83, t84, t85, t86, t87, t88, t89, t90;
    real t91, t92, t93, t94, t95, t96, t97, t98;
    real t99, t100, t101, t102, t103, t104, t105;
    real t106, t107, t108, t109, t110, t111, t112;
    real t113, t114, t115, t116, t117, t118, t119;
    real t120, t121, t122, t123, t124, t125, t126;
    real t127, t128, t129, t130, t131, t132, t133;
    real t134, t135, t136, t137, t138, t139, t140;
    real t141, t142, t143, t144, t145, t146, t147;
    real t148, t149, t150, t151, t152, t153, t154;
    real t155, t156, t157, t158, t159, t160, t161;
    real t162, t163, t164, t165, t166, t167, t168;
    real t169, t170, t171, t172, t173, t174, t175;
    real t176, t177, t178, t179, t180, t181, t182;
    real t183, t184, t185, t186, t187, t188, t189;
    real t190, t191, t192, t193, t194, t195, t196;
    real t197, t198, t199, t200, t201, t202, t203;
    real t204, t205;

    t1 = rhob+rhoa;
    t2 = rhoa-1.0*rhob;
    t3 = 1/pow(t1,2.333333333333333);
    t4 = pow(t2,2.0);
    t5 = 1/pow(t1,2.0);
    t6 = 1.0-1.0*t4*t5;
    t7 = 1/pow(t6,0.333333333333333);
    t8 = 1/pow(t1,3.0);
    t9 = 4.0*t2*t8*gradab;
    t10 = pow(grada,2.0);
    t11 = 1/t1;
    t12 = t2*t5;
    t13 = -1.0*t11;
    t14 = t13+t12;
    t15 = pow(gradb,2.0);
    t16 = -1.0*t2*t5;
    t17 = t11+t16;
    t18 = 1.0-1.0*t11*t2;
    t19 = -1.0*t10*t18*t5;
    t20 = t2*t11+1.0;
    t21 = t15*t5*t20;
    t22 = -2.0*t5*gradab+t9+t21+t19-1.0*t11*t15*t17+t10*t11*
        t14;
    t23 = 2.0*t4*t8;
    t24 = t23-2.0*t2*t5;
    t25 = 1/pow(t6,1.333333333333333);
    t26 = -2.0*t2*t5*gradab-1.0*t11*t15*t20+t10*t11*t18;
    t27 = 1/
        pow(t1,3.333333333333333);
    t28 = 1.068666666666667*t2*t27*t7*t26;
    t29 = 1/pow(t1,0.333333333333333);
    t30 = 0.1*t4*t5-0.037;
    t31 = 4.0*t6*t8*gradab;
    t32 = pow(t18,2.0);
    t33 = -2.0*t10*t32*t8;
    t34 = pow(t20,2.0);
    t35 = -2.0*t15*t34*t8;
    t36 = -2.0*t24*t5*gradab+2.0*t15*t17*t20*t5+2.0*t10*t14*
        t18*t5+t35+t33+t31;
    t37 = -0.2*t4*t8;
    t38 = 0.2*t2*t5+t37;
    t39 = -2.0*t5*t6*gradab+t15*t5*t34+t10*t5*t32;
    t40 = 1/pow(t1,1.333333333333333);
    t41 = -0.333333333333333*t30*t39*t40*t6;
    t42 = t41+t29*t24*t30*t39+t29*t38*t6*t39+t29*t6*t30*t36-
        0.458*t3*t7*t26+t28+0.152666666666667*t2*t3*t24*t25*t26-0.458*
        t2*t3*t7*t22;
    t43 = 0.004235*(t29*t6*t30*t39-0.458*t2*t3*t7*t26);
    t44 = t13+
        t16;
    t45 = t11+t12;
    t46 = 2.0*t5*gradab+t9+t10*t11*t45-1.0*t11*t15*t44+t21+
        t19;
    t47 = 2.0*t2*t5+t23;
    t48 = -2.0*t47*t5*gradab+2.0*t10*t18*t45*t5+2.0*t15*t20*
        t44*t5+t35+t33+t31;
    t49 = t37-0.2*t2*t5;
    t50 = t41+t29*t47*t30*t39+t29*t49*t6*t39+t29*t6*t30*t48+
        0.458*t3*t7*t26+t28+0.152666666666667*t2*t3*t47*t25*t26-0.458*
        t2*t3*t7*t46;
    t51 = 2.0*t3*t30*t32*t6*grada-0.916*grada*t2*t27*t7*t18;
    t52 = 2.0*
        t3*t30*t34*t6*gradb+0.916*gradb*t2*t27*t7*t20;
    t53 = 1/pow(t1,4.333333333333333);
    t54 = pow(t6,2.0);
    t55 = 0.916*t4*t53*t7-2.0*t3*t30*t54;
    t56 = 1/pow(t1,4.0);
    t57 = -12.0*t2*t56*gradab;
    t58 = 2.0*t2*t8;
    t59 = -2.0*t5;
    t60 = t59+t58;
    t61 = -2.0*t2*t8;
    t62 = 2.0*t5;
    t63 = t62+t61;
    t64 = 2.0*t10*t18*t8;
    t65 = -2.0*t15*t20*t8;
    t66 = 8.0*t8*gradab+t65+t64+t10*t11*t63-1.0*t11*t15*t60+
        t57+2.0*t15*t17*t5-2.0*t10*t14*t5;
    t67 = pow(t24,2.0);
    t68 = 1/pow(t6,2.333333333333333);
    t69 = -6.0*t4*t56;
    t70 = 8.0*t2*t8+t69+t59;
    t71 = -3.562222222222223*t2*t53*t7*t26;
    t72 = -12.0*t56*t6*gradab;
    t73 = pow(t14,2.0);
    t74 = pow(t17,2.0);
    t75 = 6.0*t10*t32*t56;
    t76 = 6.0*t15*t34*t56;
    t77 = 8.0*t24*t8*gradab-2.0*t5*t70*gradab-8.0*t15*t17*
        t20*t8-8.0*t10*t14*t18*t8+t76+t75+2.0*t15*t5*t74+2.0*t10*t5*
        t73+t72+2.0*t10*t18*t5*t63+2.0*t15*t20*t5*t60;
    t78 = 0.6*t4*t56;
    t79 = 0.2*t5;
    t80 = t79-0.8*t2*t8+t78;
    t81 = 0.444444444444444*t3*t30*t39*t6;
    t82 = t81+t29*t6*t30*t77+t71-0.458*t2*t3*t7*t66-0.666666666666667*
        t38*t39*t40*t6-0.666666666666667*t30*t36*t40*t6+2.0*t29*t36*
        t38*t6-0.666666666666667*t24*t30*t39*t40+t29*t80*t6*t39+2.0*
        t24*t29*t38*t39+t29*t70*t30*t39+2.0*t24*t29*t30*t36+2.137333333333334*
        t27*t7*t26-0.203555555555556*t2*t3*t67*t68*t26+0.152666666666667*
        t2*t3*t70*t25*t26+0.305333333333333*t3*t24*t25*t26-0.712444444444445*
        t2*t27*t24*t25*t26-0.916*t3*t7*t22+2.137333333333334*t2*t27*
        t7*t22+0.305333333333333*t2*t3*t24*t25*t22;
    t83 = t65+t64+t57-2.0*t15*t2*t56-2.0*t10*t2*t56-1.0*t10*
        t45*t5-1.0*t10*t14*t5+t15*t5*t44+t15*t5*t17;
    t84 = t62+t69;
    t85 = 1/pow(t1,5.0);
    t86 = -2.0*t5*t84*gradab+4.0*t47*t8*gradab+4.0*t24*t8*
        gradab+4.0*t15*t2*t20*t85-4.0*t10*t18*t2*t85-4.0*t10*t18*t45*
        t8-4.0*t15*t20*t44*t8-4.0*t15*t17*t20*t8-4.0*t10*t14*t18*t8+
        t76+t75+t72+2.0*t10*t14*t45*t5+2.0*t15*t17*t44*t5;
    t87 = t78-0.2*t5;
    t88 = t29*t6*t30*t86-0.458*t2*t3*t7*t83+t81+t71-0.333333333333333*
        t39*t40*t49*t6-0.333333333333333*t30*t40*t48*t6-0.333333333333333*
        t38*t39*t40*t6-0.333333333333333*t30*t36*t40*t6+t29*t38*t6*
        t48+t29*t24*t30*t48-0.333333333333333*t30*t39*t40*t47-0.458*
        t3*t7*t46+1.068666666666667*t2*t27*t7*t46+0.152666666666667*
        t2*t3*t24*t25*t46-0.333333333333333*t24*t30*t39*t40+t29*t87*
        t6*t39+t29*t24*t49*t39+t29*t38*t47*t39+t29*t84*t30*t39+t29*
        t49*t6*t36+t29*t47*t30*t36-0.203555555555556*t2*t3*t24*t47*
        t68*t26+0.152666666666667*t2*t3*t84*t25*t26+0.152666666666667*
        t3*t47*t25*t26-0.356222222222222*t2*t27*t47*t25*t26-0.152666666666667*
        t3*t24*t25*t26-0.356222222222222*t2*t27*t24*t25*t26+0.458*
        t3*t7*t22+1.068666666666667*t2*t27*t7*t22+0.152666666666667*
        t2*t3*t47*t25*t22;
    t89 = 0.004235*t51;
    t90 = 2.137333333333334*grada*t2*t53*t7*t18;
    t91 = -0.666666666666667*t27*t30*t32*t6*grada;
    t92 = -2.0*t18*t5*grada;
    t93 = 2.0*t11*t14*grada+t92;
    t94 = -4.0*t32*t8*grada;
    t95 = 4.0*t14*t18*t5*grada+t94;
    t96 = 2.0*t3*t32*t38*t6*grada+2.0*t24*t3*t30*t32*grada+
        t29*t6*t30*t95-0.458*t2*t3*t7*t93+t91+t90-0.916*grada*t27*
        t7*t18+0.305333333333333*grada*t2*t27*t24*t25*t18;
    t97 = 0.004235*t52;
    t98 = -2.137333333333334*gradb*t2*t53*t7*t20;
    t99 = -0.666666666666667*t27*t30*t34*t6*gradb;
    t100 = 2.0*t20*t5*gradb;
    t101 = t100-2.0*t11*t17*gradb;
    t102 = -4.0*t34*t8*gradb;
    t103 = 4.0*t17*t20*t5*gradb+t102;
    t104 = 2.0*t3*t34*t38*t6*gradb+2.0*t24*t3*t30*t34*gradb+
        t99+t98+0.916*gradb*t27*t7*t20-0.305333333333333*gradb*t2*
        t27*t24*t25*t20+t29*t6*t30*t103-0.458*t2*t3*t7*t101;
    t105 = 0.004235*t55;
    t106 = 1/pow(t1,5.333333333333333);
    t107 = -2.137333333333334*t4*t106*t7;
    t108 = 4.0*t2*t8;
    t109 = t59+t108;
    t110 = 0.666666666666667*t27*t30*t54;
    t111 = 4.0*t6*t8;
    t112 = t111-2.0*t24*t5;
    t113 = 0.916*t2*t53*t7-0.458*t2*t3*t109*t7-2.0*t24*t3*
        t30*t6-2.0*t3*t38*t54-0.305333333333333*t4*t53*t24*t25+t29*
        t6*t30*t112+t110+t107;
    t114 = t59+t61;
    t115 = t62+t58;
    t116 = -8.0*t8*gradab+t65+t64+t57-2.0*t10*t45*t5+2.0*
        t15*t44*t5-1.0*t11*t115*t15+t10*t11*t114;
    t117 = pow(t47,2.0);
    t118 = -8.0*t2*t8+t69+t59;
    t119 = pow(t44,2.0);
    t120 = pow(t45,2.0);
    t121 = 8.0*t47*t8*gradab-2.0*t118*t5*gradab-8.0*t10*t18*
        t45*t8-8.0*t15*t20*t44*t8+t76+t75+t72+2.0*t115*t15*t20*t5+
        2.0*t10*t114*t18*t5+2.0*t119*t15*t5+2.0*t10*t120*t5;
    t122 = t79+0.8*t2*t8+t78;
    t123 = t81+t71+2.0*t29*t48*t49*t6-0.666666666666667*t39*
        t40*t49*t6-0.666666666666667*t30*t40*t48*t6+2.0*t29*t39*t47*
        t49+2.0*t29*t30*t47*t48-0.666666666666667*t30*t39*t40*t47+
        0.916*t3*t7*t46+2.137333333333334*t2*t27*t7*t46+0.305333333333333*
        t2*t3*t47*t25*t46+t29*t122*t6*t39+t29*t118*t30*t39-2.137333333333334*
        t27*t7*t26-0.203555555555556*t2*t3*t117*t68*t26-0.305333333333333*
        t3*t47*t25*t26-0.712444444444445*t2*t27*t47*t25*t26+0.152666666666667*
        t2*t3*t118*t25*t26+t29*t6*t30*t121-0.458*t2*t3*t7*t116;
    t124 = 2.0*t11*t45*grada+t92;
    t125 = 4.0*t18*t45*t5*grada+t94;
    t126 = 2.0*t3*t32*t49*t6*grada+2.0*t3*t30*t32*t47*grada+
        t91+t90+0.916*grada*t27*t7*t18+0.305333333333333*grada*t2*
        t27*t47*t25*t18+t29*t6*t30*t125-0.458*t2*t3*t7*t124;
    t127 = t100-2.0*t11*t44*gradb;
    t128 = 4.0*t20*t44*t5*gradb+t102;
    t129 = 2.0*t3*t34*t49*t6*gradb+2.0*t3*t30*t34*t47*gradb+
        t99+t98-0.916*gradb*t27*t7*t20-0.305333333333333*gradb*t2*
        t27*t47*t25*t20+t29*t6*t30*t128-0.458*t2*t3*t7*t127;
    t130 = t62+t108;
    t131 = t111-2.0*t47*t5;
    t132 = -0.916*t2*t53*t7-0.458*t2*t3*t130*t7-2.0*t3*t30*
        t47*t6-2.0*t3*t49*t54-0.305333333333333*t4*t53*t47*t25+t29*
        t6*t30*t131+t110+t107;
    t133 = 2.0*t3*t30*t32*t6-0.916*t2*t27*t7*t18;
    t134 = 2.0*t3*t30*t34*t6+0.916*t2*t27*t7*t20;
    t135 = 48.0*t2*t85*gradab;
    t136 = 6.0*t2*t56;
    t137 = -6.0*t8;
    t138 = t137+t136;
    t139 = -6.0*t2*t56;
    t140 = 6.0*t8;
    t141 = t140+t139;
    t142 = -6.0*t10*t18*t56;
    t143 = 6.0*t15*t20*t56;
    t144 = 1/pow(t6,3.333333333333333);
    t145 = 24.0*t4*t85;
    t146 = 12.0*t8;
    t147 = -36.0*t2*t56+t146+t145;
    t148 = 15.4362962962963*t2*t106*t7*t26;
    t149 = 48.0*t6*t85*gradab;
    t150 = -24.0*t10*t32*t85;
    t151 = -24.0*t15*t34*t85;
    t152 = -2.4*t4*t85;
    t153 = -1.2*t8;
    t154 = -1.037037037037037*t27*t30*t39*t6;
    t155 = t136-2.0*t8;
    t156 = 2.0*t8+t139;
    t157 = 2.137333333333334*t2*t27*t7*t83;
    t158 = -0.712444444444445*t2*t27*t47*t25*t22;
    t159 = -0.712444444444445*t2*t27*t24*t25*t46;
    t160 = 0.949925925925926*t2*t27*t24*t47*t68*t26;
    t161 = -12.0*t2*t56;
    t162 = -4.0*t8;
    t163 = t162+t161+t145;
    t164 = -0.712444444444445*t2*t27*t84*t25*t26;
    t165 = 8.0*t8*t84*gradab;
    t166 = -8.0*t15*t17*t44*t8;
    t167 = -8.0*t10*t14*t45*t8;
    t168 = 1/pow(t1,6.0);
    t169 = -0.666666666666667*t30*t40*t6*t86;
    t170 = -0.666666666666667*t38*t40*t48*t6;
    t171 = -0.666666666666667*t24*t30*t40*t48;
    t172 = -0.666666666666667*t36*t40*t49*t6;
    t173 = -0.666666666666667*t30*t36*t40*t47;
    t174 = -0.666666666666667*t24*t39*t40*t49;
    t175 = -0.666666666666667*t38*t39*t40*t47;
    t176 = 0.4*t8;
    t177 = -0.666666666666667*t39*t40*t6*t87;
    t178 = -0.666666666666667*t30*t39*t40*t84;
    t179 = 0.00847*t88;
    t180 = -7.124444444444446*grada*t2*t106*t7*t18;
    t181 = 0.888888888888889*t30*t32*t53*t6*grada;
    t182 = 4.0*t18*t8*grada;
    t183 = 12.0*t32*t56*grada;
    t184 = 7.124444444444446*gradb*t2*t106*t7*t20;
    t185 = 0.888888888888889*t30*t34*t53*t6*gradb;
    t186 = -4.0*t20*t8*gradb;
    t187 = 12.0*t34*t56*gradb;
    t188 = 1/pow(t1,6.333333333333333);
    t189 = 7.124444444444446*t4*t188*t7;
    t190 = -0.888888888888889*t30*t53*t54;
    t191 = -12.0*t56*t6;
    t192 = 12.0*t2*t56+t162+t145;
    t193 = 0.004235*t133;
    t194 = 2.137333333333334*t2*t53*t7*t18;
    t195 = -0.666666666666667*t27*t30*t32*t6;
    t196 = -2.0*t18*t5;
    t197 = -4.0*t32*t8;
    t198 = 0.004235*t134;
    t199 = -2.137333333333334*t2*t53*t7*t20;
    t200 = -0.666666666666667*t27*t30*t34*t6;
    t201 = 2.0*t20*t5;
    t202 = -4.0*t34*t8;
    t203 = t137+t139;
    t204 = t140+t136;
    t205 = 36.0*t2*t56+t146+t145;

   /* code */
    dfdra = t43+0.004235*t1*t42;
    dfdrb = t43+0.004235*t1*t50;
    dfdga = 0.004235*t1*t51;
    dfdgb = 0.004235*t1*t52;
    dfdgab = 0.004235*t1*t55;
    d2fdrara = 0.00847*t42+0.004235*t1*t82;
    d2fdrarb = 0.004235*t50+0.004235*t42+0.004235*t1*t88;
    d2fdraga = 0.004235*
        t1*t96+t89;
    d2fdragb = 0.004235*t1*t104+t97;
    d2fdraab = 0.004235*t1*t113+t105;
    d2fdrbrb = 0.00847*t50+0.004235*t1*t123;
    d2fdrbga = 0.004235*t1*t126+t89;
    d2fdrbgb = 0.004235*t1*t129+t97;
    d2fdrbgab = 0.004235*t1*t132+t105;
    d2fdgaga = 0.004235*t1*t133;
    d2fdgagb = 0.0;
    d2fdgagab = 0.0;
    d2fdgbgb = 0.004235*t1*t134;
    d2fdgbgab = 0.0;
    d2fdgabgab = 0.0;
    d3fdrarara = 0.004235*t1*(t29*t30*t6*(12.0*t70*t8*gradab-
        36.0*t24*t56*gradab-2.0*t147*t5*gradab-12.0*t15*t74*t8-12.0*
        t10*t73*t8-12.0*t10*t18*t63*t8-12.0*t15*t20*t60*t8+6.0*t10*
        t14*t5*t63+6.0*t15*t17*t5*t60+36.0*t15*t17*t20*t56+36.0*t10*
        t14*t18*t56+2.0*t141*t15*t20*t5+2.0*t10*t138*t18*t5+t151+t150+
        t149)-0.458*t2*t3*t7*(-36.0*t56*gradab-6.0*t15*t17*t8+6.0*
        t10*t14*t8-3.0*t10*t5*t63+3.0*t15*t5*t60-1.0*t11*t141*t15+
        t143+t142+t10*t11*t138+t135)-1.0*t39*t40*t6*t80+3.0*t29*t36*
        t6*t80+3.0*t24*t29*t39*t80-1.0*t30*t40*t6*t77+3.0*t29*t38*
        t6*t77+3.0*t24*t29*t30*t77-1.0*t30*t39*t40*t70+3.0*t29*t38*
        t39*t70+3.0*t29*t30*t36*t70-1.374*t3*t7*t66+3.206*t2*t27*t7*
        t66+0.458*t2*t3*t24*t25*t66-2.0*t36*t38*t40*t6+1.333333333333333*
        t3*t38*t39*t6+1.333333333333333*t3*t30*t36*t6-2.0*t24*t38*
        t39*t40-2.0*t24*t30*t36*t40+t29*(t153+3.6*t2*t56+t152)*t6*
        t39+1.333333333333333*t24*t3*t30*t39+t29*t147*t30*t39+6.0*
        t24*t29*t36*t38+0.474962962962963*t144*t2*pow(t24,3.0)*t26*
        t3-10.68666666666667*t53*t7*t26-0.610666666666667*t3*t67*t68*
        t26+1.424888888888889*t2*t27*t67*t68*t26-0.610666666666667*
        t2*t3*t70*t24*t68*t26+0.458*t3*t70*t25*t26-1.068666666666667*
        t2*t27*t70*t25*t26+3.562222222222223*t2*t53*t24*t25*t26-2.137333333333334*
        t27*t24*t25*t26+0.152666666666667*t2*t3*t147*t25*t26-10.68666666666667*
        t2*t53*t7*t22+6.412*t27*t7*t22-0.610666666666667*t2*t3*t67*
        t68*t22+0.458*t2*t3*t70*t25*t22+0.916*t3*t24*t25*t22-2.137333333333334*
        t2*t27*t24*t25*t22+t154+t148)+0.012705*t82;
    d3fdrararb = 0.004235*t1*(t29*t30*t6*(4.0*t70*t8*gradab-
        12.0*t47*t56*gradab-24.0*t24*t56*gradab-2.0*t163*t5*gradab+
        8.0*t15*t17*t2*t85-8.0*t10*t14*t2*t85-4.0*t15*t74*t8-4.0*t10*
        t73*t8-4.0*t10*t18*t63*t8-4.0*t15*t20*t60*t8+2.0*t10*t45*t5*
        t63+2.0*t15*t44*t5*t60+12.0*t10*t18*t45*t56+12.0*t15*t20*t44*
        t56+24.0*t15*t17*t20*t56+24.0*t10*t14*t18*t56+2.0*t15*t156*
        t20*t5+2.0*t10*t155*t18*t5-16.0*t15*t168*t2*t20+16.0*t10*t168*
        t18*t2+t167+t166+t165+t151+t150+t149)-0.458*t2*t3*t7*(-12.0*
        t56*gradab+4.0*t15*t2*t85+4.0*t10*t2*t85+2.0*t10*t45*t8-2.0*
        t15*t44*t8-4.0*t15*t17*t8+4.0*t10*t14*t8-1.0*t10*t5*t63+t15*
        t5*t60-1.0*t11*t15*t156+t10*t11*t155+t143+t142+t135)+2.0*t29*
        t36*t6*t87+2.0*t24*t29*t39*t87+2.0*t29*t38*t6*t86+2.0*t24*
        t29*t30*t86+2.0*t29*t38*t39*t84+2.0*t29*t30*t36*t84-0.916*
        t3*t7*t83+0.305333333333333*t2*t3*t24*t25*t83-0.333333333333333*
        t39*t40*t6*t80+t29*t49*t6*t77-0.333333333333333*t30*t40*t6*
        t77+t29*t47*t30*t77-0.333333333333333*t30*t39*t40*t70+0.458*
        t3*t7*t66+1.068666666666667*t2*t27*t7*t66+0.152666666666667*
        t2*t3*t47*t25*t66+0.444444444444444*t3*t39*t49*t6+0.444444444444444*
        t3*t30*t48*t6-0.666666666666667*t36*t38*t40*t6+0.888888888888889*
        t3*t38*t39*t6+0.888888888888889*t3*t30*t36*t6+2.0*t24*t29*
        t36*t49+t29*t80*t6*t48+2.0*t24*t29*t38*t48+t29*t70*t30*t48+
        0.444444444444444*t3*t30*t39*t47+2.0*t29*t36*t38*t47-3.562222222222223*
        t2*t53*t7*t46+2.137333333333334*t27*t7*t46-0.203555555555556*
        t2*t3*t67*t68*t46+0.152666666666667*t2*t3*t70*t25*t46+0.305333333333333*
        t3*t24*t25*t46-0.666666666666667*t24*t38*t39*t40-0.666666666666667*
        t24*t30*t36*t40+t29*(t176+1.2*t2*t56+t152)*t6*t39+t29*t70*
        t49*t39+t29*t80*t47*t39+0.888888888888889*t24*t3*t30*t39+t29*
        t163*t30*t39-3.562222222222223*t53*t7*t26+0.203555555555556*
        t3*t67*t68*t26+0.474962962962963*t2*t27*t67*t68*t26-0.203555555555556*
        t2*t3*t70*t47*t68*t26-0.407111111111111*t3*t24*t47*t68*t26-
        0.407111111111111*t2*t3*t84*t24*t68*t26+0.305333333333333*
        t3*t84*t25*t26-0.152666666666667*t3*t70*t25*t26-0.356222222222222*
        t2*t27*t70*t25*t26+1.187407407407408*t2*t53*t47*t25*t26-0.712444444444445*
        t27*t47*t25*t26+2.374814814814815*t2*t53*t24*t25*t26+0.152666666666667*
        t2*t3*t163*t25*t26+0.474962962962963*t2*t3*t67*t47*t144*t26-
        7.124444444444446*t2*t53*t7*t22-0.407111111111111*t2*t3*t24*
        t47*t68*t22+0.305333333333333*t2*t3*t84*t25*t22+0.305333333333333*
        t3*t47*t25*t22-0.305333333333333*t3*t24*t25*t22-0.712444444444445*
        t2*t27*t24*t25*t22+t178+t177+t175+t174+t173+t172+t171+t170+
        t169+t164+t160+t159+t158+t157+t154+t148)+0.004235*t82+t179;
    d3fdraraga = 0.004235*
        t1*(t29*t30*t6*(-16.0*t14*t18*t8*grada+4.0*t5*t73*grada+4.0*
        t18*t5*t63*grada+t183)-0.458*t2*t3*t7*(2.0*t11*t63*grada-4.0*
        t14*t5*grada+t182)+2.0*t3*t32*t6*t80*grada+2.0*t3*t30*t32*
        t70*grada-1.333333333333333*t27*t32*t38*t6*grada+4.0*t24*t3*
        t32*t38*grada-1.333333333333333*t24*t27*t30*t32*grada-0.666666666666667*
        t30*t40*t6*t95+2.0*t29*t38*t6*t95+2.0*t24*t29*t30*t95-0.916*
        t3*t7*t93+2.137333333333334*t2*t27*t7*t93+0.305333333333333*
        t2*t3*t24*t25*t93+t181+t180+4.274666666666668*grada*t53*t7*
        t18-0.407111111111111*grada*t2*t27*t67*t68*t18+0.305333333333333*
        grada*t2*t27*t70*t25*t18-1.424888888888889*grada*t2*t53*t24*
        t25*t18+0.610666666666667*grada*t27*t24*t25*t18)+0.00847*t96;
    d3fdraragb = 0.004235*
        t1*(t29*t30*t6*(-16.0*t17*t20*t8*gradb+4.0*t5*t74*gradb+4.0*
        t20*t5*t60*gradb+t187)-0.458*t2*t3*t7*(-2.0*t11*t60*gradb+
        4.0*t17*t5*gradb+t186)+2.0*t3*t34*t6*t80*gradb+2.0*t3*t30*
        t34*t70*gradb-1.333333333333333*t27*t34*t38*t6*gradb+4.0*t24*
        t3*t34*t38*gradb-1.333333333333333*t24*t27*t30*t34*gradb-0.666666666666667*
        t103*t30*t40*t6+2.0*t103*t29*t38*t6+2.0*t103*t24*t29*t30-4.274666666666668*
        gradb*t53*t7*t20+0.407111111111111*gradb*t2*t27*t67*t68*t20-
        0.305333333333333*gradb*t2*t27*t70*t25*t20+1.424888888888889*
        gradb*t2*t53*t24*t25*t20-0.610666666666667*gradb*t27*t24*t25*
        t20+t185+t184-0.916*t3*t7*t101+2.137333333333334*t2*t27*t7*
        t101+0.305333333333333*t2*t3*t24*t25*t101)+0.00847*t104;
    d3fdraraab = 0.004235*
        t1*(-2.0*t3*t54*t80+t29*t30*t6*(8.0*t24*t8-2.0*t5*t70+t191)-
        0.458*t2*t3*t7*(8.0*t8+t161)-2.0*t3*t30*t6*t70-0.916*t3*t109*
        t7+2.137333333333334*t2*t27*t109*t7-4.274666666666668*t2*t106*
        t7+0.407111111111111*t4*t53*t67*t68-0.666666666666667*t112*
        t30*t40*t6-4.0*t24*t3*t38*t6+2.0*t112*t29*t38*t6+1.333333333333333*
        t24*t27*t30*t6+1.333333333333333*t27*t38*t54+2.0*t112*t24*
        t29*t30-0.305333333333333*t4*t53*t70*t25-0.610666666666667*
        t2*t53*t24*t25+0.305333333333333*t2*t3*t109*t24*t25+1.424888888888889*
        t4*t106*t24*t25+t190+t189)+0.00847*t113;
    d3fdrarbrb = 0.004235*t1*(t29*t30*t6*(4.0*t118*t8*gradab-
        24.0*t47*t56*gradab-12.0*t24*t56*gradab-2.0*t192*t5*gradab-
        8.0*t10*t2*t45*t85+8.0*t15*t2*t44*t85-4.0*t15*t20*t85+4.0*
        t10*t18*t85-4.0*t115*t15*t20*t8-4.0*t10*t114*t18*t8-4.0*t119*
        t15*t8-4.0*t10*t120*t8+24.0*t10*t18*t45*t56+24.0*t15*t20*t44*
        t56+12.0*t15*t17*t20*t56+12.0*t10*t14*t18*t56+2.0*t115*t15*
        t17*t5+2.0*t10*t114*t14*t5-28.0*t15*t168*t2*t20+28.0*t10*t168*
        t18*t2+t167+t166+t165+t151+t150+t149)-0.458*t2*t3*t7*(12.0*
        t56*gradab+10.0*t15*t2*t85+10.0*t10*t2*t85+4.0*t10*t45*t8-
        4.0*t15*t44*t8-2.0*t15*t17*t8+2.0*t10*t14*t8+2.0*t15*t56+2.0*
        t10*t56-1.0*t10*t114*t5+t143+t142+t135+t15*t5*t115)+2.0*t29*
        t48*t6*t87+2.0*t29*t39*t47*t87+2.0*t29*t49*t6*t86+2.0*t29*
        t30*t47*t86+2.0*t29*t39*t49*t84+2.0*t29*t30*t48*t84+0.916*
        t3*t7*t83+0.305333333333333*t2*t3*t47*t25*t83-0.666666666666667*
        t40*t48*t49*t6+0.888888888888889*t3*t39*t49*t6+0.888888888888889*
        t3*t30*t48*t6-0.333333333333333*t122*t39*t40*t6-0.333333333333333*
        t121*t30*t40*t6+0.444444444444444*t3*t38*t39*t6+0.444444444444444*
        t3*t30*t36*t6+2.0*t24*t29*t48*t49-0.666666666666667*t39*t40*
        t47*t49+2.0*t29*t36*t47*t49-0.666666666666667*t30*t40*t47*
        t48+2.0*t29*t38*t47*t48+0.888888888888889*t3*t30*t39*t47-7.124444444444446*
        t2*t53*t7*t46-0.407111111111111*t2*t3*t24*t47*t68*t46+0.305333333333333*
        t2*t3*t84*t25*t46+0.305333333333333*t3*t47*t25*t46-0.712444444444445*
        t2*t27*t47*t25*t46-0.305333333333333*t3*t24*t25*t46-0.333333333333333*
        t118*t30*t39*t40+t29*(t176-1.2*t2*t56+t152)*t6*t39+t29*t118*
        t38*t39+0.444444444444444*t24*t3*t30*t39+t29*t192*t30*t39+
        t29*t122*t24*t39+t29*t122*t6*t36+t29*t118*t30*t36+3.562222222222223*
        t53*t7*t26-0.407111111111111*t2*t3*t84*t47*t68*t26+0.407111111111111*
        t3*t24*t47*t68*t26-0.203555555555556*t2*t3*t118*t24*t68*t26-
        0.203555555555556*t3*t117*t68*t26+0.474962962962963*t2*t27*
        t117*t68*t26-0.305333333333333*t3*t84*t25*t26+2.374814814814815*
        t2*t53*t47*t25*t26+1.187407407407408*t2*t53*t24*t25*t26+0.712444444444445*
        t27*t24*t25*t26+0.152666666666667*t2*t3*t192*t25*t26+0.152666666666667*
        t3*t118*t25*t26-0.356222222222222*t2*t27*t118*t25*t26+0.474962962962963*
        t2*t3*t24*t117*t144*t26-3.562222222222223*t2*t53*t7*t22-2.137333333333334*
        t27*t7*t22-0.203555555555556*t2*t3*t117*t68*t22-0.305333333333333*
        t3*t47*t25*t22+0.152666666666667*t2*t3*t118*t25*t22+t178+t177+
        t175+t174+t173+t172+t171+t170+t169+t164+t160+t159+t158+t157+
        t154+t148+t29*t38*t6*t121+t29*t24*t30*t121-0.458*t3*t7*t116+
        1.068666666666667*t2*t27*t7*t116+0.152666666666667*t2*t3*t24*
        t25*t116)+t179+0.004235*t123;
    d3fdrarbga = 0.004235*t1*(t29*t30*t6*(-8.0*t18*t2*t85*
        grada-8.0*t18*t45*t8*grada-8.0*t14*t18*t8*grada+4.0*t14*t45*
        t5*grada+t183)-0.458*t2*t3*t7*(-4.0*t2*t56*grada-2.0*t45*t5*
        grada-2.0*t14*t5*grada+t182)+2.0*t3*t32*t6*t87*grada+2.0*t3*
        t30*t32*t84*grada-0.666666666666667*t27*t32*t49*t6*grada-0.666666666666667*
        t27*t32*t38*t6*grada+2.0*t24*t3*t32*t49*grada+2.0*t3*t32*t38*
        t47*grada-0.666666666666667*t27*t30*t32*t47*grada-0.666666666666667*
        t24*t27*t30*t32*grada+t29*t49*t6*t95-0.333333333333333*t30*
        t40*t6*t95+t29*t47*t30*t95+0.458*t3*t7*t93+1.068666666666667*
        t2*t27*t7*t93+0.152666666666667*t2*t3*t47*t25*t93-0.333333333333333*
        t125*t30*t40*t6+t181+t180-0.407111111111111*grada*t2*t27*t24*
        t47*t68*t18+0.305333333333333*grada*t2*t27*t84*t25*t18-0.712444444444445*
        grada*t2*t53*t47*t25*t18+0.305333333333333*grada*t27*t47*t25*
        t18-0.712444444444445*grada*t2*t53*t24*t25*t18-0.305333333333333*
        grada*t27*t24*t25*t18+t29*t38*t6*t125+t29*t24*t30*t125-0.458*
        t3*t7*t124+1.068666666666667*t2*t27*t7*t124+0.152666666666667*
        t2*t3*t24*t25*t124)+0.004235*t96+0.004235*t126;
    d3fdrarbgb = 0.004235*t1*(t29*t30*t6*(8.0*t2*t20*t85*
        gradb-8.0*t20*t44*t8*gradb-8.0*t17*t20*t8*gradb+4.0*t17*t44*
        t5*gradb+t187)-0.458*t2*t3*t7*(-4.0*t2*t56*gradb+2.0*t44*t5*
        gradb+2.0*t17*t5*gradb+t186)+2.0*t3*t34*t6*t87*gradb+2.0*t3*
        t30*t34*t84*gradb-0.666666666666667*t27*t34*t49*t6*gradb-0.666666666666667*
        t27*t34*t38*t6*gradb+2.0*t24*t3*t34*t49*gradb+2.0*t3*t34*t38*
        t47*gradb-0.666666666666667*t27*t30*t34*t47*gradb-0.666666666666667*
        t24*t27*t30*t34*gradb-0.333333333333333*t128*t30*t40*t6-0.333333333333333*
        t103*t30*t40*t6+0.407111111111111*gradb*t2*t27*t24*t47*t68*
        t20-0.305333333333333*gradb*t2*t27*t84*t25*t20+0.712444444444445*
        gradb*t2*t53*t47*t25*t20-0.305333333333333*gradb*t27*t47*t25*
        t20+0.712444444444445*gradb*t2*t53*t24*t25*t20+0.305333333333333*
        gradb*t27*t24*t25*t20+t185+t184+t29*t38*t6*t128+t29*t24*t30*
        t128-0.458*t3*t7*t127+1.068666666666667*t2*t27*t7*t127+0.152666666666667*
        t2*t3*t24*t25*t127+t29*t49*t6*t103+t29*t47*t30*t103+0.458*
        t3*t7*t101+1.068666666666667*t2*t27*t7*t101+0.152666666666667*
        t2*t3*t47*t25*t101)+0.004235*t129+0.004235*t104;
    d3fdrarbab = 0.004235*t1*(-2.0*t3*t54*t87+t29*t30*t6*
        (-2.0*t5*t84+4.0*t47*t8+4.0*t24*t8+t191)-2.0*t3*t30*t6*t84+
        12.62044444444445*t4*t188*t7-0.458*t3*t130*t7+1.068666666666667*
        t2*t27*t130*t7+0.458*t3*t109*t7+1.068666666666667*t2*t27*t109*
        t7+0.407111111111111*t4*t53*t24*t47*t68-2.0*t24*t3*t49*t6-
        2.0*t3*t38*t47*t6+0.666666666666667*t27*t30*t47*t6-0.333333333333333*
        t131*t30*t40*t6-0.333333333333333*t112*t30*t40*t6+0.666666666666667*
        t24*t27*t30*t6+0.666666666666667*t27*t49*t54+0.666666666666667*
        t27*t38*t54-0.305333333333333*t4*t53*t84*t25-0.305333333333333*
        t2*t53*t47*t25+0.152666666666667*t2*t3*t109*t47*t25+0.712444444444445*
        t4*t106*t47*t25+0.305333333333333*t2*t53*t24*t25+0.152666666666667*
        t2*t3*t130*t24*t25+0.712444444444445*t4*t106*t24*t25+t190+
        t29*t38*t6*t131+t29*t24*t30*t131+t29*t49*t6*t112+t29*t47*t30*
        t112)+0.004235*t132+0.004235*t113;
    d3fdragaga = 0.004235*t1*(-0.458*(t196+2.0*t11*t14)*t2*
        t3*t7+t29*t30*(4.0*t14*t18*t5+t197)*t6+2.0*t3*t32*t38*t6+2.0*
        t24*t3*t30*t32+t195+t194-0.916*t27*t7*t18+0.305333333333333*
        t2*t27*t24*t25*t18)+t193;
    d3fdragagb = 0.0;
    d3fdragaab = 0.0;
    d3fdragbgb = 0.004235*t1*(-0.458*t2*(t201-2.0*t11*t17)*
        t3*t7+t29*t30*(4.0*t17*t20*t5+t202)*t6+2.0*t3*t34*t38*t6+2.0*
        t24*t3*t30*t34+t200+0.916*t27*t7*t20-0.305333333333333*t2*
        t27*t24*t25*t20+t199)+t198;
    d3fdragbab = 0.0;
    d3fdraabab = 0.0;
    d3fdrbrbrb = 0.004235*t1*(t29*t30*t6*(12.0*t118*t8*gradab-
        36.0*t47*t56*gradab-2.0*t205*t5*gradab-12.0*t115*t15*t20*t8-
        12.0*t10*t114*t18*t8-12.0*t119*t15*t8-12.0*t10*t120*t8+36.0*
        t10*t18*t45*t56+36.0*t15*t20*t44*t56+6.0*t10*t114*t45*t5+6.0*
        t115*t15*t44*t5+2.0*t10*t18*t204*t5+2.0*t15*t20*t203*t5+t151+
        t150+t149)-0.458*t2*t3*t7*(36.0*t56*gradab+6.0*t10*t45*t8-
        6.0*t15*t44*t8+3.0*t115*t15*t5-3.0*t10*t114*t5+t10*t11*t204-
        1.0*t11*t15*t203+t143+t142+t135)-2.0*t40*t48*t49*t6+1.333333333333333*
        t3*t39*t49*t6+3.0*t121*t29*t49*t6+1.333333333333333*t3*t30*
        t48*t6+3.0*t122*t29*t48*t6-1.0*t122*t39*t40*t6-1.0*t121*t30*
        t40*t6+6.0*t29*t47*t48*t49-2.0*t39*t40*t47*t49+3.0*t118*t29*
        t39*t49-2.0*t30*t40*t47*t48+3.0*t118*t29*t30*t48+0.474962962962963*
        t144*t2*t26*t3*pow(t47,3.0)+1.333333333333333*t3*t30*t39*t47+
        3.0*t122*t29*t39*t47+3.0*t121*t29*t30*t47-10.68666666666667*
        t2*t53*t7*t46-6.412*t27*t7*t46-0.610666666666667*t2*t3*t117*
        t68*t46-0.916*t3*t47*t25*t46-2.137333333333334*t2*t27*t47*
        t25*t46+0.458*t2*t3*t118*t25*t46-1.0*t118*t30*t39*t40+t29*
        (t153-3.6*t2*t56+t152)*t6*t39+t29*t205*t30*t39+10.68666666666667*
        t53*t7*t26-0.610666666666667*t2*t3*t118*t47*t68*t26+0.610666666666667*
        t3*t117*t68*t26+1.424888888888889*t2*t27*t117*t68*t26+3.562222222222223*
        t2*t53*t47*t25*t26+2.137333333333334*t27*t47*t25*t26+0.152666666666667*
        t2*t3*t205*t25*t26-0.458*t3*t118*t25*t26-1.068666666666667*
        t2*t27*t118*t25*t26+t154+t148+1.374*t3*t7*t116+3.206*t2*t27*
        t7*t116+0.458*t2*t3*t47*t25*t116)+0.012705*t123;
    d3fdrbrbga = 0.004235*t1*(t29*t30*t6*(-16.0*t18*t45*t8*
        grada+4.0*t114*t18*t5*grada+4.0*t120*t5*grada+t183)-0.458*
        t2*t3*t7*(-4.0*t45*t5*grada+2.0*t11*t114*grada+t182)-1.333333333333333*
        t27*t32*t49*t6*grada+2.0*t122*t3*t32*t6*grada+4.0*t3*t32*t47*
        t49*grada-1.333333333333333*t27*t30*t32*t47*grada+2.0*t118*
        t3*t30*t32*grada+2.0*t125*t29*t49*t6-0.666666666666667*t125*
        t30*t40*t6+2.0*t125*t29*t30*t47+t181+t180-4.274666666666668*
        grada*t53*t7*t18-0.407111111111111*grada*t2*t27*t117*t68*t18-
        1.424888888888889*grada*t2*t53*t47*t25*t18-0.610666666666667*
        grada*t27*t47*t25*t18+0.305333333333333*grada*t2*t27*t118*
        t25*t18+0.916*t3*t7*t124+2.137333333333334*t2*t27*t7*t124+
        0.305333333333333*t2*t3*t47*t25*t124)+0.00847*t126;
    d3fdrbrbgb = 0.004235*t1*(t29*t30*t6*(-16.0*t20*t44*t8*
        gradb+4.0*t115*t20*t5*gradb+4.0*t119*t5*gradb+t187)-0.458*
        t2*t3*t7*(4.0*t44*t5*gradb-2.0*t11*t115*gradb+t186)-1.333333333333333*
        t27*t34*t49*t6*gradb+2.0*t122*t3*t34*t6*gradb+4.0*t3*t34*t47*
        t49*gradb-1.333333333333333*t27*t30*t34*t47*gradb+2.0*t118*
        t3*t30*t34*gradb+2.0*t128*t29*t49*t6-0.666666666666667*t128*
        t30*t40*t6+2.0*t128*t29*t30*t47+4.274666666666668*gradb*t53*
        t7*t20+0.407111111111111*gradb*t2*t27*t117*t68*t20+1.424888888888889*
        gradb*t2*t53*t47*t25*t20+0.610666666666667*gradb*t27*t47*t25*
        t20-0.305333333333333*gradb*t2*t27*t118*t25*t20+t185+t184+
        0.916*t3*t7*t127+2.137333333333334*t2*t27*t7*t127+0.305333333333333*
        t2*t3*t47*t25*t127)+0.00847*t129;
    d3fdrbrbab = 0.004235*t1*(t29*t30*t6*(8.0*t47*t8-2.0*
        t118*t5+t191)-0.458*t2*t3*t7*(t161-8.0*t8)+0.916*t3*t130*t7+
        2.137333333333334*t2*t27*t130*t7+4.274666666666668*t2*t106*
        t7+0.407111111111111*t4*t53*t117*t68-4.0*t3*t47*t49*t6+2.0*
        t131*t29*t49*t6+1.333333333333333*t27*t30*t47*t6-0.666666666666667*
        t131*t30*t40*t6-2.0*t118*t3*t30*t6+1.333333333333333*t27*t49*
        t54-2.0*t122*t3*t54+2.0*t131*t29*t30*t47+0.610666666666667*
        t2*t53*t47*t25+0.305333333333333*t2*t3*t130*t47*t25+1.424888888888889*
        t4*t106*t47*t25-0.305333333333333*t4*t53*t118*t25+t190+t189)+
        0.00847*t132;
    d3fdrbgaga = 0.004235*t1*(-0.458*t2*t3*(2.0*t11*t45+t196)*
        t7+t29*t30*(4.0*t18*t45*t5+t197)*t6+2.0*t3*t32*t49*t6+2.0*
        t3*t30*t32*t47+t195+t194+0.916*t27*t7*t18+0.305333333333333*
        t2*t27*t47*t25*t18)+t193;
    d3fdrbgagb = 0.0;
    d3fdrbgaab = 0.0;
    d3fdrbgbgb = 0.004235*t1*(-0.458*t2*t3*(t201-2.0*t11*
        t44)*t7+t29*t30*(4.0*t20*t44*t5+t202)*t6+2.0*t3*t34*t49*t6+
        2.0*t3*t30*t34*t47+t200-0.916*t27*t7*t20-0.305333333333333*
        t2*t27*t47*t25*t20+t199)+t198;
    d3fdrbgbab = 0.0;
    d3fdrbabab = 0.0;
    d3fdgagaga = 0.0;
    d3fdgagagb = 0.0;
    d3fdgagaab = 0.0;
    d3fdgagbgb = 0.0;
    d3fdgagbab = 0.0;
    d3fdgaabab = 0.0;
    d3fdgbgbgb = 0.0;
    d3fdgbgbab = 0.0;
    d3fdgbabab = 0.0;
    d3fdababab = 0.0;


    ds->df1000 += factor*dfdra;
    ds->df0100 += factor*dfdrb;
    ds->df0010 += factor*dfdga;
    ds->df0001 += factor*dfdgb;
    ds->df00001+= factor*dfdgab;

    ds->df2000 += factor*d2fdrara;
    ds->df1100 += factor*d2fdrarb;
    ds->df1010 += factor*d2fdraga;
    ds->df1001 += factor*d2fdragb;
    ds->df10001+= factor*d2fdraab;
    ds->df0200 += factor*d2fdrbrb;
    ds->df0110 += factor*d2fdrbga;
    ds->df0101 += factor*d2fdrbgb;
    ds->df01001+= factor*d2fdrbgab;
    ds->df0020 += factor*d2fdgaga;
    ds->df0011 += factor*d2fdgagb;
    ds->df00101+= factor*d2fdgagab;
    ds->df0002 += factor*d2fdgbgb;
    ds->df00011+= factor*d2fdgbgab;
    ds->df00002+= factor*d2fdgabgab;

    ds->df3000 += factor*d3fdrarara;
    ds->df2100  += factor*d3fdrararb;
    ds->df2010  += factor*d3fdraraga;
    ds->df2001  += factor*d3fdraragb;
    ds->df20001 += factor*d3fdraraab;
    ds->df1200  += factor*d3fdrarbrb;
    ds->df1110  += factor*d3fdrarbga;
    ds->df1101  += factor*d3fdrarbgb;
    ds->df11001 += factor*d3fdrarbab;
    ds->df1020  += factor*d3fdragaga;
    ds->df1011  += factor*d3fdragagb;
    ds->df10101 += factor*d3fdragaab;
    ds->df1002  += factor*d3fdragbgb;
    ds->df10011 += factor*d3fdragbab;
    ds->df10002 += factor*d3fdraabab;
    ds->df0300  += factor*d3fdrbrbrb;
    ds->df0210  += factor*d3fdrbrbga;
    ds->df0201  += factor*d3fdrbrbgb;
    ds->df02001 += factor*d3fdrbrbab;
    ds->df0120  += factor*d3fdrbgaga;
    ds->df0111  += factor*d3fdrbgagb;
    ds->df01101 += factor*d3fdrbgaab;
    ds->df0102  += factor*d3fdrbgbgb;
    ds->df01011 += factor*d3fdrbgbab;
    ds->df01002 += factor*d3fdrbabab;
    ds->df0030  += factor*d3fdgagaga;
    ds->df0021  += factor*d3fdgagagb;
    ds->df00201 += factor*d3fdgagaab;
    ds->df0012  += factor*d3fdgagbgb;
    ds->df00111 += factor*d3fdgagbab;
    ds->df00102 += factor*d3fdgaabab;
    ds->df0003  += factor*d3fdgbgbgb;
    ds->df00021 += factor*d3fdgbgbab;
    ds->df00012 += factor*d3fdgbabab;
    ds->df00003 += factor*d3fdababab;

}

static void
pw92ac_fourth(FunFourthFuncDrv *ds, real factor, const FunDensProp* dp)
{
    real dfdra, dfdrb, dfdga, dfdgb, dfdgab;
    real d2fdrara, d2fdrarb, d2fdraga, d2fdragb, d2fdraab, d2fdrbrb,
        d2fdrbga, d2fdrbgb, d2fdrbgab, d2fdgaga, d2fdgagb, d2fdgagab,
        d2fdgbgb, d2fdgbgab, d2fdgabgab;
    real d3fdrarara, d3fdrararb, d3fdraraga, d3fdraragb, d3fdraraab,
         d3fdrarbrb, d3fdrarbga, d3fdrarbgb, d3fdrarbab, d3fdragaga,
         d3fdragagb, d3fdragaab, d3fdragbgb, d3fdragbab, d3fdraabab,
         d3fdrbrbrb, d3fdrbrbga, d3fdrbrbgb, d3fdrbrbab, d3fdrbgaga,
         d3fdrbgagb, d3fdrbgaab, d3fdrbgbgb, d3fdrbgbab, d3fdrbabab,
         d3fdgagaga, d3fdgagagb, d3fdgagaab, d3fdgagbgb, d3fdgagbab,
         d3fdgaabab, d3fdgbgbgb, d3fdgbgbab, d3fdgbabab, d3fdababab;
    real d4fdrararara, d4fdrarararb, d4fdrararaga, d4fdrararagb,
         d4fdrararaab, d4fdrararbrb, d4fdrararbga, d4fdrararbgb, d4fdrararbab,
         d4fdraragaga, d4fdraragagb, d4fdraragaab, d4fdraragbgb, d4fdraragbab,
         d4fdraraabab, d4fdrarbrbrb, d4fdrarbrbga, d4fdrarbrbgb, d4fdrarbrbab,
         d4fdrarbgaga, d4fdrarbgagb, d4fdrarbgaab, d4fdrarbgbgb, d4fdrarbgbab,
         d4fdrarbabab, d4fdragagaga, d4fdragagagb, d4fdragagaab, d4fdragagbgb,
         d4fdragagbab, d4fdragaabab, d4fdragbgbgb, d4fdragbgbab, d4fdragbabab,
         d4fdraababab, d4fdrbrbrbrb, d4fdrbrbrbga, d4fdrbrbrbgb, d4fdrbrbrbab,
         d4fdrbrbgaga, d4fdrbrbgagb, d4fdrbrbgaab, d4fdrbrbgbgb, d4fdrbrbgbab,
         d4fdrbrbabab, d4fdrbgagaga, d4fdrbgagagb, d4fdrbgagaab, d4fdrbgagbgb,
         d4fdrbgagbab, d4fdrbgaabab, d4fdrbgbgbgb, d4fdrbgbgbab, d4fdrbgbabab,
         d4fdrbababab, d4fdgagagaga, d4fdgagagagb, d4fdgagagaab, d4fdgagagbgb,
         d4fdgagagbab, d4fdgagaabab, d4fdgagbgbgb, d4fdgagbgbab, d4fdgagbabab,
         d4fdgaababab, d4fdgbgbgbgb, d4fdgbgbgbab, d4fdgbgbabab, d4fdgbababab,
         d4fdabababab;
    real rhoa = dp->rhoa, rhob = dp->rhob;
    real grada = dp->grada, gradb = dp->gradb, gradab = dp->gradab;

    real t1, t2, t3, t4, t5, t6, t7, t8, t9, t10;
    real t11, t12, t13, t14, t15, t16, t17, t18;
    real t19, t20, t21, t22, t23, t24, t25, t26;
    real t27, t28, t29, t30, t31, t32, t33, t34;
    real t35, t36, t37, t38, t39, t40, t41, t42;
    real t43, t44, t45, t46, t47, t48, t49, t50;
    real t51, t52, t53, t54, t55, t56, t57, t58;
    real t59, t60, t61, t62, t63, t64, t65, t66;
    real t67, t68, t69, t70, t71, t72, t73, t74;
    real t75, t76, t77, t78, t79, t80, t81, t82;
    real t83, t84, t85, t86, t87, t88, t89, t90;
    real t91, t92, t93, t94, t95, t96, t97, t98;
    real t99, t100, t101, t102, t103, t104, t105;
    real t106, t107, t108, t109, t110, t111, t112;
    real t113, t114, t115, t116, t117, t118, t119;
    real t120, t121, t122, t123, t124, t125, t126;
    real t127, t128, t129, t130, t131, t132, t133;
    real t134, t135, t136, t137, t138, t139, t140;
    real t141, t142, t143, t144, t145, t146, t147;
    real t148, t149, t150, t151, t152, t153, t154;
    real t155, t156, t157, t158, t159, t160, t161;
    real t162, t163, t164, t165, t166, t167, t168;
    real t169, t170, t171, t172, t173, t174, t175;
    real t176, t177, t178, t179, t180, t181, t182;
    real t183, t184, t185, t186, t187, t188, t189;
    real t190, t191, t192, t193, t194, t195, t196;
    real t197, t198, t199, t200, t201, t202, t203;
    real t204, t205, t206, t207, t208, t209, t210;
    real t211, t212, t213, t214, t215, t216, t217;
    real t218, t219, t220, t221, t222, t223, t224;
    real t225, t226, t227, t228, t229, t230, t231;
    real t232, t233, t234, t235, t236, t237, t238;
    real t239, t240, t241, t242, t243, t244, t245;
    real t246, t247, t248, t249, t250, t251, t252;
    real t253, t254, t255, t256, t257, t258, t259;
    real t260, t261, t262, t263, t264, t265, t266;
    real t267, t268, t269, t270, t271, t272, t273;
    real t274, t275, t276, t277, t278, t279, t280;
    real t281, t282, t283, t284, t285, t286, t287;
    real t288, t289, t290, t291, t292, t293, t294;
    real t295, t296, t297, t298, t299, t300, t301;
    real t302, t303, t304, t305, t306, t307, t308;
    real t309, t310, t311, t312, t313, t314, t315;
    real t316, t317, t318, t319, t320, t321, t322;
    real t323, t324, t325, t326, t327, t328, t329;
    real t330, t331, t332, t333, t334, t335, t336;
    real t337, t338, t339, t340, t341, t342, t343;
    real t344, t345, t346, t347, t348, t349, t350;
    real t351, t352, t353, t354, t355, t356, t357;
    real t358, t359, t360, t361, t362, t363, t364;
    real t365, t366, t367, t368, t369, t370, t371;
    real t372, t373, t374, t375, t376, t377, t378;
    real t379, t380, t381, t382;

    t1 = rhob+rhoa;
    t2 = rhoa-1.0*rhob;
    t3 = 1/pow(t1,2.333333333333333);
    t4 = pow(t2,2.0);
    t5 = 1/pow(t1,2.0);
    t6 = 1.0-1.0*t4*t5;
    t7 = 1/pow(t6,0.333333333333333);
    t8 = 1/pow(t1,3.0);
    t9 = 4.0*t2*t8*gradab;
    t10 = pow(grada,2.0);
    t11 = 1/t1;
    t12 = t2*t5;
    t13 = -1.0*t11;
    t14 = t13+t12;
    t15 = pow(gradb,2.0);
    t16 = -1.0*t2*t5;
    t17 = t11+t16;
    t18 = 1.0-1.0*t11*t2;
    t19 = -1.0*t10*t18*t5;
    t20 = t2*t11+1.0;
    t21 = t15*t5*t20;
    t22 = -2.0*t5*gradab+t9+t21+t19-1.0*t11*t15*t17+t10*t11*
        t14;
    t23 = 2.0*t4*t8;
    t24 = t23-2.0*t2*t5;
    t25 = 1/pow(t6,1.333333333333333);
    t26 = -2.0*t2*t5*gradab-1.0*t11*t15*t20+t10*t11*t18;
    t27 = 1/
        pow(t1,3.333333333333333);
    t28 = 1.068666666666667*t2*t27*t7*t26;
    t29 = 1/pow(t1,0.333333333333333);
    t30 = 0.1*t4*t5-0.037;
    t31 = 4.0*t6*t8*gradab;
    t32 = pow(t18,2.0);
    t33 = -2.0*t10*t32*t8;
    t34 = pow(t20,2.0);
    t35 = -2.0*t15*t34*t8;
    t36 = -2.0*t24*t5*gradab+2.0*t15*t17*t20*t5+2.0*t10*t14*
        t18*t5+t35+t33+t31;
    t37 = -0.2*t4*t8;
    t38 = 0.2*t2*t5+t37;
    t39 = -2.0*t5*t6*gradab+t15*t5*t34+t10*t5*t32;
    t40 = 1/pow(t1,1.333333333333333);
    t41 = -0.333333333333333*t30*t39*t40*t6;
    t42 = t41+t29*t24*t30*t39+t29*t38*t6*t39+t29*t6*t30*t36-
        0.458*t3*t7*t26+t28+0.152666666666667*t2*t3*t24*t25*t26-0.458*
        t2*t3*t7*t22;
    t43 = 0.004235*(t29*t6*t30*t39-0.458*t2*t3*t7*t26);
    t44 = t13+
        t16;
    t45 = t11+t12;
    t46 = 2.0*t5*gradab+t9+t10*t11*t45-1.0*t11*t15*t44+t21+
        t19;
    t47 = 2.0*t2*t5+t23;
    t48 = -2.0*t47*t5*gradab+2.0*t10*t18*t45*t5+2.0*t15*t20*
        t44*t5+t35+t33+t31;
    t49 = t37-0.2*t2*t5;
    t50 = t41+t29*t47*t30*t39+t29*t49*t6*t39+t29*t6*t30*t48+
        0.458*t3*t7*t26+t28+0.152666666666667*t2*t3*t47*t25*t26-0.458*
        t2*t3*t7*t46;
    t51 = 2.0*t3*t30*t32*t6*grada-0.916*grada*t2*t27*t7*t18;
    t52 = 2.0*
        t3*t30*t34*t6*gradb+0.916*gradb*t2*t27*t7*t20;
    t53 = 1/pow(t1,4.333333333333333);
    t54 = pow(t6,2.0);
    t55 = 0.916*t4*t53*t7-2.0*t3*t30*t54;
    t56 = 1/pow(t1,4.0);
    t57 = -12.0*t2*t56*gradab;
    t58 = 2.0*t2*t8;
    t59 = -2.0*t5;
    t60 = t59+t58;
    t61 = -2.0*t2*t8;
    t62 = 2.0*t5;
    t63 = t62+t61;
    t64 = 2.0*t10*t18*t8;
    t65 = -2.0*t15*t20*t8;
    t66 = 8.0*t8*gradab+t65+t64+t10*t11*t63-1.0*t11*t15*t60+
        t57+2.0*t15*t17*t5-2.0*t10*t14*t5;
    t67 = pow(t24,2.0);
    t68 = 1/pow(t6,2.333333333333333);
    t69 = -6.0*t4*t56;
    t70 = 8.0*t2*t8+t69+t59;
    t71 = -3.562222222222223*t2*t53*t7*t26;
    t72 = -12.0*t56*t6*gradab;
    t73 = pow(t14,2.0);
    t74 = pow(t17,2.0);
    t75 = 6.0*t10*t32*t56;
    t76 = 6.0*t15*t34*t56;
    t77 = 8.0*t24*t8*gradab-2.0*t5*t70*gradab-8.0*t15*t17*
        t20*t8-8.0*t10*t14*t18*t8+t76+t75+2.0*t15*t5*t74+2.0*t10*t5*
        t73+t72+2.0*t10*t18*t5*t63+2.0*t15*t20*t5*t60;
    t78 = 0.6*t4*t56;
    t79 = 0.2*t5;
    t80 = t79-0.8*t2*t8+t78;
    t81 = 0.444444444444444*t3*t30*t39*t6;
    t82 = t81+t29*t6*t30*t77+t71-0.458*t2*t3*t7*t66-0.666666666666667*
        t38*t39*t40*t6-0.666666666666667*t30*t36*t40*t6+2.0*t29*t36*
        t38*t6-0.666666666666667*t24*t30*t39*t40+t29*t80*t6*t39+2.0*
        t24*t29*t38*t39+t29*t70*t30*t39+2.0*t24*t29*t30*t36+2.137333333333334*
        t27*t7*t26-0.203555555555556*t2*t3*t67*t68*t26+0.152666666666667*
        t2*t3*t70*t25*t26+0.305333333333333*t3*t24*t25*t26-0.712444444444445*
        t2*t27*t24*t25*t26-0.916*t3*t7*t22+2.137333333333334*t2*t27*
        t7*t22+0.305333333333333*t2*t3*t24*t25*t22;
    t83 = t65+t64+t57-2.0*t15*t2*t56-2.0*t10*t2*t56-1.0*t10*
        t45*t5-1.0*t10*t14*t5+t15*t5*t44+t15*t5*t17;
    t84 = t62+t69;
    t85 = 1/pow(t1,5.0);
    t86 = -2.0*t5*t84*gradab+4.0*t47*t8*gradab+4.0*t24*t8*
        gradab+4.0*t15*t2*t20*t85-4.0*t10*t18*t2*t85-4.0*t10*t18*t45*
        t8-4.0*t15*t20*t44*t8-4.0*t15*t17*t20*t8-4.0*t10*t14*t18*t8+
        t76+t75+t72+2.0*t10*t14*t45*t5+2.0*t15*t17*t44*t5;
    t87 = t78-0.2*t5;
    t88 = t29*t6*t30*t86-0.458*t2*t3*t7*t83+t81+t71-0.333333333333333*
        t39*t40*t49*t6-0.333333333333333*t30*t40*t48*t6-0.333333333333333*
        t38*t39*t40*t6-0.333333333333333*t30*t36*t40*t6+t29*t38*t6*
        t48+t29*t24*t30*t48-0.333333333333333*t30*t39*t40*t47-0.458*
        t3*t7*t46+1.068666666666667*t2*t27*t7*t46+0.152666666666667*
        t2*t3*t24*t25*t46-0.333333333333333*t24*t30*t39*t40+t29*t87*
        t6*t39+t29*t24*t49*t39+t29*t38*t47*t39+t29*t84*t30*t39+t29*
        t49*t6*t36+t29*t47*t30*t36-0.203555555555556*t2*t3*t24*t47*
        t68*t26+0.152666666666667*t2*t3*t84*t25*t26+0.152666666666667*
        t3*t47*t25*t26-0.356222222222222*t2*t27*t47*t25*t26-0.152666666666667*
        t3*t24*t25*t26-0.356222222222222*t2*t27*t24*t25*t26+0.458*
        t3*t7*t22+1.068666666666667*t2*t27*t7*t22+0.152666666666667*
        t2*t3*t47*t25*t22;
    t89 = 0.004235*t51;
    t90 = 2.137333333333334*grada*t2*t53*t7*t18;
    t91 = -0.666666666666667*t27*t30*t32*t6*grada;
    t92 = -2.0*t18*t5*grada;
    t93 = 2.0*t11*t14*grada+t92;
    t94 = -4.0*t32*t8*grada;
    t95 = 4.0*t14*t18*t5*grada+t94;
    t96 = 2.0*t3*t32*t38*t6*grada+2.0*t24*t3*t30*t32*grada+
        t29*t6*t30*t95-0.458*t2*t3*t7*t93+t91+t90-0.916*grada*t27*
        t7*t18+0.305333333333333*grada*t2*t27*t24*t25*t18;
    t97 = 0.004235*t52;
    t98 = -2.137333333333334*gradb*t2*t53*t7*t20;
    t99 = -0.666666666666667*t27*t30*t34*t6*gradb;
    t100 = 2.0*t20*t5*gradb;
    t101 = t100-2.0*t11*t17*gradb;
    t102 = -4.0*t34*t8*gradb;
    t103 = 4.0*t17*t20*t5*gradb+t102;
    t104 = 2.0*t3*t34*t38*t6*gradb+2.0*t24*t3*t30*t34*gradb+
        t99+t98+0.916*gradb*t27*t7*t20-0.305333333333333*gradb*t2*
        t27*t24*t25*t20+t29*t6*t30*t103-0.458*t2*t3*t7*t101;
    t105 = 0.004235*t55;
    t106 = 1/pow(t1,5.333333333333333);
    t107 = -2.137333333333334*t4*t106*t7;
    t108 = 4.0*t2*t8;
    t109 = t59+t108;
    t110 = 0.666666666666667*t27*t30*t54;
    t111 = 4.0*t6*t8;
    t112 = t111-2.0*t24*t5;
    t113 = 0.916*t2*t53*t7-0.458*t2*t3*t109*t7-2.0*t24*t3*
        t30*t6-2.0*t3*t38*t54-0.305333333333333*t4*t53*t24*t25+t29*
        t6*t30*t112+t110+t107;
    t114 = t59+t61;
    t115 = t62+t58;
    t116 = -8.0*t8*gradab+t65+t64+t57-2.0*t10*t45*t5+2.0*
        t15*t44*t5-1.0*t11*t115*t15+t10*t11*t114;
    t117 = pow(t47,2.0);
    t118 = -8.0*t2*t8+t69+t59;
    t119 = pow(t44,2.0);
    t120 = pow(t45,2.0);
    t121 = 8.0*t47*t8*gradab-2.0*t118*t5*gradab-8.0*t10*t18*
        t45*t8-8.0*t15*t20*t44*t8+t76+t75+t72+2.0*t115*t15*t20*t5+
        2.0*t10*t114*t18*t5+2.0*t119*t15*t5+2.0*t10*t120*t5;
    t122 = t79+0.8*t2*t8+t78;
    t123 = t81+t71+2.0*t29*t48*t49*t6-0.666666666666667*t39*
        t40*t49*t6-0.666666666666667*t30*t40*t48*t6+2.0*t29*t39*t47*
        t49+2.0*t29*t30*t47*t48-0.666666666666667*t30*t39*t40*t47+
        0.916*t3*t7*t46+2.137333333333334*t2*t27*t7*t46+0.305333333333333*
        t2*t3*t47*t25*t46+t29*t122*t6*t39+t29*t118*t30*t39-2.137333333333334*
        t27*t7*t26-0.203555555555556*t2*t3*t117*t68*t26-0.305333333333333*
        t3*t47*t25*t26-0.712444444444445*t2*t27*t47*t25*t26+0.152666666666667*
        t2*t3*t118*t25*t26+t29*t6*t30*t121-0.458*t2*t3*t7*t116;
    t124 = 2.0*t11*t45*grada+t92;
    t125 = 4.0*t18*t45*t5*grada+t94;
    t126 = 2.0*t3*t32*t49*t6*grada+2.0*t3*t30*t32*t47*grada+
        t91+t90+0.916*grada*t27*t7*t18+0.305333333333333*grada*t2*
        t27*t47*t25*t18+t29*t6*t30*t125-0.458*t2*t3*t7*t124;
    t127 = t100-2.0*t11*t44*gradb;
    t128 = 4.0*t20*t44*t5*gradb+t102;
    t129 = 2.0*t3*t34*t49*t6*gradb+2.0*t3*t30*t34*t47*gradb+
        t99+t98-0.916*gradb*t27*t7*t20-0.305333333333333*gradb*t2*
        t27*t47*t25*t20+t29*t6*t30*t128-0.458*t2*t3*t7*t127;
    t130 = t62+t108;
    t131 = t111-2.0*t47*t5;
    t132 = -0.916*t2*t53*t7-0.458*t2*t3*t130*t7-2.0*t3*t30*
        t47*t6-2.0*t3*t49*t54-0.305333333333333*t4*t53*t47*t25+t29*
        t6*t30*t131+t110+t107;
    t133 = 2.0*t3*t30*t32*t6-0.916*t2*t27*t7*t18;
    t134 = 2.0*t3*t30*t34*t6+0.916*t2*t27*t7*t20;
    t135 = 48.0*t2*t85*gradab;
    t136 = 6.0*t2*t56;
    t137 = -6.0*t8;
    t138 = t137+t136;
    t139 = -6.0*t2*t56;
    t140 = 6.0*t8;
    t141 = t140+t139;
    t142 = -6.0*t10*t18*t56;
    t143 = 6.0*t15*t20*t56;
    t144 = -36.0*t56*gradab-6.0*t15*t17*t8+6.0*t10*t14*t8-
        3.0*t10*t5*t63+3.0*t15*t5*t60-1.0*t11*t141*t15+t143+t142+t10*
        t11*t138+t135;
    t145 = pow(t24,3.0);
    t146 = 1/pow(t6,3.333333333333333);
    t147 = 24.0*t4*t85;
    t148 = 12.0*t8;
    t149 = -36.0*t2*t56+t148+t147;
    t150 = 15.4362962962963*t2*t106*t7*t26;
    t151 = 48.0*t6*t85*gradab;
    t152 = -24.0*t10*t32*t85;
    t153 = -24.0*t15*t34*t85;
    t154 = 12.0*t70*t8*gradab-36.0*t24*t56*gradab-2.0*t149*
        t5*gradab-12.0*t15*t74*t8-12.0*t10*t73*t8-12.0*t10*t18*t63*
        t8-12.0*t15*t20*t60*t8+6.0*t10*t14*t5*t63+6.0*t15*t17*t5*t60+
        36.0*t15*t17*t20*t56+36.0*t10*t14*t18*t56+2.0*t141*t15*t20*
        t5+2.0*t10*t138*t18*t5+t153+t152+t151;
    t155 = -2.4*t4*t85;
    t156 = -1.2*t8;
    t157 = t156+3.6*t2*t56+t155;
    t158 = -1.037037037037037*t27*t30*t39*t6;
    t159 = -1.0*t39*t40*t6*t80+3.0*t29*t36*t6*t80+3.0*t24*
        t29*t39*t80-1.0*t30*t40*t6*t77+3.0*t29*t38*t6*t77+3.0*t24*
        t29*t30*t77-1.0*t30*t39*t40*t70+3.0*t29*t38*t39*t70+3.0*t29*
        t30*t36*t70-1.374*t3*t7*t66+3.206*t2*t27*t7*t66+0.458*t2*t3*
        t24*t25*t66-2.0*t36*t38*t40*t6+1.333333333333333*t3*t38*t39*
        t6+1.333333333333333*t3*t30*t36*t6-2.0*t24*t38*t39*t40-2.0*
        t24*t30*t36*t40+t29*t157*t6*t39+1.333333333333333*t24*t3*t30*
        t39+t29*t149*t30*t39+6.0*t24*t29*t36*t38-10.68666666666667*
        t53*t7*t26-0.610666666666667*t3*t67*t68*t26+1.424888888888889*
        t2*t27*t67*t68*t26-0.610666666666667*t2*t3*t70*t24*t68*t26+
        0.458*t3*t70*t25*t26-1.068666666666667*t2*t27*t70*t25*t26+
        3.562222222222223*t2*t53*t24*t25*t26-2.137333333333334*t27*
        t24*t25*t26+0.152666666666667*t2*t3*t149*t25*t26+0.474962962962963*
        t2*t3*t145*t146*t26-10.68666666666667*t2*t53*t7*t22+6.412*
        t27*t7*t22-0.610666666666667*t2*t3*t67*t68*t22+0.458*t2*t3*
        t70*t25*t22+0.916*t3*t24*t25*t22-2.137333333333334*t2*t27*
        t24*t25*t22+t158+t29*t6*t30*t154+t150-0.458*t2*t3*t7*t144;
    t160 = t136-
        2.0*t8;
    t161 = 2.0*t8+t139;
    t162 = -12.0*t56*gradab+4.0*t15*t2*t85+4.0*t10*t2*t85+
        2.0*t10*t45*t8-2.0*t15*t44*t8-4.0*t15*t17*t8+4.0*t10*t14*t8-
        1.0*t10*t5*t63+t15*t5*t60-1.0*t11*t15*t161+t10*t11*t160+t143+
        t142+t135;
    t163 = 2.137333333333334*t2*t27*t7*t83;
    t164 = -0.712444444444445*t2*t27*t47*t25*t22;
    t165 = -0.712444444444445*t2*t27*t24*t25*t46;
    t166 = 0.949925925925926*t2*t27*t24*t47*t68*t26;
    t167 = -12.0*t2*t56;
    t168 = -4.0*t8;
    t169 = t168+t167+t147;
    t170 = -0.712444444444445*t2*t27*t84*t25*t26;
    t171 = 8.0*t8*t84*gradab;
    t172 = -8.0*t15*t17*t44*t8;
    t173 = -8.0*t10*t14*t45*t8;
    t174 = 1/pow(t1,6.0);
    t175 = 4.0*t70*t8*gradab-12.0*t47*t56*gradab-24.0*t24*
        t56*gradab-2.0*t169*t5*gradab+8.0*t15*t17*t2*t85-8.0*t10*t14*
        t2*t85-4.0*t15*t74*t8-4.0*t10*t73*t8-4.0*t10*t18*t63*t8-4.0*
        t15*t20*t60*t8+2.0*t10*t45*t5*t63+2.0*t15*t44*t5*t60+12.0*
        t10*t18*t45*t56+12.0*t15*t20*t44*t56+24.0*t15*t17*t20*t56+
        24.0*t10*t14*t18*t56+2.0*t15*t161*t20*t5+2.0*t10*t160*t18*
        t5-16.0*t15*t174*t2*t20+16.0*t10*t174*t18*t2+t173+t172+t171+
        t153+t152+t151;
    t176 = -0.666666666666667*t30*t40*t6*t86;
    t177 = -0.666666666666667*t38*t40*t48*t6;
    t178 = -0.666666666666667*t24*t30*t40*t48;
    t179 = -0.666666666666667*t36*t40*t49*t6;
    t180 = -0.666666666666667*t30*t36*t40*t47;
    t181 = -0.666666666666667*t24*t39*t40*t49;
    t182 = -0.666666666666667*t38*t39*t40*t47;
    t183 = 0.4*t8;
    t184 = t183+1.2*t2*t56+t155;
    t185 = -0.666666666666667*t39*t40*t6*t87;
    t186 = -0.666666666666667*t30*t39*t40*t84;
    t187 = 2.0*t29*t36*t6*t87+2.0*t24*t29*t39*t87+2.0*t29*
        t38*t6*t86+2.0*t24*t29*t30*t86+2.0*t29*t38*t39*t84+2.0*t29*
        t30*t36*t84-0.916*t3*t7*t83+0.305333333333333*t2*t3*t24*t25*
        t83-0.333333333333333*t39*t40*t6*t80+t29*t49*t6*t77-0.333333333333333*
        t30*t40*t6*t77+t29*t47*t30*t77-0.333333333333333*t30*t39*t40*
        t70+0.458*t3*t7*t66+1.068666666666667*t2*t27*t7*t66+0.152666666666667*
        t2*t3*t47*t25*t66+0.444444444444444*t3*t39*t49*t6+0.444444444444444*
        t3*t30*t48*t6-0.666666666666667*t36*t38*t40*t6+0.888888888888889*
        t3*t38*t39*t6+0.888888888888889*t3*t30*t36*t6+2.0*t24*t29*
        t36*t49+t29*t80*t6*t48+2.0*t24*t29*t38*t48+t29*t70*t30*t48+
        0.444444444444444*t3*t30*t39*t47+2.0*t29*t36*t38*t47-3.562222222222223*
        t2*t53*t7*t46+2.137333333333334*t27*t7*t46-0.203555555555556*
        t2*t3*t67*t68*t46+0.152666666666667*t2*t3*t70*t25*t46+0.305333333333333*
        t3*t24*t25*t46-0.666666666666667*t24*t38*t39*t40-0.666666666666667*
        t24*t30*t36*t40+t29*t184*t6*t39+t29*t70*t49*t39+t29*t80*t47*
        t39+0.888888888888889*t24*t3*t30*t39+t29*t169*t30*t39-3.562222222222223*
        t53*t7*t26+0.203555555555556*t3*t67*t68*t26+0.474962962962963*
        t2*t27*t67*t68*t26-0.203555555555556*t2*t3*t70*t47*t68*t26-
        0.407111111111111*t3*t24*t47*t68*t26-0.407111111111111*t2*
        t3*t84*t24*t68*t26+0.305333333333333*t3*t84*t25*t26-0.152666666666667*
        t3*t70*t25*t26-0.356222222222222*t2*t27*t70*t25*t26+1.187407407407408*
        t2*t53*t47*t25*t26-0.712444444444445*t27*t47*t25*t26+2.374814814814815*
        t2*t53*t24*t25*t26+0.152666666666667*t2*t3*t169*t25*t26+0.474962962962963*
        t2*t3*t67*t47*t146*t26-7.124444444444446*t2*t53*t7*t22-0.407111111111111*
        t2*t3*t24*t47*t68*t22+0.305333333333333*t2*t3*t84*t25*t22+
        0.305333333333333*t3*t47*t25*t22-0.305333333333333*t3*t24*
        t25*t22-0.712444444444445*t2*t27*t24*t25*t22+t186+t185+t182+
        t181+t180+t179+t178+t177+t176+t29*t6*t30*t175+t170+t166+t165+
        t164+t163-0.458*t2*t3*t7*t162+t158+t150;
    t188 = 0.00847*t88;
    t189 = -7.124444444444446*grada*t2*t106*t7*t18;
    t190 = 0.888888888888889*t30*t32*t53*t6*grada;
    t191 = 4.0*t18*t8*grada;
    t192 = 2.0*t11*t63*grada-4.0*t14*t5*grada+t191;
    t193 = 12.0*t32*t56*grada;
    t194 = -16.0*t14*t18*t8*grada+4.0*t5*t73*grada+4.0*t18*
        t5*t63*grada+t193;
    t195 = 2.0*t3*t32*t6*t80*grada+2.0*t3*t30*t32*t70*grada-
        1.333333333333333*t27*t32*t38*t6*grada+4.0*t24*t3*t32*t38*
        grada-1.333333333333333*t24*t27*t30*t32*grada-0.666666666666667*
        t30*t40*t6*t95+2.0*t29*t38*t6*t95+2.0*t24*t29*t30*t95-0.916*
        t3*t7*t93+2.137333333333334*t2*t27*t7*t93+0.305333333333333*
        t2*t3*t24*t25*t93+t29*t6*t30*t194-0.458*t2*t3*t7*t192+t190+
        t189+4.274666666666668*grada*t53*t7*t18-0.407111111111111*
        grada*t2*t27*t67*t68*t18+0.305333333333333*grada*t2*t27*t70*
        t25*t18-1.424888888888889*grada*t2*t53*t24*t25*t18+0.610666666666667*
        grada*t27*t24*t25*t18;
    t196 = 7.124444444444446*gradb*t2*t106*t7*t20;
    t197 = 0.888888888888889*t30*t34*t53*t6*gradb;
    t198 = -4.0*t20*t8*gradb;
    t199 = -2.0*t11*t60*gradb+4.0*t17*t5*gradb+t198;
    t200 = 12.0*t34*t56*gradb;
    t201 = -16.0*t17*t20*t8*gradb+4.0*t5*t74*gradb+4.0*t20*
        t5*t60*gradb+t200;
    t202 = 2.0*t3*t34*t6*t80*gradb+2.0*t3*t30*t34*t70*gradb-
        1.333333333333333*t27*t34*t38*t6*gradb+4.0*t24*t3*t34*t38*
        gradb-1.333333333333333*t24*t27*t30*t34*gradb-0.666666666666667*
        t103*t30*t40*t6+2.0*t103*t29*t38*t6+2.0*t103*t24*t29*t30+t29*
        t6*t30*t201-4.274666666666668*gradb*t53*t7*t20+0.407111111111111*
        gradb*t2*t27*t67*t68*t20-0.305333333333333*gradb*t2*t27*t70*
        t25*t20+1.424888888888889*gradb*t2*t53*t24*t25*t20-0.610666666666667*
        gradb*t27*t24*t25*t20-0.458*t2*t3*t7*t199+t197+t196-0.916*
        t3*t7*t101+2.137333333333334*t2*t27*t7*t101+0.305333333333333*
        t2*t3*t24*t25*t101;
    t203 = 1/pow(t1,6.333333333333333);
    t204 = 7.124444444444446*t4*t203*t7;
    t205 = 8.0*t8+t167;
    t206 = -0.888888888888889*t30*t53*t54;
    t207 = -12.0*t56*t6;
    t208 = 8.0*t24*t8-2.0*t5*t70+t207;
    t209 = -2.0*t3*t54*t80-2.0*t3*t30*t6*t70-0.458*t2*t3*
        t205*t7-0.916*t3*t109*t7+2.137333333333334*t2*t27*t109*t7-
        4.274666666666668*t2*t106*t7+0.407111111111111*t4*t53*t67*
        t68-0.666666666666667*t112*t30*t40*t6-4.0*t24*t3*t38*t6+2.0*
        t112*t29*t38*t6+1.333333333333333*t24*t27*t30*t6+1.333333333333333*
        t27*t38*t54+2.0*t112*t24*t29*t30-0.305333333333333*t4*t53*
        t70*t25-0.610666666666667*t2*t53*t24*t25+0.305333333333333*
        t2*t3*t109*t24*t25+1.424888888888889*t4*t106*t24*t25+t29*t6*
        t30*t208+t206+t204;
    t210 = 12.0*t56*gradab+10.0*t15*t2*t85+10.0*t10*t2*t85+
        4.0*t10*t45*t8-4.0*t15*t44*t8-2.0*t15*t17*t8+2.0*t10*t14*t8+
        2.0*t15*t56+2.0*t10*t56-1.0*t10*t114*t5+t143+t142+t135+t15*
        t5*t115;
    t211 = 12.0*t2*t56+t168+t147;
    t212 = 4.0*t118*t8*gradab-24.0*t47*t56*gradab-12.0*t24*
        t56*gradab-2.0*t211*t5*gradab-8.0*t10*t2*t45*t85+8.0*t15*t2*
        t44*t85-4.0*t15*t20*t85+4.0*t10*t18*t85-4.0*t115*t15*t20*t8-
        4.0*t10*t114*t18*t8-4.0*t119*t15*t8-4.0*t10*t120*t8+24.0*t10*
        t18*t45*t56+24.0*t15*t20*t44*t56+12.0*t15*t17*t20*t56+12.0*
        t10*t14*t18*t56+2.0*t115*t15*t17*t5+2.0*t10*t114*t14*t5-28.0*
        t15*t174*t2*t20+28.0*t10*t174*t18*t2+t173+t172+t171+t153+t152+
        t151;
    t213 = t183-1.2*t2*t56+t155;
    t214 = 2.0*t29*t48*t6*t87+2.0*t29*t39*t47*t87+2.0*t29*
        t49*t6*t86+2.0*t29*t30*t47*t86+2.0*t29*t39*t49*t84+2.0*t29*
        t30*t48*t84+0.916*t3*t7*t83+0.305333333333333*t2*t3*t47*t25*
        t83-0.666666666666667*t40*t48*t49*t6+0.888888888888889*t3*
        t39*t49*t6+0.888888888888889*t3*t30*t48*t6-0.333333333333333*
        t122*t39*t40*t6-0.333333333333333*t121*t30*t40*t6+0.444444444444444*
        t3*t38*t39*t6+0.444444444444444*t3*t30*t36*t6+2.0*t24*t29*
        t48*t49-0.666666666666667*t39*t40*t47*t49+2.0*t29*t36*t47*
        t49-0.666666666666667*t30*t40*t47*t48+2.0*t29*t38*t47*t48+
        0.888888888888889*t3*t30*t39*t47-7.124444444444446*t2*t53*
        t7*t46-0.407111111111111*t2*t3*t24*t47*t68*t46+0.305333333333333*
        t2*t3*t84*t25*t46+0.305333333333333*t3*t47*t25*t46-0.712444444444445*
        t2*t27*t47*t25*t46-0.305333333333333*t3*t24*t25*t46-0.333333333333333*
        t118*t30*t39*t40+t29*t213*t6*t39+t29*t118*t38*t39+0.444444444444444*
        t24*t3*t30*t39+t29*t211*t30*t39+t29*t122*t24*t39+t29*t122*
        t6*t36+t29*t118*t30*t36+3.562222222222223*t53*t7*t26-0.407111111111111*
        t2*t3*t84*t47*t68*t26+0.407111111111111*t3*t24*t47*t68*t26-
        0.203555555555556*t2*t3*t118*t24*t68*t26-0.203555555555556*
        t3*t117*t68*t26+0.474962962962963*t2*t27*t117*t68*t26-0.305333333333333*
        t3*t84*t25*t26+2.374814814814815*t2*t53*t47*t25*t26+1.187407407407408*
        t2*t53*t24*t25*t26+0.712444444444445*t27*t24*t25*t26+0.152666666666667*
        t2*t3*t211*t25*t26+0.152666666666667*t3*t118*t25*t26-0.356222222222222*
        t2*t27*t118*t25*t26+0.474962962962963*t2*t3*t24*t117*t146*
        t26-3.562222222222223*t2*t53*t7*t22-2.137333333333334*t27*
        t7*t22-0.203555555555556*t2*t3*t117*t68*t22-0.305333333333333*
        t3*t47*t25*t22+0.152666666666667*t2*t3*t118*t25*t22+t29*t6*
        t30*t212-0.458*t2*t3*t7*t210+t186+t185+t182+t181+t180+t179+
        t178+t177+t176+t170+t166+t165+t164+t163+t158+t150+t29*t38*
        t6*t121+t29*t24*t30*t121-0.458*t3*t7*t116+1.068666666666667*
        t2*t27*t7*t116+0.152666666666667*t2*t3*t24*t25*t116;
    t215 = -4.0*t2*t56*grada-2.0*t45*t5*grada-2.0*t14*t5*
        grada+t191;
    t216 = -8.0*t18*t2*t85*grada-8.0*t18*t45*t8*grada-8.0*
        t14*t18*t8*grada+4.0*t14*t45*t5*grada+t193;
    t217 = 2.0*t3*t32*t6*t87*grada+2.0*t3*t30*t32*t84*grada-
        0.666666666666667*t27*t32*t49*t6*grada-0.666666666666667*t27*
        t32*t38*t6*grada+2.0*t24*t3*t32*t49*grada+2.0*t3*t32*t38*t47*
        grada-0.666666666666667*t27*t30*t32*t47*grada-0.666666666666667*
        t24*t27*t30*t32*grada+t29*t49*t6*t95-0.333333333333333*t30*
        t40*t6*t95+t29*t47*t30*t95+0.458*t3*t7*t93+1.068666666666667*
        t2*t27*t7*t93+0.152666666666667*t2*t3*t47*t25*t93-0.333333333333333*
        t125*t30*t40*t6+t29*t6*t30*t216-0.458*t2*t3*t7*t215+t190+t189-
        0.407111111111111*grada*t2*t27*t24*t47*t68*t18+0.305333333333333*
        grada*t2*t27*t84*t25*t18-0.712444444444445*grada*t2*t53*t47*
        t25*t18+0.305333333333333*grada*t27*t47*t25*t18-0.712444444444445*
        grada*t2*t53*t24*t25*t18-0.305333333333333*grada*t27*t24*t25*
        t18+t29*t38*t6*t125+t29*t24*t30*t125-0.458*t3*t7*t124+1.068666666666667*
        t2*t27*t7*t124+0.152666666666667*t2*t3*t24*t25*t124;
    t218 = -4.0*t2*t56*gradb+2.0*t44*t5*gradb+2.0*t17*t5*
        gradb+t198;
    t219 = 8.0*t2*t20*t85*gradb-8.0*t20*t44*t8*gradb-8.0*
        t17*t20*t8*gradb+4.0*t17*t44*t5*gradb+t200;
    t220 = 2.0*t3*t34*t6*t87*gradb+2.0*t3*t30*t34*t84*gradb-
        0.666666666666667*t27*t34*t49*t6*gradb-0.666666666666667*t27*
        t34*t38*t6*gradb+2.0*t24*t3*t34*t49*gradb+2.0*t3*t34*t38*t47*
        gradb-0.666666666666667*t27*t30*t34*t47*gradb-0.666666666666667*
        t24*t27*t30*t34*gradb-0.333333333333333*t128*t30*t40*t6-0.333333333333333*
        t103*t30*t40*t6+t29*t6*t30*t219-0.458*t2*t3*t7*t218+0.407111111111111*
        gradb*t2*t27*t24*t47*t68*t20-0.305333333333333*gradb*t2*t27*
        t84*t25*t20+0.712444444444445*gradb*t2*t53*t47*t25*t20-0.305333333333333*
        gradb*t27*t47*t25*t20+0.712444444444445*gradb*t2*t53*t24*t25*
        t20+0.305333333333333*gradb*t27*t24*t25*t20+t197+t196+t29*
        t38*t6*t128+t29*t24*t30*t128-0.458*t3*t7*t127+1.068666666666667*
        t2*t27*t7*t127+0.152666666666667*t2*t3*t24*t25*t127+t29*t49*
        t6*t103+t29*t47*t30*t103+0.458*t3*t7*t101+1.068666666666667*
        t2*t27*t7*t101+0.152666666666667*t2*t3*t47*t25*t101;
    t221 = -2.0*t5*t84+4.0*t47*t8+4.0*t24*t8+t207;
    t222 = -2.0*t3*t54*t87-2.0*t3*t30*t6*t84+12.62044444444445*
        t4*t203*t7-0.458*t3*t130*t7+1.068666666666667*t2*t27*t130*
        t7+0.458*t3*t109*t7+1.068666666666667*t2*t27*t109*t7+0.407111111111111*
        t4*t53*t24*t47*t68-2.0*t24*t3*t49*t6-2.0*t3*t38*t47*t6+0.666666666666667*
        t27*t30*t47*t6-0.333333333333333*t131*t30*t40*t6-0.333333333333333*
        t112*t30*t40*t6+0.666666666666667*t24*t27*t30*t6+0.666666666666667*
        t27*t49*t54+0.666666666666667*t27*t38*t54-0.305333333333333*
        t4*t53*t84*t25-0.305333333333333*t2*t53*t47*t25+0.152666666666667*
        t2*t3*t109*t47*t25+0.712444444444445*t4*t106*t47*t25+0.305333333333333*
        t2*t53*t24*t25+0.152666666666667*t2*t3*t130*t24*t25+0.712444444444445*
        t4*t106*t24*t25+t29*t6*t30*t221+t206+t29*t38*t6*t131+t29*t24*
        t30*t131+t29*t49*t6*t112+t29*t47*t30*t112;
    t223 = 0.004235*t133;
    t224 = 2.137333333333334*t2*t53*t7*t18;
    t225 = -0.666666666666667*t27*t30*t32*t6;
    t226 = -2.0*t18*t5;
    t227 = t226+2.0*t11*t14;
    t228 = -4.0*t32*t8;
    t229 = 4.0*t14*t18*t5+t228;
    t230 = 2.0*t3*t32*t38*t6+2.0*t24*t3*t30*t32+t29*t6*t30*
        t229-0.458*t2*t3*t7*t227+t225+t224-0.916*t27*t7*t18+0.305333333333333*
        t2*t27*t24*t25*t18;
    t231 = 0.004235*t134;
    t232 = -2.137333333333334*t2*t53*t7*t20;
    t233 = -0.666666666666667*t27*t30*t34*t6;
    t234 = 2.0*t20*t5;
    t235 = t234-2.0*t11*t17;
    t236 = -4.0*t34*t8;
    t237 = 4.0*t17*t20*t5+t236;
    t238 = 2.0*t3*t34*t38*t6+2.0*t24*t3*t30*t34+t29*t6*t30*
        t237-0.458*t2*t3*t7*t235+t233+t232+0.916*t27*t7*t20-0.305333333333333*
        t2*t27*t24*t25*t20;
    t239 = t137+t139;
    t240 = t140+t136;
    t241 = 36.0*t56*gradab+6.0*t10*t45*t8-6.0*t15*t44*t8+
        3.0*t115*t15*t5-3.0*t10*t114*t5+t10*t11*t240-1.0*t11*t15*t239+
        t143+t142+t135;
    t242 = pow(t47,3.0);
    t243 = 36.0*t2*t56+t148+t147;
    t244 = 12.0*t118*t8*gradab-36.0*t47*t56*gradab-2.0*t243*
        t5*gradab-12.0*t115*t15*t20*t8-12.0*t10*t114*t18*t8-12.0*t119*
        t15*t8-12.0*t10*t120*t8+36.0*t10*t18*t45*t56+36.0*t15*t20*
        t44*t56+6.0*t10*t114*t45*t5+6.0*t115*t15*t44*t5+2.0*t10*t18*
        t240*t5+2.0*t15*t20*t239*t5+t153+t152+t151;
    t245 = t156-3.6*t2*t56+t155;
    t246 = -2.0*t40*t48*t49*t6+1.333333333333333*t3*t39*t49*
        t6+3.0*t121*t29*t49*t6+1.333333333333333*t3*t30*t48*t6+3.0*
        t122*t29*t48*t6-1.0*t122*t39*t40*t6-1.0*t121*t30*t40*t6+6.0*
        t29*t47*t48*t49-2.0*t39*t40*t47*t49+3.0*t118*t29*t39*t49-2.0*
        t30*t40*t47*t48+3.0*t118*t29*t30*t48+1.333333333333333*t3*
        t30*t39*t47+3.0*t122*t29*t39*t47+3.0*t121*t29*t30*t47-10.68666666666667*
        t2*t53*t7*t46-6.412*t27*t7*t46-0.610666666666667*t2*t3*t117*
        t68*t46-0.916*t3*t47*t25*t46-2.137333333333334*t2*t27*t47*
        t25*t46+0.458*t2*t3*t118*t25*t46-1.0*t118*t30*t39*t40+t29*
        t245*t6*t39+t29*t243*t30*t39+10.68666666666667*t53*t7*t26-
        0.610666666666667*t2*t3*t118*t47*t68*t26+0.610666666666667*
        t3*t117*t68*t26+1.424888888888889*t2*t27*t117*t68*t26+3.562222222222223*
        t2*t53*t47*t25*t26+2.137333333333334*t27*t47*t25*t26+0.152666666666667*
        t2*t3*t243*t25*t26-0.458*t3*t118*t25*t26-1.068666666666667*
        t2*t27*t118*t25*t26+0.474962962962963*t2*t3*t242*t146*t26+
        t29*t6*t30*t244-0.458*t2*t3*t7*t241+t158+t150+1.374*t3*t7*
        t116+3.206*t2*t27*t7*t116+0.458*t2*t3*t47*t25*t116;
    t247 = -4.0*t45*t5*grada+2.0*t11*t114*grada+t191;
    t248 = -16.0*t18*t45*t8*grada+4.0*t114*t18*t5*grada+4.0*
        t120*t5*grada+t193;
    t249 = -1.333333333333333*t27*t32*t49*t6*grada+2.0*t122*
        t3*t32*t6*grada+4.0*t3*t32*t47*t49*grada-1.333333333333333*
        t27*t30*t32*t47*grada+2.0*t118*t3*t30*t32*grada+2.0*t125*t29*
        t49*t6-0.666666666666667*t125*t30*t40*t6+2.0*t125*t29*t30*
        t47+t29*t6*t30*t248-0.458*t2*t3*t7*t247+t190+t189-4.274666666666668*
        grada*t53*t7*t18-0.407111111111111*grada*t2*t27*t117*t68*t18-
        1.424888888888889*grada*t2*t53*t47*t25*t18-0.610666666666667*
        grada*t27*t47*t25*t18+0.305333333333333*grada*t2*t27*t118*
        t25*t18+0.916*t3*t7*t124+2.137333333333334*t2*t27*t7*t124+
        0.305333333333333*t2*t3*t47*t25*t124;
    t250 = 4.0*t44*t5*gradb-2.0*t11*t115*gradb+t198;
    t251 = -16.0*t20*t44*t8*gradb+4.0*t115*t20*t5*gradb+4.0*
        t119*t5*gradb+t200;
    t252 = -1.333333333333333*t27*t34*t49*t6*gradb+2.0*t122*
        t3*t34*t6*gradb+4.0*t3*t34*t47*t49*gradb-1.333333333333333*
        t27*t30*t34*t47*gradb+2.0*t118*t3*t30*t34*gradb+2.0*t128*t29*
        t49*t6-0.666666666666667*t128*t30*t40*t6+2.0*t128*t29*t30*
        t47+t29*t6*t30*t251-0.458*t2*t3*t7*t250+4.274666666666668*
        gradb*t53*t7*t20+0.407111111111111*gradb*t2*t27*t117*t68*t20+
        1.424888888888889*gradb*t2*t53*t47*t25*t20+0.610666666666667*
        gradb*t27*t47*t25*t20-0.305333333333333*gradb*t2*t27*t118*
        t25*t20+t197+t196+0.916*t3*t7*t127+2.137333333333334*t2*t27*
        t7*t127+0.305333333333333*t2*t3*t47*t25*t127;
    t253 = t167-8.0*t8;
    t254 = 8.0*t47*t8-2.0*t118*t5+t207;
    t255 = -0.458*t2*t3*t253*t7+0.916*t3*t130*t7+2.137333333333334*
        t2*t27*t130*t7+4.274666666666668*t2*t106*t7+0.407111111111111*
        t4*t53*t117*t68-4.0*t3*t47*t49*t6+2.0*t131*t29*t49*t6+1.333333333333333*
        t27*t30*t47*t6-0.666666666666667*t131*t30*t40*t6-2.0*t118*
        t3*t30*t6+1.333333333333333*t27*t49*t54-2.0*t122*t3*t54+2.0*
        t131*t29*t30*t47+t29*t6*t30*t254+0.610666666666667*t2*t53*
        t47*t25+0.305333333333333*t2*t3*t130*t47*t25+1.424888888888889*
        t4*t106*t47*t25-0.305333333333333*t4*t53*t118*t25+t206+t204;
    t256 = 2.0*
        t11*t45+t226;
    t257 = 4.0*t18*t45*t5+t228;
    t258 = 2.0*t3*t32*t49*t6+2.0*t3*t30*t32*t47+t29*t6*t30*
        t257-0.458*t2*t3*t7*t256+t225+t224+0.916*t27*t7*t18+0.305333333333333*
        t2*t27*t47*t25*t18;
    t259 = t234-2.0*t11*t44;
    t260 = 4.0*t20*t44*t5+t236;
    t261 = 2.0*t3*t34*t49*t6+2.0*t3*t30*t34*t47+t29*t6*t30*
        t260-0.458*t2*t3*t7*t259+t233+t232-0.916*t27*t7*t20-0.305333333333333*
        t2*t27*t47*t25*t20;
    t262 = -240.0*t174*t2*gradab;
    t263 = 24.0*t2*t85;
    t264 = -24.0*t56;
    t265 = t264+t263;
    t266 = -24.0*t2*t85;
    t267 = 24.0*t56;
    t268 = t267+t266;
    t269 = 24.0*t10*t18*t85;
    t270 = -24.0*t15*t20*t85;
    t271 = 1/pow(t6,4.333333333333333);
    t272 = -120.0*t174*t4;
    t273 = -72.0*t56;
    t274 = 192.0*t2*t85+t273+t272;
    t275 = -82.32691358024692*t2*t203*t7*t26;
    t276 = -240.0*t174*t6*gradab;
    t277 = 120.0*t10*t174*t32;
    t278 = 120.0*t15*t174*t34;
    t279 = 12.0*t4*t174;
    t280 = 7.200000000000001*t56;
    t281 = 3.45679012345679*t30*t39*t53*t6;
    t282 = -12.0*t56;
    t283 = t282+t263;
    t284 = 12.0*t56;
    t285 = t284+t266;
    t286 = -10.68666666666667*t2*t53*t7*t83;
    t287 = 3.562222222222223*t2*t53*t47*t25*t22;
    t288 = 3.562222222222223*t2*t53*t24*t25*t46;
    t289 = -4.74962962962963*t2*t53*t24*t47*t68*t26;
    t290 = 96.0*t2*t85+t272;
    t291 = 3.562222222222223*t2*t53*t84*t25*t26;
    t292 = -36.0*t56*t84*gradab;
    t293 = 36.0*t15*t17*t44*t56;
    t294 = 36.0*t10*t14*t45*t56;
    t295 = 1/pow(t1,7.0);
    t296 = 1.333333333333333*t3*t30*t6*t86;
    t297 = 1.333333333333333*t3*t38*t48*t6;
    t298 = 1.333333333333333*t24*t3*t30*t48;
    t299 = 1.333333333333333*t3*t36*t49*t6;
    t300 = 1.333333333333333*t3*t30*t36*t47;
    t301 = 1.333333333333333*t24*t3*t39*t49;
    t302 = 1.333333333333333*t3*t38*t39*t47;
    t303 = 1.333333333333333*t3*t39*t6*t87;
    t304 = 1.333333333333333*t3*t30*t39*t84;
    t305 = 30.8725925925926*grada*t2*t203*t7*t18;
    t306 = -2.074074074074074*t106*t30*t32*t6*grada;
    t307 = -12.0*t18*t56*grada;
    t308 = -48.0*t32*t85*grada;
    t309 = -30.8725925925926*gradb*t2*t203*t7*t20;
    t310 = -2.074074074074074*t106*t30*t34*t6*gradb;
    t311 = 12.0*t20*t56*gradb;
    t312 = -48.0*t34*t85*gradb;
    t313 = 1/pow(t1,7.333333333333333);
    t314 = -30.8725925925926*t4*t313*t7;
    t315 = 48.0*t2*t85;
    t316 = 2.074074074074074*t106*t30*t54;
    t317 = 48.0*t6*t85;
    t318 = t267+t272;
    t319 = 1/pow(t1,8.0);
    t320 = -192.0*t10*t18*t2*t295;
    t321 = 192.0*t15*t2*t20*t295;
    t322 = 1.899851851851852*grada*t2*t53*t24*t47*t68*t18;
    t323 = -
        1.424888888888889*grada*t2*t53*t84*t25*t18;
    t324 = -1.333333333333333*t24*t27*t32*t49*grada;
    t325 = -1.333333333333333*t27*t32*t38*t47*grada;
    t326 = -1.333333333333333*t27*t32*t6*t87*grada;
    t327 = -1.333333333333333*t27*t30*t32*t84*grada;
    t328 = 2.137333333333334*t2*t27*t7*t215;
    t329 = -0.712444444444445*t2*t27*t47*t25*t93;
    t330 = -0.712444444444445*t2*t27*t24*t25*t124;
    t331 = -16.0*t14*t45*t8*grada;
    t332 = -0.666666666666667*t216*t30*t40*t6;
    t333 = -0.666666666666667*t40*t49*t6*t95;
    t334 = -0.666666666666667*t30*t40*t47*t95;
    t335 = -0.666666666666667*t125*t38*t40*t6;
    t336 = -0.666666666666667*t125*t24*t30*t40;
    t337 = 0.00847*t217;
    t338 = -1.899851851851852*gradb*t2*t53*t24*t47*t68*t20;
    t339 = 1.424888888888889*
        gradb*t2*t53*t84*t25*t20;
    t340 = -1.333333333333333*t24*t27*t34*t49*gradb;
    t341 = -1.333333333333333*t27*t34*t38*t47*gradb;
    t342 = -1.333333333333333*t27*t34*t6*t87*gradb;
    t343 = -1.333333333333333*t27*t30*t34*t84*gradb;
    t344 = 2.137333333333334*t2*t27*t7*t218;
    t345 = -0.712444444444445*t2*t27*t24*t25*t127;
    t346 = -0.712444444444445*t2*t27*t47*t25*t101;
    t347 = -16.0*t17*t44*t8*gradb;
    t348 = -0.666666666666667*t219*t30*t40*t6;
    t349 = -0.666666666666667*t128*t38*t40*t6;
    t350 = -0.666666666666667*t128*t24*t30*t40;
    t351 = -0.666666666666667*t103*t40*t49*t6;
    t352 = -0.666666666666667*t103*t30*t40*t47;
    t353 = 0.00847*t220;
    t354 = -1.899851851851852*t4*t106*t24*t47*t68;
    t355 = 1.424888888888889*t4*t106*t84*t25;
    t356 = -0.712444444444445*t2*t27*t130*t24*t25;
    t357 = -0.712444444444445*t2*t27*t109*t47*t25;
    t358 = -56.5205925925926*t4*t313*t7;
    t359 = 1.333333333333333*t24*t27*t49*t6;
    t360 = 1.333333333333333*t27*t38*t47*t6;
    t361 = 1.333333333333333*t27*t54*t87;
    t362 = 1.333333333333333*t27*t30*t6*t84;
    t363 = 8.0*t8*t84;
    t364 = -0.666666666666667*t221*t30*t40*t6;
    t365 = -0.666666666666667*t112*t40*t49*t6;
    t366 = -0.666666666666667*t112*t30*t40*t47;
    t367 = -0.666666666666667*t131*t38*t40*t6;
    t368 = -0.666666666666667*t131*t24*t30*t40;
    t369 = 0.00847*t222;
    t370 = -7.124444444444446*t2*t106*t7*t18;
    t371 = 0.888888888888889*t30*t32*t53*t6;
    t372 = 4.0*t18*t8;
    t373 = 12.0*t32*t56;
    t374 = 7.124444444444446*t2*t106*t7*t20;
    t375 = 0.888888888888889*t30*t34*t53*t6;
    t376 = -4.0*t20*t8;
    t377 = 12.0*t34*t56;
    t378 = t272-96.0*t2*t85;
    t379 = -4.0*t2*t56;
    t380 = t264+t266;
    t381 = t267+t263;
    t382 = -192.0*t2*t85+t273+t272;

   /* code */
    dfdra = t43+0.004235*t1*t42;
    dfdrb = t43+0.004235*t1*t50;
    dfdga = 0.004235*t1*t51;
    dfdgb = 0.004235*t1*t52;
    dfdgab = 0.004235*t1*t55;
    d2fdrara = 0.00847*t42+0.004235*t1*t82;
    d2fdrarb = 0.004235*t50+0.004235*t42+0.004235*t1*t88;
    d2fdraga = 0.004235*
        t1*t96+t89;
    d2fdragb = 0.004235*t1*t104+t97;
    d2fdraab = 0.004235*t1*t113+t105;
    d2fdrbrb = 0.00847*t50+0.004235*t1*t123;
    d2fdrbga = 0.004235*t1*t126+t89;
    d2fdrbgb = 0.004235*t1*t129+t97;
    d2fdrbgab = 0.004235*t1*t132+t105;
    d2fdgaga = 0.004235*t1*t133;
    d2fdgagb = 0.0;
    d2fdgagab = 0.0;
    d2fdgbgb = 0.004235*t1*t134;
    d2fdgbgab = 0.0;
    d2fdgabgab = 0.0;
    d3fdrarara = 0.012705*t82+0.004235*t1*t159;
    d3fdrararb = t188+0.004235*t82+0.004235*t1*t187;
    d3fdraraga = 0.00847*t96+0.004235*t1*t195;
    d3fdraragb = 0.00847*t104+0.004235*t1*t202;
    d3fdraraab = 0.00847*t113+0.004235*t1*t209;
    d3fdrarbrb = t188+0.004235*t123+0.004235*t1*t214;
    d3fdrarbga = 0.004235*t126+0.004235*t1*t217+0.004235*
        t96;
    d3fdrarbgb = 0.004235*t104+0.004235*t1*t220+0.004235*
        t129;
    d3fdrarbab = 0.004235*t132+0.004235*t1*t222+0.004235*
        t113;
    d3fdragaga = 0.004235*t1*t230+t223;
    d3fdragagb = 0.0;
    d3fdragaab = 0.0;
    d3fdragbgb = 0.004235*t1*t238+t231;
    d3fdragbab = 0.0;
    d3fdraabab = 0.0;
    d3fdrbrbrb = 0.012705*t123+0.004235*t1*t246;
    d3fdrbrbga = 0.00847*t126+0.004235*t1*t249;
    d3fdrbrbgb = 0.00847*t129+0.004235*t1*t252;
    d3fdrbrbab = 0.00847*t132+0.004235*t1*t255;
    d3fdrbgaga = 0.004235*t1*t258+t223;
    d3fdrbgagb = 0.0;
    d3fdrbgaab = 0.0;
    d3fdrbgbgb = 0.004235*t1*t261+t231;
    d3fdrbgbab = 0.0;
    d3fdrbabab = 0.0;
    d3fdgagaga = 0.0;
    d3fdgagagb = 0.0;
    d3fdgagaab = 0.0;
    d3fdgagbgb = 0.0;
    d3fdgagbab = 0.0;
    d3fdgaabab = 0.0;
    d3fdgbgbgb = 0.0;
    d3fdgbgbab = 0.0;
    d3fdgbabab = 0.0;
    d3fdababab = 0.0;
    d4fdrararara = 0.004235*t1*(t29*t30*t6*(192.0*t24*t85*
        gradab+16.0*t149*t8*gradab-72.0*t56*t70*gradab-2.0*t274*t5*
        gradab-192.0*t15*t17*t20*t85-192.0*t10*t14*t18*t85-48.0*t10*
        t14*t63*t8-48.0*t15*t17*t60*t8-16.0*t141*t15*t20*t8-16.0*t10*
        t138*t18*t8+72.0*t15*t56*t74+72.0*t10*t56*t73+6.0*t10*t5*pow(t63,
        2.0)+72.0*t10*t18*t56*t63+6.0*t15*t5*pow(t60,2.0)+72.0*t15*
        t20*t56*t60+2.0*t10*t18*t268*t5+2.0*t15*t20*t265*t5+8.0*t141*
        t15*t17*t5+8.0*t10*t138*t14*t5+t278+t277+t276)-0.458*t2*t3*
        t7*(192.0*t85*gradab+12.0*t10*t63*t8-12.0*t15*t60*t8+24.0*
        t15*t17*t56-24.0*t10*t14*t56+4.0*t141*t15*t5-4.0*t10*t138*
        t5+t270+t269+t10*t11*t268-1.0*t11*t15*t265+t262)+6.0*t29*t6*
        t77*t80+6.0*t29*t39*t70*t80-4.0*t36*t40*t6*t80+2.666666666666667*
        t3*t39*t6*t80-4.0*t24*t39*t40*t80+12.0*t24*t29*t36*t80+6.0*
        t29*t30*t70*t77-4.0*t38*t40*t6*t77+2.666666666666667*t3*t30*
        t6*t77-4.0*t24*t30*t40*t77+12.0*t24*t29*t38*t77-0.610666666666667*
        t2*t26*t3*t68*pow(t70,2.0)-4.0*t38*t39*t40*t70-4.0*t30*t36*
        t40*t70+2.666666666666667*t3*t30*t39*t70+12.0*t29*t36*t38*
        t70-21.37333333333334*t2*t53*t7*t66+12.824*t27*t7*t66-1.221333333333333*
        t2*t3*t67*t68*t66+0.916*t2*t3*t70*t25*t66+1.832*t3*t24*t25*
        t66-4.274666666666667*t2*t27*t24*t25*t66-1.333333333333333*
        t157*t39*t40*t6-1.333333333333333*t154*t30*t40*t6-4.148148148148148*
        t27*t38*t39*t6+5.333333333333333*t3*t36*t38*t6+4.0*t154*t29*
        t38*t6-4.148148148148148*t27*t30*t36*t6+4.0*t157*t29*t36*t6-
        1.333333333333333*t149*t30*t39*t40-8.0*t24*t36*t38*t40+t29*
        (t280-19.2*t2*t85+t279)*t6*t39+5.333333333333333*t24*t3*t38*
        t39+4.0*t149*t29*t38*t39+t29*t274*t30*t39-4.148148148148148*
        t24*t27*t30*t39+4.0*t157*t24*t29*t39+5.333333333333333*t24*
        t3*t30*t36+4.0*t149*t29*t30*t36+4.0*t154*t24*t29*t30-1.58320987654321*
        t2*pow(t24,4.0)*t26*t271*t3+t281+t275+61.74518518518519*t106*
        t7*t26-9.499259259259262*t2*t53*t67*t68*t26+5.699555555555556*
        t27*t67*t68*t26-2.442666666666667*t3*t70*t24*t68*t26+5.699555555555556*
        t2*t27*t70*t24*t68*t26-0.814222222222222*t2*t3*t149*t24*t68*
        t26+7.124444444444446*t2*t53*t70*t25*t26-4.274666666666668*
        t27*t70*t25*t26+0.152666666666667*t2*t3*t274*t25*t26+14.24888888888889*
        t53*t24*t25*t26-20.58172839506173*t2*t106*t24*t25*t26+0.610666666666667*
        t3*t149*t25*t26-1.424888888888889*t2*t27*t149*t25*t26+2.849777777777778*
        t2*t3*t70*t67*t146*t26+1.899851851851852*t3*t145*t146*t26-
        4.432987654320988*t2*t27*t145*t146*t26-42.74666666666668*t53*
        t7*t22+61.74518518518519*t2*t106*t7*t22-2.442666666666667*
        t3*t67*t68*t22+5.699555555555556*t2*t27*t67*t68*t22-2.442666666666667*
        t2*t3*t70*t24*t68*t22+1.832*t3*t70*t25*t22-4.274666666666668*
        t2*t27*t70*t25*t22+14.24888888888889*t2*t53*t24*t25*t22-8.549333333333335*
        t27*t24*t25*t22+0.610666666666667*t2*t3*t149*t25*t22+1.899851851851852*
        t2*t3*t145*t146*t22-1.832*t3*t7*t144+4.274666666666668*t2*
        t27*t7*t144+0.610666666666667*t2*t3*t24*t25*t144)+0.01694*
        t159;
    d4fdrarararb = 0.004235*t1*(t29*t30*t6*(48.0*t47*t85*
        gradab+144.0*t24*t85*gradab+12.0*t169*t8*gradab+4.0*t149*t8*
        gradab-36.0*t56*t70*gradab-2.0*t290*t5*gradab-12.0*t10*t2*
        t63*t85+12.0*t15*t2*t60*t85-48.0*t10*t18*t45*t85-48.0*t15*
        t20*t44*t85-144.0*t15*t17*t20*t85-144.0*t10*t14*t18*t85-12.0*
        t10*t45*t63*t8-12.0*t10*t14*t63*t8-12.0*t15*t44*t60*t8-12.0*
        t15*t17*t60*t8-12.0*t15*t161*t20*t8-4.0*t141*t15*t20*t8-12.0*
        t10*t160*t18*t8-4.0*t10*t138*t18*t8+36.0*t15*t56*t74+36.0*
        t10*t56*t73+36.0*t10*t18*t56*t63+36.0*t15*t20*t56*t60+2.0*
        t10*t138*t45*t5+2.0*t141*t15*t44*t5+2.0*t10*t18*t285*t5+2.0*
        t15*t20*t283*t5+6.0*t15*t161*t17*t5+6.0*t10*t14*t160*t5+72.0*
        t15*t2*t20*t295-72.0*t10*t18*t2*t295+t294+t293+t292+t278+t277+
        t276-48.0*t15*t17*t174*t2+48.0*t10*t14*t174*t2)-0.458*t2*t3*
        t7*(96.0*t85*gradab+6.0*t10*t63*t8-6.0*t15*t60*t8-6.0*t10*
        t45*t56+6.0*t15*t44*t56+18.0*t15*t17*t56-18.0*t10*t14*t56+
        3.0*t15*t161*t5-3.0*t10*t160*t5-1.0*t10*t138*t5+t10*t11*t285-
        1.0*t11*t15*t283+t270+t269+t262-12.0*t15*t174*t2-12.0*t10*
        t174*t2+t15*t5*t141)+3.0*t29*t6*t77*t87+3.0*t29*t39*t70*t87-
        2.0*t36*t40*t6*t87-2.0*t24*t39*t40*t87+6.0*t24*t29*t36*t87+
        3.0*t29*t6*t80*t86+3.0*t29*t30*t70*t86-2.0*t38*t40*t6*t86-
        2.0*t24*t30*t40*t86+6.0*t24*t29*t38*t86+3.0*t29*t39*t80*t84+
        3.0*t29*t30*t77*t84-2.0*t38*t39*t40*t84-2.0*t30*t36*t40*t84+
        6.0*t29*t36*t38*t84+6.412*t27*t7*t83-0.610666666666667*t2*
        t3*t67*t68*t83+0.458*t2*t3*t70*t25*t83+0.916*t3*t24*t25*t83-
        2.137333333333334*t2*t27*t24*t25*t83-1.0*t40*t48*t6*t80-1.0*
        t36*t40*t6*t80+1.333333333333333*t3*t39*t6*t80+3.0*t24*t29*
        t48*t80-1.0*t39*t40*t47*t80+3.0*t29*t36*t47*t80-1.0*t24*t39*
        t40*t80-1.0*t40*t49*t6*t77-1.0*t38*t40*t6*t77+1.333333333333333*
        t3*t30*t6*t77+3.0*t24*t29*t49*t77-1.0*t30*t40*t47*t77+3.0*
        t29*t38*t47*t77-1.0*t24*t30*t40*t77-1.0*t39*t40*t49*t70+3.0*
        t29*t36*t49*t70-1.0*t30*t40*t48*t70+3.0*t29*t38*t48*t70-1.0*
        t38*t39*t40*t70-1.0*t30*t36*t40*t70+1.333333333333333*t3*t30*
        t39*t70-10.68666666666667*t2*t53*t7*t66-0.610666666666667*
        t2*t3*t24*t47*t68*t66+0.458*t2*t3*t84*t25*t66+0.458*t3*t47*
        t25*t66-1.068666666666667*t2*t27*t47*t25*t66-0.458*t3*t24*
        t25*t66-1.068666666666667*t2*t27*t24*t25*t66-1.037037037037037*
        t27*t39*t49*t6-1.037037037037037*t27*t30*t48*t6-1.0*t184*t39*
        t40*t6-0.333333333333333*t157*t39*t40*t6-1.0*t175*t30*t40*
        t6-0.333333333333333*t154*t30*t40*t6-3.111111111111111*t27*
        t38*t39*t6+2.666666666666667*t3*t36*t38*t6+3.0*t175*t29*t38*
        t6-3.111111111111111*t27*t30*t36*t6+3.0*t184*t29*t36*t6-2.0*
        t24*t36*t40*t49+t29*t157*t6*t48-2.0*t24*t38*t40*t48+t29*t149*
        t30*t48-2.0*t36*t38*t40*t47-1.037037037037037*t27*t30*t39*
        t47-10.68666666666667*t53*t7*t46+15.4362962962963*t2*t106*
        t7*t46-0.610666666666667*t3*t67*t68*t46+1.424888888888889*
        t2*t27*t67*t68*t46-0.610666666666667*t2*t3*t70*t24*t68*t46+
        0.458*t3*t70*t25*t46-1.068666666666667*t2*t27*t70*t25*t46-
        2.137333333333334*t27*t24*t25*t46+0.152666666666667*t2*t3*
        t149*t25*t46+0.474962962962963*t2*t3*t145*t146*t46-1.0*t169*
        t30*t39*t40-0.333333333333333*t149*t30*t39*t40-2.0*t24*t36*
        t38*t40+t29*(t279-9.600000000000001*t2*t85)*t6*t39+t29*t149*
        t49*t39+t29*t157*t47*t39+2.666666666666667*t24*t3*t38*t39+
        3.0*t169*t29*t38*t39+t29*t290*t30*t39-3.111111111111111*t24*
        t27*t30*t39+3.0*t184*t24*t29*t39+2.666666666666667*t24*t3*
        t30*t36+3.0*t169*t29*t30*t36+t304+t303+t302+t301+t300+3.0*
        t175*t24*t29*t30+t299+t298+t297+t296+t291+t289+t288+t287+t286+
        t281+t275+30.8725925925926*t106*t7*t26-0.610666666666667*t2*
        t3*t70*t84*t68*t26-4.74962962962963*t2*t53*t67*t68*t26-0.610666666666667*
        t3*t70*t47*t68*t26+1.424888888888889*t2*t27*t70*t47*t68*t26+
        2.849777777777778*t27*t24*t47*t68*t26-0.203555555555556*t2*
        t3*t149*t47*t68*t26-1.221333333333333*t3*t84*t24*t68*t26+2.849777777777778*
        t2*t27*t84*t24*t68*t26+0.610666666666667*t3*t70*t24*t68*t26+
        1.424888888888889*t2*t27*t70*t24*t68*t26-0.610666666666667*
        t2*t3*t169*t24*t68*t26-1.58320987654321*t2*t3*t145*t47*t271*
        t26-2.137333333333334*t27*t84*t25*t26+3.562222222222223*t2*
        t53*t70*t25*t26+3.562222222222223*t53*t47*t25*t26-5.145432098765433*
        t2*t106*t47*t25*t26+0.152666666666667*t2*t3*t290*t25*t26+3.562222222222223*
        t53*t24*t25*t26-15.4362962962963*t2*t106*t24*t25*t26+0.458*
        t3*t169*t25*t26-1.068666666666667*t2*t27*t169*t25*t26-0.152666666666667*
        t3*t149*t25*t26-0.356222222222222*t2*t27*t149*t25*t26+1.424888888888889*
        t2*t3*t84*t67*t146*t26+1.424888888888889*t3*t67*t47*t146*t26-
        3.324740740740741*t2*t27*t67*t47*t146*t26+1.424888888888889*
        t2*t3*t70*t24*t47*t146*t26-0.474962962962963*t3*t145*t146*
        t26-1.108246913580247*t2*t27*t145*t146*t26-10.68666666666667*
        t53*t7*t22+46.3088888888889*t2*t106*t7*t22+0.610666666666667*
        t3*t67*t68*t22+1.424888888888889*t2*t27*t67*t68*t22-0.610666666666667*
        t2*t3*t70*t47*t68*t22-1.221333333333333*t3*t24*t47*t68*t22+
        2.849777777777778*t2*t27*t24*t47*t68*t22-1.221333333333333*
        t2*t3*t84*t24*t68*t22+0.916*t3*t84*t25*t22-2.137333333333334*
        t2*t27*t84*t25*t22-0.458*t3*t70*t25*t22-1.068666666666667*
        t2*t27*t70*t25*t22-2.137333333333333*t27*t47*t25*t22+7.124444444444446*
        t2*t53*t24*t25*t22+0.458*t2*t3*t169*t25*t22+1.424888888888889*
        t2*t3*t67*t47*t146*t22-1.374*t3*t7*t162+3.206*t2*t27*t7*t162+
        0.458*t2*t3*t24*t25*t162+t29*t49*t6*t154+t29*t47*t30*t154+
        0.458*t3*t7*t144+1.068666666666667*t2*t27*t7*t144+0.152666666666667*
        t2*t3*t47*t25*t144)+0.012705*t187+0.004235*t159;
    d4fdrararaga = 0.004235*t1*(t29*t30*t6*(-24.0*t73*t8*
        grada-24.0*t18*t63*t8*grada+12.0*t14*t5*t63*grada+72.0*t14*
        t18*t56*grada+4.0*t138*t18*t5*grada+t308)-0.458*t2*t3*t7*(12.0*
        t14*t8*grada-6.0*t5*t63*grada+2.0*t11*t138*grada+t307)-2.0*
        t27*t32*t6*t80*grada+6.0*t24*t3*t32*t80*grada+6.0*t3*t32*t38*
        t70*grada-2.0*t27*t30*t32*t70*grada+2.666666666666667*t32*
        t38*t53*t6*grada+2.0*t157*t3*t32*t6*grada+2.666666666666667*
        t24*t30*t32*t53*grada-4.0*t24*t27*t32*t38*grada+2.0*t149*t3*
        t30*t32*grada+3.0*t29*t6*t80*t95+3.0*t29*t30*t70*t95-2.0*t38*
        t40*t6*t95+1.333333333333333*t3*t30*t6*t95-2.0*t24*t30*t40*
        t95+6.0*t24*t29*t38*t95-10.68666666666667*t2*t53*t7*t93+6.412*
        t27*t7*t93-0.610666666666667*t2*t3*t67*t68*t93+0.458*t2*t3*
        t70*t25*t93+0.916*t3*t24*t25*t93-2.137333333333334*t2*t27*
        t24*t25*t93-1.0*t194*t30*t40*t6+3.0*t194*t29*t38*t6+t306+t305+
        3.0*t194*t24*t29*t30-1.374*t3*t7*t192+3.206*t2*t27*t7*t192+
        0.458*t2*t3*t24*t25*t192-21.37333333333334*grada*t106*t7*t18+
        2.849777777777778*grada*t2*t53*t67*t68*t18-1.221333333333333*
        grada*t27*t67*t68*t18-1.221333333333333*grada*t2*t27*t70*t24*
        t68*t18-2.137333333333334*grada*t2*t53*t70*t25*t18+0.916*grada*
        t27*t70*t25*t18-4.274666666666668*grada*t53*t24*t25*t18+7.124444444444446*
        grada*t2*t106*t24*t25*t18+0.305333333333333*grada*t2*t27*t149*
        t25*t18+0.949925925925926*grada*t2*t27*t145*t146*t18)+0.012705*
        t195;
    d4fdrararagb = 0.004235*t1*(t29*t30*t6*(-24.0*t74*t8*
        gradb-24.0*t20*t60*t8*gradb+12.0*t17*t5*t60*gradb+72.0*t17*
        t20*t56*gradb+4.0*t141*t20*t5*gradb+t312)-0.458*t2*t3*t7*(-
        12.0*t17*t8*gradb+6.0*t5*t60*gradb-2.0*t11*t141*gradb+t311)-
        2.0*t27*t34*t6*t80*gradb+6.0*t24*t3*t34*t80*gradb+6.0*t3*t34*
        t38*t70*gradb-2.0*t27*t30*t34*t70*gradb+2.666666666666667*
        t34*t38*t53*t6*gradb+2.0*t157*t3*t34*t6*gradb+2.666666666666667*
        t24*t30*t34*t53*gradb-4.0*t24*t27*t34*t38*gradb+2.0*t149*t3*
        t30*t34*gradb+3.0*t103*t29*t6*t80+3.0*t103*t29*t30*t70-2.0*
        t103*t38*t40*t6-1.0*t201*t30*t40*t6+3.0*t201*t29*t38*t6+1.333333333333333*
        t103*t3*t30*t6-2.0*t103*t24*t30*t40+6.0*t103*t24*t29*t38+t310+
        t309+3.0*t201*t24*t29*t30+21.37333333333334*gradb*t106*t7*
        t20-2.849777777777778*gradb*t2*t53*t67*t68*t20+1.221333333333333*
        gradb*t27*t67*t68*t20+1.221333333333333*gradb*t2*t27*t70*t24*
        t68*t20+2.137333333333334*gradb*t2*t53*t70*t25*t20-0.916*gradb*
        t27*t70*t25*t20+4.274666666666668*gradb*t53*t24*t25*t20-7.124444444444446*
        gradb*t2*t106*t24*t25*t20-0.305333333333333*gradb*t2*t27*t149*
        t25*t20-0.949925925925926*gradb*t2*t27*t145*t146*t20-1.374*
        t3*t7*t199+3.206*t2*t27*t7*t199+0.458*t2*t3*t24*t25*t199-10.68666666666667*
        t2*t53*t7*t101+6.412*t27*t7*t101-0.610666666666667*t2*t3*t67*
        t68*t101+0.458*t2*t3*t70*t25*t101+0.916*t3*t24*t25*t101-2.137333333333334*
        t2*t27*t24*t25*t101)+0.012705*t202;
    d4fdrararaab = 0.004235*t1*(-6.0*t24*t3*t6*t80+3.0*t112*
        t29*t6*t80+2.0*t27*t54*t80+t29*t30*t6*(12.0*t70*t8-36.0*t24*
        t56-2.0*t149*t5+t317)-6.0*t3*t38*t6*t70+2.0*t27*t30*t6*t70+
        3.0*t112*t29*t30*t70-0.458*t2*t3*(t315-36.0*t56)*t7-1.374*
        t3*t205*t7+3.206*t2*t27*t205*t7+21.37333333333334*t2*t203*
        t7-10.68666666666667*t2*t53*t109*t7+6.412*t27*t109*t7+1.221333333333333*
        t2*t53*t67*t68-0.610666666666667*t2*t3*t109*t67*t68-2.849777777777778*
        t4*t106*t67*t68+1.221333333333333*t4*t53*t70*t24*t68-2.666666666666667*
        t24*t30*t53*t6-2.0*t112*t38*t40*t6-1.0*t208*t30*t40*t6+3.0*
        t208*t29*t38*t6+4.0*t24*t27*t38*t6-2.0*t149*t3*t30*t6+1.333333333333333*
        t112*t3*t30*t6-2.666666666666667*t38*t53*t54-2.0*t157*t3*t54-
        2.0*t112*t24*t30*t40+6.0*t112*t24*t29*t38+t316+t314+3.0*t208*
        t24*t29*t30-0.916*t2*t53*t70*t25+0.458*t2*t3*t109*t70*t25+
        2.137333333333334*t4*t106*t70*t25+0.458*t2*t3*t205*t24*t25-
        7.124444444444446*t4*t203*t24*t25+0.916*t3*t109*t24*t25-2.137333333333334*
        t2*t27*t109*t24*t25+4.274666666666668*t2*t106*t24*t25-0.305333333333333*
        t4*t53*t149*t25-0.949925925925926*t4*t53*t145*t146)+0.012705*
        t209;
    d4fdrararbrb = 0.004235*t1*(t29*t30*t6*(96.0*t47*t85*
        gradab+96.0*t24*t85*gradab-48.0*t56*t84*gradab+8.0*t211*t8*
        gradab+8.0*t169*t8*gradab-12.0*t56*t70*gradab-12.0*t118*t56*
        gradab-2.0*t318*t5*gradab-96.0*t10*t18*t45*t85-96.0*t15*t20*
        t44*t85-96.0*t15*t17*t20*t85-96.0*t10*t14*t18*t85-8.0*t15*
        t17*t85+8.0*t10*t14*t85-8.0*t10*t45*t63*t8-8.0*t15*t44*t60*
        t8-8.0*t15*t161*t20*t8-8.0*t10*t160*t18*t8-8.0*t115*t15*t17*
        t8-8.0*t10*t114*t14*t8+12.0*t15*t56*t74+12.0*t10*t56*t73+12.0*
        t10*t18*t56*t63+2.0*t10*t114*t5*t63+12.0*t15*t20*t56*t60+2.0*
        t115*t15*t5*t60+48.0*t10*t14*t45*t56+48.0*t15*t17*t44*t56+
        12.0*t115*t15*t20*t56+12.0*t10*t114*t18*t56+12.0*t119*t15*
        t56+12.0*t10*t120*t56+4.0*t10*t160*t45*t5+4.0*t15*t161*t44*
        t5+32.0*t10*t174*t2*t45-32.0*t15*t174*t2*t44+16.0*t15*t319*
        t4+16.0*t10*t319*t4+t321+t320+t278+t277+t276+16.0*t15*t174*
        t20-56.0*t15*t17*t174*t2+56.0*t10*t14*t174*t2-16.0*t10*t174*
        t18)+4.0*t29*t6*t86*t87+4.0*t29*t39*t84*t87-1.333333333333333*
        t40*t48*t6*t87-1.333333333333333*t36*t40*t6*t87+1.777777777777778*
        t3*t39*t6*t87+4.0*t24*t29*t48*t87-1.333333333333333*t39*t40*
        t47*t87+4.0*t29*t36*t47*t87-1.333333333333333*t24*t39*t40*
        t87+4.0*t29*t30*t84*t86-1.333333333333333*t40*t49*t6*t86-1.333333333333333*
        t38*t40*t6*t86+1.777777777777778*t3*t30*t6*t86+4.0*t24*t29*
        t49*t86-1.333333333333333*t30*t40*t47*t86+4.0*t29*t38*t47*
        t86-1.333333333333333*t24*t30*t40*t86-0.458*t2*t3*t7*(-4.0*
        t15*t85-4.0*t10*t85+2.0*t10*t63*t8-2.0*t15*t60*t8-2.0*t115*
        t15*t8+2.0*t10*t114*t8-12.0*t10*t45*t56+12.0*t15*t44*t56+12.0*
        t15*t17*t56-12.0*t10*t14*t56+2.0*t15*t161*t5-2.0*t10*t160*
        t5+t270+t269+t262-52.0*t15*t174*t2-52.0*t10*t174*t2)-0.407111111111111*
        t2*t26*t3*t68*pow(t84,2.0)-1.333333333333333*t39*t40*t49*t84+
        4.0*t29*t36*t49*t84-1.333333333333333*t30*t40*t48*t84+4.0*
        t29*t38*t48*t84-1.333333333333333*t38*t39*t40*t84-1.333333333333333*
        t30*t36*t40*t84+1.777777777777778*t3*t30*t39*t84-14.24888888888889*
        t2*t53*t7*t83-0.814222222222222*t2*t3*t24*t47*t68*t83+0.610666666666667*
        t2*t3*t84*t25*t83+0.610666666666667*t3*t47*t25*t83-1.424888888888889*
        t2*t27*t47*t25*t83-0.610666666666667*t3*t24*t25*t83-1.424888888888889*
        t2*t27*t24*t25*t83-0.666666666666667*t40*t48*t6*t80+0.444444444444444*
        t3*t39*t6*t80+2.0*t29*t47*t48*t80-0.666666666666667*t39*t40*
        t47*t80-0.666666666666667*t40*t49*t6*t77+0.444444444444444*
        t3*t30*t6*t77+t29*t122*t6*t77+2.0*t29*t47*t49*t77-0.666666666666667*
        t30*t40*t47*t77+t29*t118*t30*t77+2.0*t29*t48*t49*t70-0.666666666666667*
        t39*t40*t49*t70-0.666666666666667*t30*t40*t48*t70+0.444444444444444*
        t3*t30*t39*t70-3.562222222222223*t2*t53*t7*t66-2.137333333333334*
        t27*t7*t66-0.203555555555556*t2*t3*t117*t68*t66-0.305333333333333*
        t3*t47*t25*t66-0.712444444444445*t2*t27*t47*t25*t66+0.152666666666667*
        t2*t3*t118*t25*t66+0.888888888888889*t3*t48*t49*t6-2.074074074074074*
        t27*t39*t49*t6+1.777777777777778*t3*t36*t49*t6+2.0*t175*t29*
        t49*t6+1.777777777777778*t3*t38*t48*t6-2.074074074074074*t27*
        t30*t48*t6+2.0*t184*t29*t48*t6-0.666666666666667*t213*t39*
        t40*t6-0.666666666666667*t184*t39*t40*t6-0.666666666666667*
        t121*t38*t40*t6-0.666666666666667*t122*t36*t40*t6-0.666666666666667*
        t212*t30*t40*t6-0.666666666666667*t175*t30*t40*t6-2.074074074074074*
        t27*t38*t39*t6+0.444444444444444*t122*t3*t39*t6+0.888888888888889*
        t3*t36*t38*t6+2.0*t212*t29*t38*t6-2.074074074074074*t27*t30*
        t36*t6+2.0*t213*t29*t36*t6+0.444444444444444*t121*t3*t30*t6-
        1.333333333333333*t24*t40*t48*t49-1.333333333333333*t36*t40*
        t47*t49+0.888888888888889*t3*t39*t47*t49-1.333333333333333*
        t24*t36*t40*t49+1.777777777777778*t24*t3*t39*t49+2.0*t169*
        t29*t39*t49-1.333333333333333*t38*t40*t47*t48+0.888888888888889*
        t3*t30*t47*t48-1.333333333333333*t24*t38*t40*t48+1.777777777777778*
        t24*t3*t30*t48+2.0*t169*t29*t30*t48-1.333333333333333*t36*
        t38*t40*t47+1.777777777777778*t3*t38*t39*t47-2.074074074074074*
        t27*t30*t39*t47+2.0*t184*t29*t39*t47+1.777777777777778*t3*
        t30*t36*t47+2.0*t175*t29*t30*t47-7.124444444444446*t53*t7*
        t46+30.8725925925926*t2*t106*t7*t46+0.407111111111111*t3*t67*
        t68*t46+0.949925925925926*t2*t27*t67*t68*t46-0.407111111111111*
        t2*t3*t70*t47*t68*t46-0.814222222222222*t3*t24*t47*t68*t46+
        1.899851851851852*t2*t27*t24*t47*t68*t46-0.814222222222222*
        t2*t3*t84*t24*t68*t46+0.610666666666667*t3*t84*t25*t46-1.424888888888889*
        t2*t27*t84*t25*t46-0.305333333333333*t3*t70*t25*t46-0.712444444444445*
        t2*t27*t70*t25*t46+2.374814814814815*t2*t53*t47*t25*t46-1.424888888888889*
        t27*t47*t25*t46+4.74962962962963*t2*t53*t24*t25*t46+0.305333333333333*
        t2*t3*t169*t25*t46+0.949925925925926*t2*t3*t67*t47*t146*t46-
        0.666666666666667*t118*t38*t39*t40-0.666666666666667*t211*
        t30*t39*t40-0.666666666666667*t169*t30*t39*t40-0.666666666666667*
        t122*t24*t39*t40-0.666666666666667*t118*t30*t36*t40-0.666666666666667*
        t121*t24*t30*t40+t29*t118*t80*t39+t29*(t279-2.4*t56)*t6*t39+
        0.888888888888889*t24*t3*t38*t39+2.0*t211*t29*t38*t39+t29*
        t318*t30*t39+0.444444444444444*t118*t3*t30*t39-2.074074074074074*
        t24*t27*t30*t39+2.0*t213*t24*t29*t39+t29*t70*t122*t39+2.0*
        t118*t29*t36*t38+2.0*t121*t24*t29*t38+0.888888888888889*t24*
        t3*t30*t36+2.0*t211*t29*t30*t36+2.0*t122*t24*t29*t36+2.0*t212*
        t24*t29*t30+t281+t275-0.203555555555556*t2*t3*t118*t70*t68*
        t26-1.58320987654321*t2*t53*t67*t68*t26-0.949925925925926*
        t27*t67*t68*t26-0.814222222222222*t3*t84*t47*t68*t26+1.899851851851852*
        t2*t27*t84*t47*t68*t26+0.407111111111111*t3*t70*t47*t68*t26+
        0.949925925925926*t2*t27*t70*t47*t68*t26-6.33283950617284*
        t2*t53*t24*t47*t68*t26-1.110223024625157E-16*t27*t24*t47*t68*
        t26-0.407111111111111*t2*t3*t169*t47*t68*t26+0.814222222222222*
        t3*t84*t24*t68*t26+1.899851851851852*t2*t27*t84*t24*t68*t26-
        0.407111111111111*t2*t3*t211*t24*t68*t26-0.407111111111111*
        t3*t118*t24*t68*t26+0.949925925925926*t2*t27*t118*t24*t68*
        t26-1.58320987654321*t2*t53*t117*t68*t26+0.949925925925926*
        t27*t117*t68*t26-1.58320987654321*t2*t3*t67*t117*t271*t26+
        4.74962962962963*t2*t53*t84*t25*t26+1.187407407407408*t2*t53*
        t70*t25*t26+0.712444444444445*t27*t70*t25*t26+2.374814814814815*
        t53*t47*t25*t26-10.29086419753087*t2*t106*t47*t25*t26+0.152666666666667*
        t2*t3*t318*t25*t26-2.374814814814815*t53*t24*t25*t26-10.29086419753087*
        t2*t106*t24*t25*t26+0.305333333333333*t3*t211*t25*t26-0.712444444444445*
        t2*t27*t211*t25*t26-0.305333333333333*t3*t169*t25*t26-0.712444444444445*
        t2*t27*t169*t25*t26+1.187407407407408*t2*t53*t118*t25*t26-
        0.712444444444445*t27*t118*t25*t26+0.474962962962963*t2*t3*
        t118*t67*t146*t26-0.949925925925926*t3*t67*t47*t146*t26-2.216493827160494*
        t2*t27*t67*t47*t146*t26+1.899851851851852*t2*t3*t84*t24*t47*
        t146*t26+0.474962962962963*t2*t3*t70*t117*t146*t26+0.949925925925926*
        t3*t24*t117*t146*t26-2.216493827160494*t2*t27*t24*t117*t146*
        t26+7.124444444444446*t53*t7*t22+30.8725925925926*t2*t106*
        t7*t22-0.814222222222222*t2*t3*t84*t47*t68*t22+0.814222222222222*
        t3*t24*t47*t68*t22+1.899851851851852*t2*t27*t24*t47*t68*t22-
        0.407111111111111*t2*t3*t118*t24*t68*t22-0.407111111111111*
        t3*t117*t68*t22+0.949925925925926*t2*t27*t117*t68*t22-0.610666666666667*
        t3*t84*t25*t22-1.424888888888889*t2*t27*t84*t25*t22+4.74962962962963*
        t2*t53*t47*t25*t22+2.374814814814815*t2*t53*t24*t25*t22+1.424888888888889*
        t27*t24*t25*t22+0.305333333333333*t2*t3*t211*t25*t22+0.305333333333333*
        t3*t118*t25*t22-0.712444444444445*t2*t27*t118*t25*t22+0.949925925925926*
        t2*t3*t24*t117*t146*t22-0.916*t3*t7*t210+2.137333333333334*
        t2*t27*t7*t210+0.305333333333333*t2*t3*t24*t25*t210+0.916*
        t3*t7*t162+2.137333333333334*t2*t27*t7*t162+0.305333333333333*
        t2*t3*t47*t25*t162+t29*t80*t6*t121+t29*t70*t30*t121-3.562222222222223*
        t2*t53*t7*t116+2.137333333333334*t27*t7*t116-0.203555555555556*
        t2*t3*t67*t68*t116+0.152666666666667*t2*t3*t70*t25*t116+0.305333333333333*
        t3*t24*t25*t116-0.712444444444445*t2*t27*t24*t25*t116)+0.00847*
        t214+0.00847*t187;
    d4fdrararbga = 0.004235*t1*(t29*t30*t6*(-16.0*t14*t2*
        t85*grada-8.0*t73*t8*grada-8.0*t18*t63*t8*grada+4.0*t45*t5*
        t63*grada+24.0*t18*t45*t56*grada+48.0*t14*t18*t56*grada+4.0*
        t160*t18*t5*grada+32.0*t174*t18*t2*grada+t331+t308)-0.458*
        t2*t3*t7*(8.0*t2*t85*grada+4.0*t45*t8*grada+8.0*t14*t8*grada-
        2.0*t5*t63*grada+2.0*t11*t160*grada+t307)+4.0*t24*t3*t32*t87*
        grada+4.0*t3*t32*t38*t84*grada-0.666666666666667*t27*t32*t6*
        t80*grada+2.0*t3*t32*t47*t80*grada+2.0*t3*t32*t49*t70*grada-
        0.666666666666667*t27*t30*t32*t70*grada+0.888888888888889*
        t32*t49*t53*t6*grada+1.777777777777778*t32*t38*t53*t6*grada+
        2.0*t184*t3*t32*t6*grada+0.888888888888889*t30*t32*t47*t53*
        grada+1.777777777777778*t24*t30*t32*t53*grada-1.333333333333333*
        t24*t27*t32*t38*grada+2.0*t169*t3*t30*t32*grada+2.0*t29*t6*
        t87*t95+2.0*t29*t30*t84*t95-0.666666666666667*t38*t40*t6*t95+
        0.888888888888889*t3*t30*t6*t95+2.0*t24*t29*t49*t95+2.0*t29*
        t38*t47*t95-0.666666666666667*t24*t30*t40*t95-7.124444444444446*
        t2*t53*t7*t93-0.407111111111111*t2*t3*t24*t47*t68*t93+0.305333333333333*
        t2*t3*t84*t25*t93+0.305333333333333*t3*t47*t25*t93-0.305333333333333*
        t3*t24*t25*t93-0.712444444444445*t2*t27*t24*t25*t93-0.333333333333333*
        t194*t30*t40*t6+2.0*t216*t29*t38*t6+0.444444444444444*t125*
        t3*t30*t6+2.0*t125*t24*t29*t38+t336+t335+t334+t333+t332+t330+
        t329+t328+t327+t326+t325+t324+t323+t322+t306+t305+2.0*t216*
        t24*t29*t30-0.916*t3*t7*t215+0.305333333333333*t2*t3*t24*t25*
        t215+t29*t49*t6*t194+t29*t47*t30*t194+0.458*t3*t7*t192+1.068666666666667*
        t2*t27*t7*t192+0.152666666666667*t2*t3*t47*t25*t192-7.124444444444446*
        grada*t106*t7*t18+0.949925925925926*grada*t2*t53*t67*t68*t18+
        0.407111111111111*grada*t27*t67*t68*t18-0.407111111111111*
        grada*t2*t27*t70*t47*t68*t18-0.814222222222222*grada*t27*t24*
        t47*t68*t18-0.814222222222222*grada*t2*t27*t84*t24*t68*t18+
        0.610666666666667*grada*t27*t84*t25*t18-0.712444444444445*
        grada*t2*t53*t70*t25*t18-0.305333333333333*grada*t27*t70*t25*
        t18-1.424888888888889*grada*t53*t47*t25*t18+2.374814814814815*
        grada*t2*t106*t47*t25*t18+4.74962962962963*grada*t2*t106*t24*
        t25*t18+0.305333333333333*grada*t2*t27*t169*t25*t18+0.949925925925926*
        grada*t2*t27*t67*t47*t146*t18+t29*t80*t6*t125+t29*t70*t30*
        t125-3.562222222222223*t2*t53*t7*t124+2.137333333333334*t27*
        t7*t124-0.203555555555556*t2*t3*t67*t68*t124+0.152666666666667*
        t2*t3*t70*t25*t124+0.305333333333333*t3*t24*t25*t124)+t337+
        0.004235*t195;
    d4fdrararbgb = 0.004235*t1*(t29*t30*t6*(16.0*t17*t2*t85*
        gradb-8.0*t74*t8*gradb-8.0*t20*t60*t8*gradb+4.0*t44*t5*t60*
        gradb+24.0*t20*t44*t56*gradb+48.0*t17*t20*t56*gradb+4.0*t161*
        t20*t5*gradb-32.0*t174*t2*t20*gradb+t347+t312)-0.458*t2*t3*
        t7*(8.0*t2*t85*gradb-4.0*t44*t8*gradb-8.0*t17*t8*gradb+2.0*
        t5*t60*gradb-2.0*t11*t161*gradb+t311)+4.0*t24*t3*t34*t87*gradb+
        4.0*t3*t34*t38*t84*gradb-0.666666666666667*t27*t34*t6*t80*
        gradb+2.0*t3*t34*t47*t80*gradb+2.0*t3*t34*t49*t70*gradb-0.666666666666667*
        t27*t30*t34*t70*gradb+0.888888888888889*t34*t49*t53*t6*gradb+
        1.777777777777778*t34*t38*t53*t6*gradb+2.0*t184*t3*t34*t6*
        gradb+0.888888888888889*t30*t34*t47*t53*gradb+1.777777777777778*
        t24*t30*t34*t53*gradb-1.333333333333333*t24*t27*t34*t38*gradb+
        2.0*t169*t3*t30*t34*gradb+2.0*t103*t29*t6*t87+2.0*t103*t29*
        t30*t84-0.666666666666667*t103*t38*t40*t6-0.333333333333333*
        t201*t30*t40*t6+2.0*t219*t29*t38*t6+0.444444444444444*t128*
        t3*t30*t6+0.888888888888889*t103*t3*t30*t6+2.0*t103*t24*t29*
        t49+2.0*t103*t29*t38*t47-0.666666666666667*t103*t24*t30*t40+
        2.0*t128*t24*t29*t38+t352+t351+t350+t349+t348+t346+t345+t344+
        t343+t342+t341+t340+t339+t338+t310+t309+2.0*t219*t24*t29*t30-
        0.916*t3*t7*t218+0.305333333333333*t2*t3*t24*t25*t218+t29*
        t49*t6*t201+t29*t47*t30*t201+7.124444444444446*gradb*t106*
        t7*t20-0.949925925925926*gradb*t2*t53*t67*t68*t20-0.407111111111111*
        gradb*t27*t67*t68*t20+0.407111111111111*gradb*t2*t27*t70*t47*
        t68*t20+0.814222222222222*gradb*t27*t24*t47*t68*t20+0.814222222222222*
        gradb*t2*t27*t84*t24*t68*t20-0.610666666666667*gradb*t27*t84*
        t25*t20+0.712444444444445*gradb*t2*t53*t70*t25*t20+0.305333333333333*
        gradb*t27*t70*t25*t20+1.424888888888889*gradb*t53*t47*t25*
        t20-2.374814814814815*gradb*t2*t106*t47*t25*t20-4.74962962962963*
        gradb*t2*t106*t24*t25*t20-0.305333333333333*gradb*t2*t27*t169*
        t25*t20-0.949925925925926*gradb*t2*t27*t67*t47*t146*t20+0.458*
        t3*t7*t199+1.068666666666667*t2*t27*t7*t199+0.152666666666667*
        t2*t3*t47*t25*t199+t29*t80*t6*t128+t29*t70*t30*t128-3.562222222222223*
        t2*t53*t7*t127+2.137333333333334*t27*t7*t127-0.203555555555556*
        t2*t3*t67*t68*t127+0.152666666666667*t2*t3*t70*t25*t127+0.305333333333333*
        t3*t24*t25*t127-7.124444444444446*t2*t53*t7*t101-0.407111111111111*
        t2*t3*t24*t47*t68*t101+0.305333333333333*t2*t3*t84*t25*t101+
        0.305333333333333*t3*t47*t25*t101-0.305333333333333*t3*t24*
        t25*t101-0.712444444444445*t2*t27*t24*t25*t101)+t353+0.004235*
        t202;
    d4fdrararbab = 0.004235*t1*(-4.0*t24*t3*t6*t87+2.0*t112*
        t29*t6*t87-4.0*t3*t38*t6*t84+2.0*t112*t29*t30*t84-2.0*t3*t47*
        t6*t80+0.666666666666667*t27*t54*t80+t29*t30*t6*(4.0*t70*t8-
        12.0*t47*t56-24.0*t24*t56-2.0*t169*t5+t363+t317)-2.0*t3*t49*
        t6*t70+0.666666666666667*t27*t30*t6*t70-0.458*t2*t3*(t282+
        t315)*t7+0.458*t3*t205*t7+1.068666666666667*t2*t27*t205*t7+
        18.11644444444445*t2*t203*t7-3.562222222222223*t2*t53*t130*
        t7+2.137333333333334*t27*t130*t7-7.124444444444446*t2*t53*
        t109*t7-0.407111111111111*t2*t53*t67*t68-0.203555555555556*
        t2*t3*t130*t67*t68-0.949925925925926*t4*t106*t67*t68+0.407111111111111*
        t4*t53*t70*t47*t68+0.814222222222222*t2*t53*t24*t47*t68-0.407111111111111*
        t2*t3*t109*t24*t47*t68+0.814222222222222*t4*t53*t84*t24*t68-
        0.888888888888889*t30*t47*t53*t6-1.777777777777778*t24*t30*
        t53*t6-0.666666666666667*t112*t38*t40*t6-0.333333333333333*
        t208*t30*t40*t6+2.0*t221*t29*t38*t6+1.333333333333333*t24*
        t27*t38*t6-2.0*t169*t3*t30*t6+0.444444444444444*t131*t3*t30*
        t6+0.888888888888889*t112*t3*t30*t6-0.888888888888889*t49*
        t53*t54-1.777777777777778*t38*t53*t54-2.0*t184*t3*t54+2.0*
        t112*t24*t29*t49+2.0*t112*t29*t38*t47-0.666666666666667*t112*
        t24*t30*t40+2.0*t131*t24*t29*t38+t368+t367+t366+t365+t364+
        t362+t361+t360+t359+t358+t357+t356+t355+t354+t316+2.0*t221*
        t24*t29*t30-0.610666666666667*t2*t53*t84*t25+0.305333333333333*
        t2*t3*t109*t84*t25+0.305333333333333*t2*t53*t70*t25+0.712444444444445*
        t4*t106*t70*t25+0.152666666666667*t2*t3*t205*t47*t25-2.374814814814815*
        t4*t203*t47*t25+0.305333333333333*t3*t109*t47*t25+1.424888888888889*
        t2*t106*t47*t25-8.41362962962963*t4*t203*t24*t25+0.305333333333333*
        t3*t130*t24*t25-0.305333333333333*t3*t109*t24*t25-0.712444444444445*
        t2*t27*t109*t24*t25-0.305333333333333*t4*t53*t169*t25+0.152666666666667*
        t2*t3*t70*t130*t25+t29*t49*t6*t208+t29*t47*t30*t208-0.949925925925926*
        t4*t53*t67*t47*t146+t29*t80*t6*t131+t29*t70*t30*t131)+t369+
        0.004235*t209;
    d4fdraragaga = 0.004235*t1*(2.0*t3*t32*t6*t80+t29*t30*
        t6*(-16.0*t14*t18*t8+4.0*t5*t73+4.0*t18*t5*t63+t373)+2.0*t3*
        t30*t32*t70-0.458*t2*t3*(2.0*t11*t63-4.0*t14*t5+t372)*t7-0.666666666666667*
        t229*t30*t40*t6-1.333333333333333*t27*t32*t38*t6+2.0*t229*
        t29*t38*t6+4.0*t24*t3*t32*t38+t371+t370-1.333333333333333*
        t24*t27*t30*t32+2.0*t229*t24*t29*t30-0.916*t3*t7*t227+2.137333333333334*
        t2*t27*t7*t227+0.305333333333333*t2*t3*t24*t25*t227+4.274666666666668*
        t53*t7*t18-0.407111111111111*t2*t27*t67*t68*t18+0.305333333333333*
        t2*t27*t70*t25*t18-1.424888888888889*t2*t53*t24*t25*t18+0.610666666666667*
        t27*t24*t25*t18)+0.00847*t230;
    d4fdraragagb = 0.0;
    d4fdraragaab = 0.0;
    d4fdraragbgb = 0.004235*t1*(2.0*t3*t34*t6*t80+t29*t30*
        t6*(-16.0*t17*t20*t8+4.0*t5*t74+4.0*t20*t5*t60+t377)+2.0*t3*
        t30*t34*t70-0.458*t2*t3*(-2.0*t11*t60+4.0*t17*t5+t376)*t7-
        0.666666666666667*t237*t30*t40*t6-1.333333333333333*t27*t34*
        t38*t6+2.0*t237*t29*t38*t6+4.0*t24*t3*t34*t38+t375+t374-1.333333333333333*
        t24*t27*t30*t34+2.0*t237*t24*t29*t30-0.916*t3*t7*t235+2.137333333333334*
        t2*t27*t7*t235+0.305333333333333*t2*t3*t24*t25*t235-4.274666666666668*
        t53*t7*t20+0.407111111111111*t2*t27*t67*t68*t20-0.305333333333333*
        t2*t27*t70*t25*t20+1.424888888888889*t2*t53*t24*t25*t20-0.610666666666667*
        t27*t24*t25*t20)+0.00847*t238;
    d4fdraragbab = 0.0;
    d4fdraraabab = 0.0;
    d4fdrarbrbrb = 0.004235*t1*(t29*t30*t6*(144.0*t47*t85*
        gradab+48.0*t24*t85*gradab+4.0*t243*t8*gradab+12.0*t211*t8*
        gradab-36.0*t118*t56*gradab-2.0*t378*t5*gradab-144.0*t10*t18*
        t45*t85+12.0*t10*t45*t85-144.0*t15*t20*t44*t85-12.0*t15*t44*
        t85-48.0*t15*t17*t20*t85+12.0*t115*t15*t2*t85-12.0*t10*t114*
        t2*t85-48.0*t10*t14*t18*t85-12.0*t10*t114*t45*t8-12.0*t115*
        t15*t44*t8-4.0*t10*t18*t240*t8-4.0*t15*t20*t239*t8-12.0*t115*
        t15*t17*t8-12.0*t10*t114*t14*t8+36.0*t115*t15*t20*t56+36.0*
        t10*t114*t18*t56+36.0*t119*t15*t56+36.0*t10*t120*t56+2.0*t10*
        t14*t240*t5+2.0*t15*t17*t239*t5+84.0*t10*t174*t2*t45-84.0*
        t15*t174*t2*t44+t321+t320+t294+t293+t292+t278+t277+t276+48.0*
        t15*t174*t20-48.0*t10*t174*t18)-0.458*t2*t3*t7*(-96.0*t85*
        gradab-18.0*t15*t85-18.0*t10*t85-6.0*t115*t15*t8+6.0*t10*t114*
        t8-18.0*t10*t45*t56+18.0*t15*t44*t56+6.0*t15*t17*t56-6.0*t10*
        t14*t56-1.0*t10*t240*t5+t270+t269+t262+t15*t5*t239-54.0*t15*
        t174*t2-54.0*t10*t174*t2)-2.0*t40*t48*t6*t87+3.0*t121*t29*
        t6*t87+6.0*t29*t47*t48*t87-2.0*t39*t40*t47*t87+3.0*t118*t29*
        t39*t87-2.0*t40*t49*t6*t86+3.0*t122*t29*t6*t86+6.0*t29*t47*
        t49*t86-2.0*t30*t40*t47*t86+3.0*t118*t29*t30*t86+6.0*t29*t48*
        t49*t84-2.0*t39*t40*t49*t84-2.0*t30*t40*t48*t84+3.0*t122*t29*
        t39*t84+3.0*t121*t29*t30*t84-6.412*t27*t7*t83-0.610666666666667*
        t2*t3*t117*t68*t83-0.916*t3*t47*t25*t83-2.137333333333334*
        t2*t27*t47*t25*t83+0.458*t2*t3*t118*t25*t83+2.666666666666667*
        t3*t48*t49*t6-1.0*t121*t40*t49*t6-3.111111111111111*t27*t39*
        t49*t6+3.0*t212*t29*t49*t6-1.0*t122*t40*t48*t6-3.111111111111111*
        t27*t30*t48*t6+3.0*t213*t29*t48*t6-0.333333333333333*t245*
        t39*t40*t6-1.0*t213*t39*t40*t6-1.0*t121*t38*t40*t6-1.0*t122*
        t36*t40*t6-0.333333333333333*t244*t30*t40*t6-1.0*t212*t30*
        t40*t6-1.037037037037037*t27*t38*t39*t6+1.333333333333333*
        t122*t3*t39*t6-1.037037037037037*t27*t30*t36*t6+1.333333333333333*
        t121*t3*t30*t6-2.0*t40*t47*t48*t49-2.0*t24*t40*t48*t49-2.0*
        t36*t40*t47*t49+2.666666666666667*t3*t39*t47*t49-1.0*t118*
        t39*t40*t49+3.0*t211*t29*t39*t49+3.0*t118*t29*t36*t49+3.0*
        t121*t24*t29*t49-2.0*t38*t40*t47*t48+2.666666666666667*t3*
        t30*t47*t48-1.0*t118*t30*t40*t48+3.0*t118*t29*t38*t48+3.0*
        t211*t29*t30*t48+3.0*t122*t24*t29*t48-1.0*t122*t39*t40*t47-
        1.0*t121*t30*t40*t47-3.111111111111111*t27*t30*t39*t47+3.0*
        t213*t29*t39*t47+3.0*t121*t29*t38*t47+3.0*t122*t29*t36*t47+
        3.0*t212*t29*t30*t47+10.68666666666667*t53*t7*t46+46.3088888888889*
        t2*t106*t7*t46-1.221333333333333*t2*t3*t84*t47*t68*t46+1.221333333333333*
        t3*t24*t47*t68*t46+2.849777777777778*t2*t27*t24*t47*t68*t46-
        0.610666666666667*t2*t3*t118*t24*t68*t46-0.610666666666667*
        t3*t117*t68*t46+1.424888888888889*t2*t27*t117*t68*t46-0.916*
        t3*t84*t25*t46-2.137333333333334*t2*t27*t84*t25*t46+7.124444444444446*
        t2*t53*t47*t25*t46+2.137333333333334*t27*t24*t25*t46+0.458*
        t2*t3*t211*t25*t46+0.458*t3*t118*t25*t46-1.068666666666667*
        t2*t27*t118*t25*t46+1.424888888888889*t2*t3*t24*t117*t146*
        t46-1.0*t118*t38*t39*t40-0.333333333333333*t243*t30*t39*t40-
        1.0*t211*t30*t39*t40-1.0*t122*t24*t39*t40-1.0*t118*t30*t36*
        t40-1.0*t121*t24*t30*t40+t29*(9.600000000000001*t2*t85+t279)*
        t6*t39+t29*t243*t38*t39+t29*t378*t30*t39+1.333333333333333*
        t118*t3*t30*t39-1.037037037037037*t24*t27*t30*t39+t29*t245*
        t24*t39+t29*t245*t6*t36+t29*t243*t30*t36+t304+t303+t302+t301+
        t300+t299+t298+t297+t296+t291+t289+t288+t287+t286+t281+t275-
        30.8725925925926*t106*t7*t26-0.610666666666667*t2*t3*t118*
        t84*t68*t26+1.221333333333333*t3*t84*t47*t68*t26+2.849777777777778*
        t2*t27*t84*t47*t68*t26-2.849777777777778*t27*t24*t47*t68*t26-
        0.610666666666667*t2*t3*t211*t47*t68*t26-0.610666666666667*
        t3*t118*t47*t68*t26+1.424888888888889*t2*t27*t118*t47*t68*
        t26-0.203555555555556*t2*t3*t243*t24*t68*t26+0.610666666666667*
        t3*t118*t24*t68*t26+1.424888888888889*t2*t27*t118*t24*t68*
        t26-4.74962962962963*t2*t53*t117*t68*t26-5.55111512312578E-17*
        t27*t117*t68*t26-1.58320987654321*t2*t3*t24*t242*t271*t26+
        2.137333333333334*t27*t84*t25*t26-3.562222222222223*t53*t47*
        t25*t26-15.4362962962963*t2*t106*t47*t25*t26+0.152666666666667*
        t2*t3*t378*t25*t26+0.152666666666667*t3*t243*t25*t26-0.356222222222222*
        t2*t27*t243*t25*t26-3.562222222222223*t53*t24*t25*t26-5.145432098765433*
        t2*t106*t24*t25*t26-0.458*t3*t211*t25*t26-1.068666666666667*
        t2*t27*t211*t25*t26+3.562222222222223*t2*t53*t118*t25*t26+
        1.424888888888889*t2*t3*t118*t24*t47*t146*t26+0.474962962962963*
        t3*t242*t146*t26-1.108246913580247*t2*t27*t242*t146*t26+1.424888888888889*
        t2*t3*t84*t117*t146*t26-1.424888888888889*t3*t24*t117*t146*
        t26-3.324740740740741*t2*t27*t24*t117*t146*t26+t29*t38*t6*
        t244+t29*t24*t30*t244-0.458*t3*t7*t241+1.068666666666667*t2*
        t27*t7*t241+0.152666666666667*t2*t3*t24*t25*t241+10.68666666666667*
        t53*t7*t22+15.4362962962963*t2*t106*t7*t22-0.610666666666667*
        t2*t3*t118*t47*t68*t22+0.610666666666667*t3*t117*t68*t22+1.424888888888889*
        t2*t27*t117*t68*t22+2.137333333333334*t27*t47*t25*t22+0.152666666666667*
        t2*t3*t243*t25*t22-0.458*t3*t118*t25*t22-1.068666666666667*
        t2*t27*t118*t25*t22+0.474962962962963*t2*t3*t242*t146*t22+
        1.374*t3*t7*t210+3.206*t2*t27*t7*t210+0.458*t2*t3*t47*t25*
        t210-10.68666666666667*t2*t53*t7*t116-0.610666666666667*t2*
        t3*t24*t47*t68*t116+0.458*t2*t3*t84*t25*t116+0.458*t3*t47*
        t25*t116-1.068666666666667*t2*t27*t47*t25*t116-0.458*t3*t24*
        t25*t116-1.068666666666667*t2*t27*t24*t25*t116)+0.004235*t246+
        0.012705*t214;
    d4fdrarbrbga = 0.004235*t1*(t29*t30*t6*(-16.0*t2*t45*
        t85*grada+8.0*t18*t85*grada-8.0*t114*t18*t8*grada-8.0*t120*
        t8*grada+48.0*t18*t45*t56*grada+24.0*t14*t18*t56*grada+4.0*
        t114*t14*t5*grada+56.0*t174*t18*t2*grada+t331+t308)-0.458*
        t2*t3*t7*(20.0*t2*t85*grada+8.0*t45*t8*grada+4.0*t14*t8*grada+
        4.0*t56*grada-2.0*t114*t5*grada+t307)+4.0*t3*t32*t47*t87*grada+
        4.0*t3*t32*t49*t84*grada+1.777777777777778*t32*t49*t53*t6*
        grada+0.888888888888889*t32*t38*t53*t6*grada+2.0*t213*t3*t32*
        t6*grada-0.666666666666667*t122*t27*t32*t6*grada+1.777777777777778*
        t30*t32*t47*t53*grada+0.888888888888889*t24*t30*t32*t53*grada-
        1.333333333333333*t27*t32*t47*t49*grada+2.0*t118*t3*t32*t38*
        grada+2.0*t211*t3*t30*t32*grada-0.666666666666667*t118*t27*
        t30*t32*grada+2.0*t122*t24*t3*t32*grada+0.444444444444444*
        t3*t30*t6*t95+t29*t122*t6*t95+2.0*t29*t47*t49*t95+t29*t118*
        t30*t95-3.562222222222223*t2*t53*t7*t93-2.137333333333334*
        t27*t7*t93-0.203555555555556*t2*t3*t117*t68*t93-0.305333333333333*
        t3*t47*t25*t93+0.152666666666667*t2*t3*t118*t25*t93+2.0*t125*
        t29*t6*t87+2.0*t125*t29*t30*t84-0.666666666666667*t125*t40*
        t49*t6+2.0*t216*t29*t49*t6-0.333333333333333*t248*t30*t40*
        t6+0.888888888888889*t125*t3*t30*t6+2.0*t125*t24*t29*t49-0.666666666666667*
        t125*t30*t40*t47+2.0*t125*t29*t38*t47+2.0*t216*t29*t30*t47+
        t336+t335+t334+t333+t332+t330+t329+t328+t327+t326+t325+t324+
        t323+t322+t306+t305+t29*t38*t6*t248+t29*t24*t30*t248-0.458*
        t3*t7*t247+1.068666666666667*t2*t27*t7*t247+0.152666666666667*
        t2*t3*t24*t25*t247+0.916*t3*t7*t215+0.305333333333333*t2*t3*
        t47*t25*t215+7.124444444444446*grada*t106*t7*t18-0.814222222222222*
        grada*t2*t27*t84*t47*t68*t18+0.814222222222222*grada*t27*t24*
        t47*t68*t18-0.407111111111111*grada*t2*t27*t118*t24*t68*t18+
        0.949925925925926*grada*t2*t53*t117*t68*t18-0.407111111111111*
        grada*t27*t117*t68*t18-0.610666666666667*grada*t27*t84*t25*
        t18+4.74962962962963*grada*t2*t106*t47*t25*t18+1.424888888888889*
        grada*t53*t24*t25*t18+2.374814814814815*grada*t2*t106*t24*
        t25*t18+0.305333333333333*grada*t2*t27*t211*t25*t18-0.712444444444445*
        grada*t2*t53*t118*t25*t18+0.305333333333333*grada*t27*t118*
        t25*t18+0.949925925925926*grada*t2*t27*t24*t117*t146*t18-7.124444444444446*
        t2*t53*t7*t124-0.407111111111111*t2*t3*t24*t47*t68*t124+0.305333333333333*
        t2*t3*t84*t25*t124+0.305333333333333*t3*t47*t25*t124-0.712444444444445*
        t2*t27*t47*t25*t124-0.305333333333333*t3*t24*t25*t124)+t337+
        0.004235*t249;
    d4fdrarbrbgb = 0.004235*t1*(t29*t30*t6*(16.0*t2*t44*t85*
        gradb-8.0*t20*t85*gradb-8.0*t115*t20*t8*gradb-8.0*t119*t8*
        gradb+48.0*t20*t44*t56*gradb+24.0*t17*t20*t56*gradb+4.0*t115*
        t17*t5*gradb-56.0*t174*t2*t20*gradb+t347+t312)-0.458*t2*t3*
        t7*(20.0*t2*t85*gradb-8.0*t44*t8*gradb-4.0*t17*t8*gradb+4.0*
        t56*gradb+2.0*t115*t5*gradb+t311)+4.0*t3*t34*t47*t87*gradb+
        4.0*t3*t34*t49*t84*gradb+1.777777777777778*t34*t49*t53*t6*
        gradb+0.888888888888889*t34*t38*t53*t6*gradb+2.0*t213*t3*t34*
        t6*gradb-0.666666666666667*t122*t27*t34*t6*gradb+1.777777777777778*
        t30*t34*t47*t53*gradb+0.888888888888889*t24*t30*t34*t53*gradb-
        1.333333333333333*t27*t34*t47*t49*gradb+2.0*t118*t3*t34*t38*
        gradb+2.0*t211*t3*t30*t34*gradb-0.666666666666667*t118*t27*
        t30*t34*gradb+2.0*t122*t24*t3*t34*gradb+2.0*t128*t29*t6*t87+
        2.0*t128*t29*t30*t84-0.666666666666667*t128*t40*t49*t6+2.0*
        t219*t29*t49*t6-0.333333333333333*t251*t30*t40*t6+0.888888888888889*
        t128*t3*t30*t6+0.444444444444444*t103*t3*t30*t6+2.0*t103*t29*
        t47*t49+2.0*t128*t24*t29*t49-0.666666666666667*t128*t30*t40*
        t47+2.0*t128*t29*t38*t47+2.0*t219*t29*t30*t47+t352+t351+t350+
        t349+t348+t346+t345+t344+t343+t342+t341+t340+t339+t338+t310+
        t309+t29*t38*t6*t251+t29*t24*t30*t251-0.458*t3*t7*t250+1.068666666666667*
        t2*t27*t7*t250+0.152666666666667*t2*t3*t24*t25*t250+0.916*
        t3*t7*t218+0.305333333333333*t2*t3*t47*t25*t218-7.124444444444446*
        gradb*t106*t7*t20+0.814222222222222*gradb*t2*t27*t84*t47*t68*
        t20-0.814222222222222*gradb*t27*t24*t47*t68*t20+0.407111111111111*
        gradb*t2*t27*t118*t24*t68*t20-0.949925925925926*gradb*t2*t53*
        t117*t68*t20+0.407111111111111*gradb*t27*t117*t68*t20+0.610666666666667*
        gradb*t27*t84*t25*t20-4.74962962962963*gradb*t2*t106*t47*t25*
        t20-1.424888888888889*gradb*t53*t24*t25*t20-2.374814814814815*
        gradb*t2*t106*t24*t25*t20-0.305333333333333*gradb*t2*t27*t211*
        t25*t20+0.712444444444445*gradb*t2*t53*t118*t25*t20-0.305333333333333*
        gradb*t27*t118*t25*t20-0.949925925925926*gradb*t2*t27*t24*
        t117*t146*t20-7.124444444444446*t2*t53*t7*t127-0.407111111111111*
        t2*t3*t24*t47*t68*t127+0.305333333333333*t2*t3*t84*t25*t127+
        0.305333333333333*t3*t47*t25*t127-0.712444444444445*t2*t27*
        t47*t25*t127-0.305333333333333*t3*t24*t25*t127+t29*t122*t6*
        t103+t29*t118*t30*t103-3.562222222222223*t2*t53*t7*t101-2.137333333333334*
        t27*t7*t101-0.203555555555556*t2*t3*t117*t68*t101-0.305333333333333*
        t3*t47*t25*t101+0.152666666666667*t2*t3*t118*t25*t101)+t353+
        0.004235*t252;
    d4fdrarbrbab = 0.004235*t1*(-4.0*t3*t47*t6*t87+2.0*t131*
        t29*t6*t87-4.0*t3*t49*t6*t84+2.0*t131*t29*t30*t84+t29*t30*
        t6*(4.0*t118*t8-24.0*t47*t56-12.0*t24*t56-2.0*t211*t5+t363+
        t317)-0.458*t2*t3*(t284+t315)*t7-0.458*t3*t253*t7+1.068666666666667*
        t2*t27*t253*t7-18.11644444444445*t2*t203*t7-7.124444444444446*
        t2*t53*t130*t7-3.562222222222223*t2*t53*t109*t7-2.137333333333334*
        t27*t109*t7+0.814222222222222*t4*t53*t84*t47*t68-0.814222222222222*
        t2*t53*t24*t47*t68-0.407111111111111*t2*t3*t130*t24*t47*t68+
        0.407111111111111*t4*t53*t118*t24*t68+0.407111111111111*t2*
        t53*t117*t68-0.203555555555556*t2*t3*t109*t117*t68-0.949925925925926*
        t4*t106*t117*t68-1.777777777777778*t30*t47*t53*t6-0.888888888888889*
        t24*t30*t53*t6+1.333333333333333*t27*t47*t49*t6-0.666666666666667*
        t131*t40*t49*t6+2.0*t221*t29*t49*t6-0.333333333333333*t254*
        t30*t40*t6-2.0*t118*t3*t38*t6-2.0*t211*t3*t30*t6+0.888888888888889*
        t131*t3*t30*t6+0.444444444444444*t112*t3*t30*t6+0.666666666666667*
        t118*t27*t30*t6-2.0*t122*t24*t3*t6-1.777777777777778*t49*t53*
        t54-0.888888888888889*t38*t53*t54-2.0*t213*t3*t54+0.666666666666667*
        t122*t27*t54+2.0*t112*t29*t47*t49+2.0*t131*t24*t29*t49-0.666666666666667*
        t131*t30*t40*t47+2.0*t131*t29*t38*t47+2.0*t221*t29*t30*t47+
        t368+t367+t366+t365+t364+t362+t361+t360+t359+t358+t357+t356+
        t355+t354+t316+t29*t38*t6*t254+t29*t24*t30*t254+0.610666666666667*
        t2*t53*t84*t25-8.41362962962963*t4*t203*t47*t25+0.305333333333333*
        t3*t130*t47*t25-0.712444444444445*t2*t27*t130*t47*t25-0.305333333333333*
        t3*t109*t47*t25+0.152666666666667*t2*t3*t253*t24*t25-2.374814814814815*
        t4*t203*t24*t25-0.305333333333333*t3*t130*t24*t25-1.424888888888889*
        t2*t106*t24*t25-0.305333333333333*t4*t53*t211*t25+0.305333333333333*
        t2*t3*t84*t130*t25-0.305333333333333*t2*t53*t118*t25+0.712444444444445*
        t4*t106*t118*t25+0.152666666666667*t2*t3*t118*t109*t25-0.949925925925926*
        t4*t53*t24*t117*t146+t29*t122*t6*t112+t29*t118*t30*t112)+t369+
        0.004235*t255;
    d4fdrarbgaga = 0.004235*t1*(2.0*t3*t32*t6*t87+t29*t30*
        t6*(-8.0*t18*t2*t85-8.0*t18*t45*t8-8.0*t14*t18*t8+4.0*t14*
        t45*t5+t373)+2.0*t3*t30*t32*t84-0.458*t2*t3*(-2.0*t45*t5-2.0*
        t14*t5+t379+t372)*t7-0.666666666666667*t27*t32*t49*t6-0.333333333333333*
        t257*t30*t40*t6-0.333333333333333*t229*t30*t40*t6-0.666666666666667*
        t27*t32*t38*t6+2.0*t24*t3*t32*t49+2.0*t3*t32*t38*t47-0.666666666666667*
        t27*t30*t32*t47+t371+t370-0.666666666666667*t24*t27*t30*t32+
        t29*t38*t6*t257+t29*t24*t30*t257-0.458*t3*t7*t256+1.068666666666667*
        t2*t27*t7*t256+0.152666666666667*t2*t3*t24*t25*t256+t29*t49*
        t6*t229+t29*t47*t30*t229+0.458*t3*t7*t227+1.068666666666667*
        t2*t27*t7*t227+0.152666666666667*t2*t3*t47*t25*t227-0.407111111111111*
        t2*t27*t24*t47*t68*t18+0.305333333333333*t2*t27*t84*t25*t18-
        0.712444444444445*t2*t53*t47*t25*t18+0.305333333333333*t27*
        t47*t25*t18-0.712444444444445*t2*t53*t24*t25*t18-0.305333333333333*
        t27*t24*t25*t18)+0.004235*t258+0.004235*t230;
    d4fdrarbgagb = 0.0;
    d4fdrarbgaab = 0.0;
    d4fdrarbgbgb = 0.004235*t1*(2.0*t3*t34*t6*t87+t29*t30*
        t6*(8.0*t2*t20*t85-8.0*t20*t44*t8-8.0*t17*t20*t8+4.0*t17*t44*
        t5+t377)+2.0*t3*t30*t34*t84-0.458*t2*t3*(2.0*t44*t5+2.0*t17*
        t5+t379+t376)*t7-0.666666666666667*t27*t34*t49*t6-0.333333333333333*
        t260*t30*t40*t6-0.333333333333333*t237*t30*t40*t6-0.666666666666667*
        t27*t34*t38*t6+2.0*t24*t3*t34*t49+2.0*t3*t34*t38*t47-0.666666666666667*
        t27*t30*t34*t47+t375+t374-0.666666666666667*t24*t27*t30*t34+
        t29*t38*t6*t260+t29*t24*t30*t260-0.458*t3*t7*t259+1.068666666666667*
        t2*t27*t7*t259+0.152666666666667*t2*t3*t24*t25*t259+t29*t49*
        t6*t237+t29*t47*t30*t237+0.458*t3*t7*t235+1.068666666666667*
        t2*t27*t7*t235+0.152666666666667*t2*t3*t47*t25*t235+0.407111111111111*
        t2*t27*t24*t47*t68*t20-0.305333333333333*t2*t27*t84*t25*t20+
        0.712444444444445*t2*t53*t47*t25*t20-0.305333333333333*t27*
        t47*t25*t20+0.712444444444445*t2*t53*t24*t25*t20+0.305333333333333*
        t27*t24*t25*t20)+0.004235*t261+0.004235*t238;
    d4fdrarbgbab = 0.0;
    d4fdrarbabab = 0.0;
    d4fdragagaga = 0.0;
    d4fdragagagb = 0.0;
    d4fdragagaab = 0.0;
    d4fdragagbgb = 0.0;
    d4fdragagbab = 0.0;
    d4fdragaabab = 0.0;
    d4fdragbgbgb = 0.0;
    d4fdragbgbab = 0.0;
    d4fdragbabab = 0.0;
    d4fdraababab = 0.0;
    d4fdrbrbrbrb = 0.004235*t1*(t29*t30*t6*(192.0*t47*t85*
        gradab+16.0*t243*t8*gradab-72.0*t118*t56*gradab-2.0*t382*t5*
        gradab-192.0*t10*t18*t45*t85-192.0*t15*t20*t44*t85-48.0*t10*
        t114*t45*t8-48.0*t115*t15*t44*t8-16.0*t10*t18*t240*t8-16.0*
        t15*t20*t239*t8+72.0*t115*t15*t20*t56+72.0*t10*t114*t18*t56+
        72.0*t119*t15*t56+72.0*t10*t120*t56+8.0*t10*t240*t45*t5+8.0*
        t15*t239*t44*t5+2.0*t15*t20*t381*t5+2.0*t10*t18*t380*t5+6.0*
        pow(t115,2.0)*t15*t5+6.0*t10*pow(t114,2.0)*t5+t278+t277+t276)-
        0.458*t2*t3*t7*(-192.0*t85*gradab-12.0*t115*t15*t8+12.0*t10*
        t114*t8-24.0*t10*t45*t56+24.0*t15*t44*t56-4.0*t10*t240*t5+
        4.0*t15*t239*t5-1.0*t11*t15*t381+t10*t11*t380+t270+t269+t262)-
        0.610666666666667*pow(t118,2.0)*t2*t26*t3*t68+5.333333333333333*
        t3*t48*t49*t6-4.0*t121*t40*t49*t6-4.148148148148148*t27*t39*
        t49*t6+4.0*t244*t29*t49*t6-4.0*t122*t40*t48*t6-4.148148148148148*
        t27*t30*t48*t6+4.0*t245*t29*t48*t6-1.333333333333333*t245*
        t39*t40*t6-1.333333333333333*t244*t30*t40*t6+2.666666666666667*
        t122*t3*t39*t6+2.666666666666667*t121*t3*t30*t6+6.0*t121*t122*
        t29*t6-8.0*t40*t47*t48*t49+12.0*t118*t29*t48*t49+5.333333333333333*
        t3*t39*t47*t49+12.0*t121*t29*t47*t49-4.0*t118*t39*t40*t49+
        4.0*t243*t29*t39*t49+5.333333333333333*t3*t30*t47*t48+12.0*
        t122*t29*t47*t48-4.0*t118*t30*t40*t48+4.0*t243*t29*t30*t48-
        1.58320987654321*t2*t26*t271*t3*pow(t47,4.0)-4.0*t122*t39*
        t40*t47-4.0*t121*t30*t40*t47-4.148148148148148*t27*t30*t39*
        t47+4.0*t245*t29*t39*t47+4.0*t244*t29*t30*t47+42.74666666666668*
        t53*t7*t46+61.74518518518519*t2*t106*t7*t46-2.442666666666667*
        t2*t3*t118*t47*t68*t46+2.442666666666667*t3*t117*t68*t46+5.699555555555556*
        t2*t27*t117*t68*t46+14.24888888888889*t2*t53*t47*t25*t46+8.549333333333335*
        t27*t47*t25*t46+0.610666666666667*t2*t3*t243*t25*t46-1.832*
        t3*t118*t25*t46-4.274666666666668*t2*t27*t118*t25*t46+1.899851851851852*
        t2*t3*t242*t146*t46-1.333333333333333*t243*t30*t39*t40+t29*
        (t280+19.2*t2*t85+t279)*t6*t39+t29*t382*t30*t39+2.666666666666667*
        t118*t3*t30*t39+6.0*t118*t122*t29*t39+6.0*t118*t121*t29*t30+
        t281+t275-61.74518518518519*t106*t7*t26-0.814222222222222*
        t2*t3*t243*t47*t68*t26+2.442666666666667*t3*t118*t47*t68*t26+
        5.699555555555556*t2*t27*t118*t47*t68*t26-9.499259259259262*
        t2*t53*t117*t68*t26-5.699555555555556*t27*t117*t68*t26-14.24888888888889*
        t53*t47*t25*t26-20.58172839506173*t2*t106*t47*t25*t26+0.152666666666667*
        t2*t3*t382*t25*t26-0.610666666666667*t3*t243*t25*t26-1.424888888888889*
        t2*t27*t243*t25*t26+7.124444444444446*t2*t53*t118*t25*t26+
        4.274666666666668*t27*t118*t25*t26-1.899851851851852*t3*t242*
        t146*t26-4.432987654320988*t2*t27*t242*t146*t26+2.849777777777778*
        t2*t3*t118*t117*t146*t26+1.832*t3*t7*t241+4.274666666666668*
        t2*t27*t7*t241+0.610666666666667*t2*t3*t47*t25*t241-21.37333333333334*
        t2*t53*t7*t116-12.824*t27*t7*t116-1.221333333333333*t2*t3*
        t117*t68*t116-1.832*t3*t47*t25*t116-4.274666666666667*t2*t27*
        t47*t25*t116+0.916*t2*t3*t118*t25*t116)+0.01694*t246;
    d4fdrbrbrbga = 0.004235*t1*(-0.458*t2*t3*t7*(12.0*t45*
        t8*grada-6.0*t114*t5*grada+2.0*t11*t240*grada+t307)+t29*t30*
        t6*(-24.0*t114*t18*t8*grada-24.0*t120*t8*grada+72.0*t18*t45*
        t56*grada+12.0*t114*t45*t5*grada+4.0*t18*t240*t5*grada+t308)+
        2.666666666666667*t32*t49*t53*t6*grada+2.0*t245*t3*t32*t6*
        grada-2.0*t122*t27*t32*t6*grada+2.666666666666667*t30*t32*
        t47*t53*grada-4.0*t27*t32*t47*t49*grada+6.0*t118*t3*t32*t49*
        grada+6.0*t122*t3*t32*t47*grada+2.0*t243*t3*t30*t32*grada-
        2.0*t118*t27*t30*t32*grada-2.0*t125*t40*t49*t6+3.0*t248*t29*
        t49*t6-1.0*t248*t30*t40*t6+1.333333333333333*t125*t3*t30*t6+
        3.0*t122*t125*t29*t6+6.0*t125*t29*t47*t49-2.0*t125*t30*t40*
        t47+3.0*t248*t29*t30*t47+t306+t305+3.0*t118*t125*t29*t30+1.374*
        t3*t7*t247+3.206*t2*t27*t7*t247+0.458*t2*t3*t47*t25*t247+21.37333333333334*
        grada*t106*t7*t18-1.221333333333333*grada*t2*t27*t118*t47*
        t68*t18+2.849777777777778*grada*t2*t53*t117*t68*t18+1.221333333333333*
        grada*t27*t117*t68*t18+4.274666666666668*grada*t53*t47*t25*
        t18+7.124444444444446*grada*t2*t106*t47*t25*t18+0.305333333333333*
        grada*t2*t27*t243*t25*t18-2.137333333333334*grada*t2*t53*t118*
        t25*t18-0.916*grada*t27*t118*t25*t18+0.949925925925926*grada*
        t2*t27*t242*t146*t18-10.68666666666667*t2*t53*t7*t124-6.412*
        t27*t7*t124-0.610666666666667*t2*t3*t117*t68*t124-0.916*t3*
        t47*t25*t124-2.137333333333334*t2*t27*t47*t25*t124+0.458*t2*
        t3*t118*t25*t124)+0.012705*t249;
    d4fdrbrbrbgb = 0.004235*t1*(-0.458*t2*t3*t7*(-12.0*t44*
        t8*gradb+6.0*t115*t5*gradb-2.0*t11*t239*gradb+t311)+t29*t30*
        t6*(-24.0*t115*t20*t8*gradb-24.0*t119*t8*gradb+72.0*t20*t44*
        t56*gradb+12.0*t115*t44*t5*gradb+4.0*t20*t239*t5*gradb+t312)+
        2.666666666666667*t34*t49*t53*t6*gradb+2.0*t245*t3*t34*t6*
        gradb-2.0*t122*t27*t34*t6*gradb+2.666666666666667*t30*t34*
        t47*t53*gradb-4.0*t27*t34*t47*t49*gradb+6.0*t118*t3*t34*t49*
        gradb+6.0*t122*t3*t34*t47*gradb+2.0*t243*t3*t30*t34*gradb-
        2.0*t118*t27*t30*t34*gradb-2.0*t128*t40*t49*t6+3.0*t251*t29*
        t49*t6-1.0*t251*t30*t40*t6+1.333333333333333*t128*t3*t30*t6+
        3.0*t122*t128*t29*t6+6.0*t128*t29*t47*t49-2.0*t128*t30*t40*
        t47+3.0*t251*t29*t30*t47+t310+t309+3.0*t118*t128*t29*t30+1.374*
        t3*t7*t250+3.206*t2*t27*t7*t250+0.458*t2*t3*t47*t25*t250-21.37333333333334*
        gradb*t106*t7*t20+1.221333333333333*gradb*t2*t27*t118*t47*
        t68*t20-2.849777777777778*gradb*t2*t53*t117*t68*t20-1.221333333333333*
        gradb*t27*t117*t68*t20-4.274666666666668*gradb*t53*t47*t25*
        t20-7.124444444444446*gradb*t2*t106*t47*t25*t20-0.305333333333333*
        gradb*t2*t27*t243*t25*t20+2.137333333333334*gradb*t2*t53*t118*
        t25*t20+0.916*gradb*t27*t118*t25*t20-0.949925925925926*gradb*
        t2*t27*t242*t146*t20-10.68666666666667*t2*t53*t7*t127-6.412*
        t27*t7*t127-0.610666666666667*t2*t3*t117*t68*t127-0.916*t3*
        t47*t25*t127-2.137333333333334*t2*t27*t47*t25*t127+0.458*t2*
        t3*t118*t25*t127)+0.012705*t252;
    d4fdrbrbrbab = 0.004235*t1*(t29*t30*t6*(12.0*t118*t8-
        36.0*t47*t56-2.0*t243*t5+t317)-0.458*t2*t3*(36.0*t56+t315)*
        t7+1.374*t3*t253*t7+3.206*t2*t27*t253*t7-21.37333333333334*
        t2*t203*t7-10.68666666666667*t2*t53*t130*t7-6.412*t27*t130*
        t7+1.221333333333333*t4*t53*t118*t47*t68-1.221333333333333*
        t2*t53*t117*t68-0.610666666666667*t2*t3*t130*t117*t68-2.849777777777778*
        t4*t106*t117*t68-2.666666666666667*t30*t47*t53*t6+4.0*t27*
        t47*t49*t6-2.0*t131*t40*t49*t6-6.0*t118*t3*t49*t6+3.0*t254*
        t29*t49*t6-6.0*t122*t3*t47*t6-1.0*t254*t30*t40*t6-2.0*t243*
        t3*t30*t6+1.333333333333333*t131*t3*t30*t6+2.0*t118*t27*t30*
        t6+3.0*t122*t131*t29*t6-2.666666666666667*t49*t53*t54-2.0*
        t245*t3*t54+2.0*t122*t27*t54+6.0*t131*t29*t47*t49-2.0*t131*
        t30*t40*t47+3.0*t254*t29*t30*t47+t316+t314+3.0*t118*t131*t29*
        t30+0.458*t2*t3*t253*t47*t25-7.124444444444446*t4*t203*t47*
        t25-0.916*t3*t130*t47*t25-2.137333333333334*t2*t27*t130*t47*
        t25-4.274666666666668*t2*t106*t47*t25-0.305333333333333*t4*
        t53*t243*t25+0.458*t2*t3*t118*t130*t25+0.916*t2*t53*t118*t25+
        2.137333333333334*t4*t106*t118*t25-0.949925925925926*t4*t53*
        t242*t146)+0.012705*t255;
    d4fdrbrbgaga = 0.004235*t1*(t29*t30*t6*(-16.0*t18*t45*
        t8+4.0*t114*t18*t5+4.0*t120*t5+t373)-0.458*t2*t3*(-4.0*t45*
        t5+t372+2.0*t11*t114)*t7-1.333333333333333*t27*t32*t49*t6+
        2.0*t257*t29*t49*t6-0.666666666666667*t257*t30*t40*t6+2.0*
        t122*t3*t32*t6+4.0*t3*t32*t47*t49-1.333333333333333*t27*t30*
        t32*t47+2.0*t257*t29*t30*t47+t371+t370+2.0*t118*t3*t30*t32+
        0.916*t3*t7*t256+2.137333333333334*t2*t27*t7*t256+0.305333333333333*
        t2*t3*t47*t25*t256-4.274666666666668*t53*t7*t18-0.407111111111111*
        t2*t27*t117*t68*t18-1.424888888888889*t2*t53*t47*t25*t18-0.610666666666667*
        t27*t47*t25*t18+0.305333333333333*t2*t27*t118*t25*t18)+0.00847*
        t258;
    d4fdrbrbgagb = 0.0;
    d4fdrbrbgaab = 0.0;
    d4fdrbrbgbgb = 0.004235*t1*(t29*t30*t6*(-16.0*t20*t44*
        t8+4.0*t115*t20*t5+4.0*t119*t5+t377)-0.458*t2*t3*(4.0*t44*
        t5+t376-2.0*t11*t115)*t7-1.333333333333333*t27*t34*t49*t6+
        2.0*t260*t29*t49*t6-0.666666666666667*t260*t30*t40*t6+2.0*
        t122*t3*t34*t6+4.0*t3*t34*t47*t49-1.333333333333333*t27*t30*
        t34*t47+2.0*t260*t29*t30*t47+t375+t374+2.0*t118*t3*t30*t34+
        0.916*t3*t7*t259+2.137333333333334*t2*t27*t7*t259+0.305333333333333*
        t2*t3*t47*t25*t259+4.274666666666668*t53*t7*t20+0.407111111111111*
        t2*t27*t117*t68*t20+1.424888888888889*t2*t53*t47*t25*t20+0.610666666666667*
        t27*t47*t25*t20-0.305333333333333*t2*t27*t118*t25*t20)+0.00847*
        t261;
    d4fdrbrbgbab = 0.0;
    d4fdrbrbabab = 0.0;
    d4fdrbgagaga = 0.0;
    d4fdrbgagagb = 0.0;
    d4fdrbgagaab = 0.0;
    d4fdrbgagbgb = 0.0;
    d4fdrbgagbab = 0.0;
    d4fdrbgaabab = 0.0;
    d4fdrbgbgbgb = 0.0;
    d4fdrbgbgbab = 0.0;
    d4fdrbgbabab = 0.0;
    d4fdrbababab = 0.0;
    d4fdgagagaga = 0.0;
    d4fdgagagagb = 0.0;
    d4fdgagagaab = 0.0;
    d4fdgagagbgb = 0.0;
    d4fdgagagbab = 0.0;
    d4fdgagaabab = 0.0;
    d4fdgagbgbgb = 0.0;
    d4fdgagbgbab = 0.0;
    d4fdgagbabab = 0.0;
    d4fdgaababab = 0.0;
    d4fdgbgbgbgb = 0.0;
    d4fdgbgbgbab = 0.0;
    d4fdgbgbabab = 0.0;
    d4fdgbababab = 0.0;
    d4fdabababab = 0.0;


    ds->df1000 += factor*dfdra;
    ds->df0100 += factor*dfdrb;
    ds->df0010 += factor*dfdga;
    ds->df0001 += factor*dfdgb;
    ds->df00001+= factor*dfdgab;

    ds->df2000 += factor*d2fdrara;
    ds->df1100 += factor*d2fdrarb;
    ds->df1010 += factor*d2fdraga;
    ds->df1001 += factor*d2fdragb;
    ds->df10001+= factor*d2fdraab;
    ds->df0200 += factor*d2fdrbrb;
    ds->df0110 += factor*d2fdrbga;
    ds->df0101 += factor*d2fdrbgb;
    ds->df01001+= factor*d2fdrbgab;
    ds->df0020 += factor*d2fdgaga;
    ds->df0011 += factor*d2fdgagb;
    ds->df00101+= factor*d2fdgagab;
    ds->df0002 += factor*d2fdgbgb;
    ds->df00011+= factor*d2fdgbgab;
    ds->df00002+= factor*d2fdgabgab;

    ds->df3000 += factor*d3fdrarara;
    ds->df2100  += factor*d3fdrararb;
    ds->df2010  += factor*d3fdraraga;
    ds->df2001  += factor*d3fdraragb;
    ds->df20001 += factor*d3fdraraab;
    ds->df1200  += factor*d3fdrarbrb;
    ds->df1110  += factor*d3fdrarbga;
    ds->df1101  += factor*d3fdrarbgb;
    ds->df11001 += factor*d3fdrarbab;
    ds->df1020  += factor*d3fdragaga;
    ds->df1011  += factor*d3fdragagb;
    ds->df10101 += factor*d3fdragaab;
    ds->df1002  += factor*d3fdragbgb;
    ds->df10011 += factor*d3fdragbab;
    ds->df10002 += factor*d3fdraabab;
    ds->df0300  += factor*d3fdrbrbrb;
    ds->df0210  += factor*d3fdrbrbga;
    ds->df0201  += factor*d3fdrbrbgb;
    ds->df02001 += factor*d3fdrbrbab;
    ds->df0120  += factor*d3fdrbgaga;
    ds->df0111  += factor*d3fdrbgagb;
    ds->df01101 += factor*d3fdrbgaab;
    ds->df0102  += factor*d3fdrbgbgb;
    ds->df01011 += factor*d3fdrbgbab;
    ds->df01002 += factor*d3fdrbabab;
    ds->df0030  += factor*d3fdgagaga;
    ds->df0021  += factor*d3fdgagagb;
    ds->df00201 += factor*d3fdgagaab;
    ds->df0012  += factor*d3fdgagbgb;
    ds->df00111 += factor*d3fdgagbab;
    ds->df00102 += factor*d3fdgaabab;
    ds->df0003  += factor*d3fdgbgbgb;
    ds->df00021 += factor*d3fdgbgbab;
    ds->df00012 += factor*d3fdgbabab;
    ds->df00003 += factor*d3fdababab;

    ds->df4000  += factor*d4fdrararara;
    ds->df3100  += factor*d4fdrarararb;
    ds->df3010  += factor*d4fdrararaga;
    ds->df3001  += factor*d4fdrararagb;
    ds->df30001 += factor*d4fdrararaab;
    ds->df2200  += factor*d4fdrararbrb;
    ds->df2110  += factor*d4fdrararbga;
    ds->df2101  += factor*d4fdrararbgb;
    ds->df21001 += factor*d4fdrararbab;
    ds->df2020  += factor*d4fdraragaga;
    ds->df2011  += factor*d4fdraragagb;
    ds->df20101 += factor*d4fdraragaab;
    ds->df2002  += factor*d4fdraragbgb;
    ds->df20011 += factor*d4fdraragbab;
    ds->df20002 += factor*d4fdraraabab;
    ds->df1300  += factor*d4fdrarbrbrb;
    ds->df1210  += factor*d4fdrarbrbga;
    ds->df1201  += factor*d4fdrarbrbgb;
    ds->df12001 += factor*d4fdrarbrbab;
    ds->df1120  += factor*d4fdrarbgaga;
    ds->df1111  += factor*d4fdrarbgagb;
    ds->df11101 += factor*d4fdrarbgaab;
    ds->df1102  += factor*d4fdrarbgbgb;
    ds->df11011 += factor*d4fdrarbgbab;
    ds->df11002 += factor*d4fdrarbabab;
    ds->df1030  += factor*d4fdragagaga;
    ds->df1021  += factor*d4fdragagagb;
    ds->df10201 += factor*d4fdragagaab;
    ds->df1012  += factor*d4fdragagbgb;
    ds->df10111 += factor*d4fdragagbab;
    ds->df10102 += factor*d4fdragaabab;
    ds->df1003  += factor*d4fdragbgbgb;
    ds->df10021 += factor*d4fdragbgbab;
    ds->df10012 += factor*d4fdragbabab;
    ds->df10003 += factor*d4fdraababab;
    ds->df0400  += factor*d4fdrbrbrbrb;
    ds->df0310  += factor*d4fdrbrbrbga;
    ds->df0301  += factor*d4fdrbrbrbgb;
    ds->df03001 += factor*d4fdrbrbrbab;
    ds->df0220  += factor*d4fdrbrbgaga;
    ds->df0211  += factor*d4fdrbrbgagb;
    ds->df02101 += factor*d4fdrbrbgaab;
    ds->df0202  += factor*d4fdrbrbgbgb;
    ds->df02011 += factor*d4fdrbrbgbab;
    ds->df02002 += factor*d4fdrbrbabab;
    ds->df0130  += factor*d4fdrbgagaga;
    ds->df0121  += factor*d4fdrbgagagb;
    ds->df01201 += factor*d4fdrbgagaab;
    ds->df0112  += factor*d4fdrbgagbgb;
    ds->df01111 += factor*d4fdrbgagbab;
    ds->df01102 += factor*d4fdrbgaabab;
    ds->df0103  += factor*d4fdrbgbgbgb;
    ds->df01021 += factor*d4fdrbgbgbab;
    ds->df01012 += factor*d4fdrbgbabab;
    ds->df01003 += factor*d4fdrbababab;
    ds->df0040  += factor*d4fdgagagaga;
    ds->df0031  += factor*d4fdgagagagb;
    ds->df00301 += factor*d4fdgagagaab;
    ds->df0022  += factor*d4fdgagagbgb;
    ds->df00211 += factor*d4fdgagagbab;
    ds->df00202 += factor*d4fdgagaabab;
    ds->df0013  += factor*d4fdgagbgbgb;
    ds->df00121 += factor*d4fdgagbgbab;
    ds->df00112 += factor*d4fdgagbabab;
    ds->df00103 += factor*d4fdgaababab;
    ds->df0004  += factor*d4fdgbgbgbgb;
    ds->df00031 += factor*d4fdgbgbgbab;
    ds->df00022 += factor*d4fdgbgbabab;
    ds->df00013 += factor*d4fdgbababab;
    ds->df00004 += factor*d4fdabababab;

}
